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The following problems are considered: 


I. Paraxial ray tracing by the focal plane method.—For spherical refracting and reflecting 
surfaces it is shown that by drawing two secondary axes, one parallel to a ray incident at an 
arbitrary angle and the other parallel to the path of the same ray after refraction or reflection, 
the positions of the focal planes are readily established and their necessary properties deduced. 
The same is done for a lens by a method based on the fact that a focal plane of the first surface 
is conjugate to a focal plane of the second surface. 

II. The variation in the distance L between the conjugate foci of a thin lens (focal length f) with 
the distance | of one focus from the lens—A general formula is established, and a simple con- 
struction for plotting the graphs between L/f and 1/f is described. 

IIT. Refraction by prism combinations and the prism refractometer—A graphical method is 
used to discuss the design of the Amici prism, the direct vision spectroscope, and achromatic 
prisms giving appreciable deviations. A graphical method is also described for finding the re- 


fractive index of a liquid from results obtained with a prism refractometer. 


I. PARAXIAL RAY TRACING BY THE 
FOCAL PLANE METHOD 


HE method is commonplace; its justifica- 
tion is not. Accordingly the positions and 
properties of focal planes are investigated for 
spherical refracting and reflecting surfaces, and 
the results obtained for the former are applied 
to justify the method for lenses. The diagrams 
are arranged to indicate the correspondence be- 
tween the effect of a converging system on a 
diverging beam and a diverging system on a 
converging beam; and, as usual, distances per- 
pendicular to the principal axes are magnified. 
In each diagram of Figs. 1 and 2 two rays 
[marked (2) ] of a paraxial pencil are shown com- 
ing from the tip of an object which is so distant 
that the rays may be regarded as parallel. One 
of these rays is directed towards the center of 


curvature C of the surface (and is, in conse- 
quence, refracted without deviation or reflected 
back along its own path), while the deviation 
of the other depends on the angle it makes with 
the normal, i.e., the radius, at the point of 
incidence. These two rays finally converge to or 
diverge from R». Two paraxial rays ‘[marked 
(1)], originally diverging from or converging 
towards an object Ri, and subsequently rendered 
parallel by refraction or reflection, are also shown; 
one of these also passes through C. For con- 
venience the ray not associated with C is common 
to both pencils. 

For the refracting surfaces (Fig. 1), if RiQ 
and CR: show the direction’ of two incident 
parallel rays from a distant object, and QR: and 
RC the direction! of two rays rendered parallel 


1 But not the sense for the surface of negative power. 
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(a) Converging surface [Concave to medium of greater 
refractive index (”2,>m)]. Power is positive because a 
parallel beam converges after refraction to form a real 
image in the second focal plane. Power = +72/f2 (by defini- 
tion) =m/f [by Eq. (1)1=(m—m)/Us—fr) (by algebra) 
=(n2—m)/r [by Eq. (2) ]. If the paths of the rays are 
reversed, Ri F; and R2F2 become, respectively, the second 
and first focal planes, the power (=-+m/PF,) being un- 
changed. 
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(b) Diverging surface. (Convex to medium of greater 
refractive index. (m2>m)). Power is negative because a 
parallel beam diverges after refraction to form a virtual 
image in the second focal plane. Power = —n2/fz (by defini- 
tion) = —m/f; [by Eq. (1)]= —(mz—m)/(fe—fi) (by alge- 
bra) = —(m.—m)/r [by Eq. (2)]. If the paths of the rays 
are reversed, RiFi and R2F:; become, respectively, the 
second and first focal planes, the power (= —/PF,) being 
unchanged. 


Fic. 1. Focal planes of refracting surfaces. 


by refraction, then R2CR,Q is a parallelogram 


with CR,=QR,2 and CR.=QR,. It follows that, 
in AQR,C, 


QR:/CRi( = QR:/QR2) = sini,/sint; = n1/N>. 


Moreover, as cosa,;= PF2/QR: and cosa,=PF;/ 
QR, 


PF,/PF2=(QR:/QR:2) (cosa2/cose:) 

= (1/n2)(cosa2/cosa;), 
where RF, and R2F:2 are perpendicular to the 
axis F,PF,. But as a;—0, cosa;—1, and R;—>F;; 
and as a2—0, cosa,—1 and R:—-F», so that F; 
and F; must be the focal points and 


fi/fo=mi/ne. (1) 


Moreover, as CR, and QR: are equal and in- 
clined at a; to the axis, PF,;=CF,, or 


fr—fir=r. (2) 


Thus, the distances of the focal planes from the 
surfaces are calculable from the values of 7, 
No, and fr. 


Again, for paraxial rays, 
b,=CF, tana; =f: tana; and 
bo= CF, tanar=f; tana. (3) 
Thus, the distances of R; and R: from the axis 
are independent of 7; and depend only on a; 
and a, respectively. It follows that all rays (of 


paraxial pencils) parallel to those shown are 
associated with R, or with Re. 

For the reflecting surfaces (Fig. 2), the same 
arguments, simplified by the fact that the angles 
of incidence and of reflection are equal, may be 
applied; they lead to similar conclusions includ- 
ing the results f:=/f, say, fet+fi=r, bi=f tana 
and b:=f tanaz corresponding to Eqs. (1), (2), 
and (3), respectively. 

Thus, to find the subsequent path of any 
paraxial ray incident at a point Q on a spherical 
refracting or reflecting surface, either (A) mark 
the point of intersection R,; of the incident ray 
(produced if necessary) with the first focal plane, 
join R, to C, the center of curvature of the sur- 
face, and then draw a line from Q parallel to the 
secondary axis R,C (or CR,); or (B) draw the 
secondary axis parallel to the incident ray and 
mark its point of intersection R2 with the second 
focal plane; then QR» (or R2Q for a surface of 
negative power) is the required direction. 

The justification for the process of ray tracing 
through a Jens is based on procedures A and B 
for a refracting surface. Inasmuch as a lens has, 
in general, finite thickness and terminates in 
two coaxial spherical refracting surfaces, the 
procedures may be applied either (a) directly to 
each lens surface in turn, or (b) to locate the 
focal planes of the lens and derive a new pro- 
cedure based upon the properties of these planes. 
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(a) Converging (concave) surface. Power is positive as (b) Diverging (convex) surface. Power is negative as a 
a parallel beam converges after reflection to form a real parallel beam diverges after reflection to form a virtual 
image in the second focal plane. Power=+n/f=+2n/r. image in the second focal plane. Power= —n/f=—2n/r. 


Fic. 2. Focal planes of reflecting surfaces. 


Fic. 3. Ray tracing through lens of positive power having refracting surfaces of positive powers P; and P2. 
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Method (a) is more direct when values of the 
refractive indices and radii of curvature are 
known, but method (b) has to be applied when 
only the focal and principal points of the lens 
are specified. The discussion of method (b) is 
based on method (a) (which does not, therefore, 
require a separate treatment) and includes the 
derivation of the usual ‘‘thick lens’ formulas. 
This discussion is limited to a lens whose power 
and that of its refracting surfaces are all positive. 
The limitation seems justified as the results 
obtained for this one system are shown to apply 
to the simpler systems normally encountered, 
and the ‘‘thick lens” formulas are better derived 
by an analytical method without any assumption 
concerning the signs of the powers. 

In Fig. 3, the surfaces (radii of curvature r, 
r’; centers of curvature C, C’) separate media 
having refractive indices n, mo, and n’ such that 
ny>n'>n. The positions of the first focal point 
F, of the first surface, and the second focal point 
F,’ of the second surface are known by inserting 
the values of , mo, n’, r, and 7’ in the formulas 
appended to Fig. 1; and it is required to trace 
through the system a ray R,K which passes 
through R, at a height b; above F;. Accordingly 
KL is drawn parallel to RiC and L is joined to 
Re, where the line C’R2, which is parallel to KL, 
cuts the focal plane F.2’R.. Remembering that 
all rays from R, are rendered parallel to each 
other by the first refraction and therefore di- 
rected towards R. by the second refraction, the 
first and second focal planes of the lens are 
located at F, and F2, respectively, by rays which 
also determine the positions of the principal 
planes containing H, and A. Finally, from R, 
is drawn a ray RN; such that, on emerging from 
the lens, it proceeds in a parallel direction N2Ro. 

It is required to derive expressions for a, B, 
fi, fo, HiNi, and H2N> (see Fig. 3) and establish 
the proposition that LR», passes through M (the 
point on the second focal plane of the lens such 
that N2M is parallel to R,K) and is parallel to 
R.N,, where R, is the point of intersection of 
the ray with the first focal plane of the lens. 

As the triangle having sides \ and a (subse- 
quently called Ada) is similar to AA’BR,, and 
using a similar notation, Ad is similar to 


AAB’Rz, 
(b1+b2)/A= (fita)/a=(fe'+8)/B, 
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so that, having regard to the legend of Fig. 1, 
a/B = fi/fe' =(n/P1)(P2/n’). (4) 


By construction R,C, C’R:, and KL are parallel, 
and therefore, 


bi/(r+f1) =b2/(r'+ fo’) =d/t. 


These results, taken in conjunction with that 
contained in Eq. (2), and additional relations 
derived from the similar triangles \a and F,H,B, 
and 48 and F.H.A show that 


a= fi(d/b2) =fit/(r' + fo’) _ ; 
=fit/fi' =fiP2t/no 
and 
B= fa(d/bi) =fat/(r+ fi) : 
= fot/f2= f2Pit/no, 
where f;' and f2 are the (unmarked on Fig. 3) 
first focal length of the second surface and the 
second focal length of the first surface. Hence, 
a/8 = (P2/P1) (f:/fe), and, by comparison with 
Eq. (4), fi/fe=n/n’. Thus, 
n/fi=n'/fr=P, say, ; 
(the power of the lens, since f; and fe 
are, by convention, the focal lengths) 
a/n = (t/no)(P2/P), and 
B/n’ = (t/no)(P1/P). 
Again, putting H,N,=x and H.N2=y and 
considering the similar triangles NiR,F, and 


N2R2F*’, F.H\B and F\F,R,, and F.H.A and 
P.F2' Ro, 


b, fitate fita-fi_ 
bs fo’ +B—y fi 


whence 


(5) 


(6) 


fe 


a ’ 3? (7) 
fx’ +B—fs 


br fitat+(hr—f) 


be fr’ +B—Ga—fy) 


Thus, Eqs. (7) and (8) are satisfied by 


x= y= (fe—f;) =constant. (9) 


This indicates that N,; and N, are fixed points 
independent of the direction which the incident 
light makes with the axis; thus, e.g., a ray RiN; 
is refracted in a parallel direction N2R2 through 
No, so that N; and Ng, like the centers of curva- 
ture of refracting and reflecting surfaces, are 
nodal points. 
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Then, substituting fi=n/P1, fe’=n'/P:2, and 
the values of a, B, fi, and fe contained in Eq. 
(6) in the last two values of b,/be given in Eq. 
(7), readily leads to 


P=P,4+P2—(t/no)P1P2. (10) 


(Equations (10), (6), and (9) fix the positions of 
the cardinal points of the lens in terms of the 
refractive indices, ¢, and P,; and Pe: (or, more 
basically, r and r’); this permits ray tracing by 
method (b) when the cardinal points have not 
been located experimentally.) 
Again, by similar triangles, 
X+Y Xb Y 
b, BBs be 
fi (fit+x) B—x : 
(fita—fi) (f'+B—-x) fe't+B—x | 


since, by Eqs. (4), (5), (7), and (9), a/B= fi/ fa’ 
= (fi/f2)(b1/b2) = (frta—fi)/fe = (fita—fr)/ 
(fitx), and, therefore, Sihitx)/(fita—fi) 
= f.'. It follows that (b2—62) =b2—(X+ Y), and 
LRz is parallel to N2Rz (and also N,R,). 

The unit magnification property of the prin- 
cipal planes is established by using ray Ri:KLR:2 
and equating values for the gradients of RiKQ 
and of Q’LR:2 (or RiN;), whence 


_ a(b1—bi) Bb, 
fitea—-fi fitx 
fird(bi— 1) ba fad, . 
—— bafads he 


by Eqs. (5), (7), and (9). Thus, if the left-hand 
side of the drawing were moved towards the 
right so that H2 coincided with Hi, Q’ would 
coincide with Q and R,N, and N2R, would be- 
come collinear. Figure RN: MQ’ would become 
a parallelogram, since R,N; is parallel to Q’M, 
and by Eq. (9), the projections of these lines on 
the axis are both fe. It follows that N2M is 
parallel to R,K, the proposition is proved, and 
the path of the ray may be traced by utilizing 
the properties of the first and second focal planes 
of the lens and the nodal points. 

If the lens is thin (i.e., t=a=B=0, H, and Hp 
coincide in “‘the plane of the lens” at H, and N, 
coincides with N2 at N, Fig. 4), the path of the 


n+p 


’ 


Fic. 4. Ray tracing through thin lenses. n¥n’. (a) Lens 
with both refracting surfaces of positive power. mo >n'>n. 
(b) Lens with both refracting surfaces of negative power. 
n>n'>no. In the above diagrams the rays associated with 
R, are rendered parallel to RiC and R:2C’ by the first re- 
fraction and converge to or diverge from Rez when re- 
fracted at the second surface, i.e., Ri and Re are conjugate 


foci. (c) m>n’>n, fi>fe. Example of nodal point on 


opposite side of H to material of greater refractive index, 
but, as always, on same side of H as longer focal length. 


ray after refraction is known either by drawing 
a line parallel to R:N (or NR,) or py joining 
KM (or producing MK), where NM is parallel 
to the incident ray. Distance NH is still (f2—f;), 
as may be observed by substituting a=@=0 in 
Eqs. (7) and (8), but it should be observed that 
NH is an essentially positive quantity and equal 
to (f:—fe) if fid>fe. Thus, N (Fig. 4) and Ni 
and N~, (Fig. 3) are situated on the same side of 
H, H,, and Hz, respectively, as the longer focal 
length, and not necessarily on the same side as 
the material of greater refractive index n’. 

If, in addition to the lens being thin, n=n’, 
then f2=fi by Eq. (6), and NH=0, i.e., N co- 
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(a) Lens of positive power. 


IRONS 


(b) 


(b) Lens of negative power. 


Fic. 5. Ray tracing through thin lenses. n=n’. Similarly marked lines are parallel. 


incides with H and the drawing of an undeviated 
ray through the center of the lens is justified. 
Ray tracing through a lens may then be effected 
by one of the relatively simple processes illus- 
trated in Fig. 5. It is sometimes sought to justify 
these constructions by stating that rays directed 
to or from a point R; on a focal plane are parallel 
either before or after refraction. Strictly, this 
should be preceded by two demonstrations: one 
showing that the property of the principal focus 
for paraxial rays is shared by all points on the 
focal “plane,” and the other that JN is situated 
at the lens “‘center.”’ In this last connection it 


may be remarked that the ‘optical center of a 
lens’? O seems to have a somewhat limited 
application in ray tracing, since either the lens 
is ‘thin’? and O coincides with N,; and Ne, or 
it is “‘thick” and more direct use can be made of 
the separated nodal points. 

To trace a ray through two separated thin 
lenses, whether surrounded by the same material 
or not, it is usually simpler to use one focal plane 
of each lens, rather than to calculate the positions 
of the cardinal points of the system by sub- 
stituting the powers of the lenses for P; and P2 
in Eqs. (6) and (10). 


TABLE I. Refraction of converging and diverging beams by a thin lens. 


Converging (+-ve power) lens 


(/f>1) (Uf <1) 
b + - 


L/f y— Ys iH Yi— Ye 


Diverging (—ve power) lens 


+= 474-3 
L=I'-1 
jee 
KL+h fF 

a 

(/f>1) (1/f <1) 

_ + 
Yi— Ys 


—N+ys yi-— Yeo 





2 The point through which pass rays that are undeviated although laterally displaced, and which divides the line 
joining the centers of curvature of the lens, faces either internally or externally in the ratio of the radii of curvature. 
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Il. THE VARIATION IN THE DISTANCE BETWEEN 
THE CONJUGATE FOCI OF A THIN LENS 
WITH THE POSITION OF ONE OF THEM 


Table I is concerned with the refraction of 
converging and diverging beams incident on 
converging and diverging lenses of focal length 
f. Putting L for the distance between the initial 
and final foci, the relations shown under the 
several diagrams are summarized in the equation 


bL 1 


a(—-tba) ff 


in which the positive or negative sign is taken 
according to the power of the lens, a= +1 ac- 
cording to whether the light initially converges 
or diverges, and b=-+1 according to whether 
the light finally converges or diverges. This 
equation may be rewritten 


bL I 1 
+ 


f (1/U/f}+e 


which is a form convenient for investigating 
graphically the relation between the dimension- 
less and essentially positive ratios L/f and 1/f. 
Thus, on Fig: 6 are plotted yi=//f, ye=[{1/ 
(/f)}+1}", and ys=({1/(//f)}-—1]", the cor- 
responding coordinates of the last two curves 
being found from the simultaneous values of the 
intercepts on the axes of a series of straight lines 


6 


Fic. 6. Graph of distance between conjugate foci/focal 
length against distance of one focus from lens/focal length 
for any thin lens. 


TO OPTICAL PROBLEMS 


Fic. 7. The design of an Amici prism. 


passing through the points (—1, 1) and (1, —1), 
respectively. [This is justified since, by similar 
triangles, y2/1=(1/f)/{(//f)+1} and y3/1=(I/ 
f)/€(/f)—1} for 1/f>1 or y3/(yst+1) =(1/f)/1 
for 1/f<1.] The curves representing y, ye, and 
ys have been used to construct others corre- 
sponding to L/f=yi—ye, L/f=yi—ys (for 1/f 
>1) and L/f=—y,+43 (for 1/f<1). These last 
curves demonstrate the relations sought, and it 
may be observed that a minimum value of 
L (=4f=21=2I’) occurs with both types of lens. 


Ill. REFRACTION BY PRISM COMBINATIONS 
AND THE PRISM REFRACTOMETER 


The treatment in this section is based on a 
geometrical construction previously used to dis- 
cuss refraction by a single prism.* This con- 
struction involves the drawing of arcs whose 
radii represent refractive indices, but,. to clarify 
the diagrams, the ratios of these radii are not 
typical of the ratios of the indices. The necessity 
for this is very apparent with the diagrams of 
the direct vision spectroscope and achromatic 
prism, where, because of the small change in 
refractive index with color, arcs would otherwise 
be crowded together. Moreover, on account of 
this small variation, information derived from 
large scale drawings of Figs. 8 and 9 is un- 
likely to be very accurate, although it may well 
economize time in the initial stages of a com- 


3E. J. Irons, Phil. Mag. 34, 608 (1943). 
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i ae 


- 
= 
=- 
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putation which proceeds by 
proximations. 

(A) The design of the Amici Prism and ray 
tracing through a direct vision spectroscope.*—In 
the Amici prism arrangement (Fig. 7a) the 
angles A and B are such that light, having re- 
fractive indices mo, m1, and m2 in the media con- 
cerned, emerges parallel to its incident direction 
and perpendicular to the face of the second prism. 
To determine A in terms of B (supposed known), 
a line ORP, representing the incident ray (Fig. 
7b) is drawn from O to cut the arcs of center O 
and radii m) and m2in Rand P,and 2B= ZOPQ 
is set off on OP so that PQ (which represents the 
direction of the normal at the interface between 
the two prisms) cuts the arc of center O and 
radius m; in Q. Then ZROS= ZA, where RS, 
like OT, is perpendicular to PQ. The proof is as 
follows: In AOPQ, OP/OQ=n2/n,=sin0QS/ 
sinig by construction and by geometry, re- 
spectively, so that ZOQS= Zi,’ by the second 
law of refraction. From Fig. 7a, i9 and 7; must 
satisfy the equations A =7,+7,' =i9+t2=i9 + B 
and mp sinip =, sint;, and, by reference to Fig. 
7b, it may be observed that putting 2QRP 
= Zip and ZRQO= Z1%, fulfills these conditions. 
This shows that ZA= ZRQS; RQ represents 
the normal at the first surface; and, from the 


successive ap- 


4J. H. Dowell, Proc. Optical Convention II, 965-981 
(1926) has considered direct vision prisms by essentially 
the same method as that employed here, except that he 
deals with wave fronts instead of rays. 
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Fic. 8. Ray tracing 
through a direct vis- 


Consrauctiort Livés Arp on ee 


D. LigurT. 
C. Ligur. 
F. LiGur. 


dimensions marked on Fig. 7b, 


(m2—Mo)sinB 


tanA =, 
V (ny? —n-?? sin?B) — ny cosB 


(Note that if Q lies between S and T, the de- 
nominator of this expression is negative, and 7A 
is obtuse as in Fig. 8.) 

Smith® derived the above formula analytically 
and observed that, as light of any color other 
than that considered will suffer a deviation which 
is a function of the corresponding refractive 
indices, the Amici prism constitutes the most 


C how ~—.— 
€ lust ~~~ 


Fic. 9. Design of an achromatic prism producing a 
large deviation. 


5 C. J. Smith, Phil. Mag. 33, 775 (1942), especially p. 804. 





GRAPHICAL SOLUTIONS TO OPTICAL PROBLEMS 


Fic. 10. Graphical determi- 
nation of refractive index by the 
prism refractometer. A =72+72’. 


simple form of direct vision spectroscope. He 
also noted that an arrangement of two identical 
Amici prisms (i.e., three prisms in all, Fig. 8) 
increases the dispersion and eliminates the lateral 
displacement. 

Ray tracing through such a compound prism 
is readily performed as shown in the lower dia- 
gram of Fig. 8. Here are indicated: (a) the con- 
struction lines by which the value of the ZA 
has been derived from a knowledge of 7 B/2 and 
No, MN, and m2 for sodium light; (b) additional full 
lines to represent normals to the third and 
fourth refracting surfaces (i.e., parallels to P’Q’ 
and Q’R’, respectively); and (c) arcs representing 
refractive indices for, say, C and F light. The 
directions of rays corresponding to yellow light 
are omitted, but the successive directions of the 
rays for C and F light are indicated and num- 
bered to correspond with the upper diagram of 
Fig. 8. Thus, for the F light, direction (1) and 
ZPRQ= Zi having already been fixed, OQ 
represents direction (2) since, in the nomen- 
clature of Fig 7a, OQ/sinip=OR/sinOQR =1/ 
sini,, Again, QP’ being perpendicular to the 
second refracting surface and ZOQU= Zi’, 
OP’ /sini,’=OQ/sinOP’Q, so OP’ gives direction 
(3); directions (4) and (5) are determined 
similarly. 

In practice, an odd number of prisms in excess 
of three may be used, but, although the con- 
struction is readily extended, the diagram be- 
comes complicated. 

(B) The achromatic prism is a compound prism 
designed to produce deviation without, as far as 
possible, any dispersion (Fig. 9).—When the re- 
quired deviation (4) is small, the refracting 


angles A and B of the component prisms are 
readily estimated by assuming that the deviation 
is independent of the angle of incidence. Larger 
values of 5 require larger values of A and B and 
a specified angle of incidence; some indication 
of the values of A and B may then be derived 
from a large scale drawing of the type shown in 
the upper diagram of Fig. 9. Here the concentric 
arcs have radii representing the refractive index 
of air and the refractive indices of the two glasses 
for the two colors (say those of the C and F 
lines) for which the prism is to be achromatized. 


Fic. 11. Graphical determination of refractive index by the 
prism refractometer. A =72—7,’. 


Line OG represents the direction of the incident 
light and 6(= 2 GOH) is set off equal to the re- 
quired common angle of deviation. Line GJK 
represents the normal at the point of incidence, 
and OJ and OK the directions of the refracted 
rays in the first prism. It is then required to 
draw through J and K two parallel lines JJ’ 
and KK’ which cut the arcs representing the 
refractive indices of the glass of the second prism 
at J’ and K’ and make H, J’, and K’ collinear; 
this is readily effected by drawing a number of 
trial lines through H and: using a parallel-rule 
to test the parallelism of JJ’ and KK’. Alter- 
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natively, if one of the prism angles (say A) is 
specified, Zi) may be varied, ZA marked off at 
K, and JJ’, KK’, and HJ’K’ located by trial. 
This done, the lower figure is derived from the 
upper. The construction is justified in the same 
way as before, and may be checked for F light 
from the similarly marked angles in the figures. 
(C) The prism refractometer —Figure 10 will 
recall the variation of Wollaston’s method for 
determining the refractive index of a liquid by 
using a prism of refracting angle A. Light enters 
one refracting face of the prism (or enters and 
is scattered from the face opposite the refracting 
edge), and the direction which separates the 
light and dark fields of view is located at 72 (90 
+7,) to the second refracting face. If the values 
of 4; (say, 7, and 7, respectively) are found when 
those of 72’ correspond to the critical angles first 
for glass and air (42), and second for glass and 
liquid (72’), the refractive indices of the glass and 
liquid may be found by the method shown in the 
right-hand diagram of Fig. 10. Here OP (P being 
the point of intersection with the arc of center O 
whose radius represents ;) is the direction of the 
emergent ray at £7, to PQ, a line parallel to 


the normal. The 2 PRS= ZA is set off on PQ, 


IRONS 


and Q is located by drawing parallel to RS the 
tangent QT to the arc of radius m,;. Then, since 
the three equations OQ/sini;=m,/siniz (from 
AOPQ), 0Q/sin90 =1,/sini,’ (from AOQT) and 
A =12+47’ are satisfied, OQ represents m2 in length 
and gives the direction of the ‘‘critical ray” in - 
the glass, ZOQP= Zi2, and ZOQT= Zi’. An 
arc of radius OQ is now constructed, and OP’ is 
drawn to make 7; with P’Q’, which, representing 
a second normal to the surface, is parallel to 
PQ and cuts the arc of radius OQ at Q’. From 
Q’, Q’T’ is drawn parallel to QT to cut OT pro- 
duced in 7’, and, by similar reasoning to that 
given above, OT’ represents ; on the same scale 
as OT represents m;. Figure 11 is similar to Fig. 
10, but the emergent ray is on the opposite side 
of the normal and the construction for deter- 
mining m2 is omitted. From the magnitudes 
marked on the diagrams there results the usual 
formula 


n3=sinAy (n?—n,? sin’2,;) —cosA (m; sini;), 


in which the value of 7, and therefore sini), 
must be prefaced by a + sign according as the 
conditions are those of Fig. 10 or Fig. 11. 


The whole system of Nature is one immense series of causes and effects, whose beginning and end 
are equally hid in the depths of infinity. Only a small, a very small, portion of it comes under our 
immediate observation; being exposed alike to the sight and other senses of all mankind. Almost 
every phaenomenon is, at once, the cause of manifold effects; and one effect, among many, of a 
superior cause. The business of Science is to extend our views, by unfolding the latent causes which 
exist in Nature; and thence explaining their manifest effects. The discovery of one such real cause, 
unknown before, if it be of general or very extensive influence, as that of universal gravity, is to be 
esteemed a great advancement of natural philosophy. To undervalue such a discovery, as some have 
done, because the cause of that cause cannot yet be assigned, is highly absurd: since the same objection 
must for ever ly against all causes, except primary ones; which are certainly removed far beyond 
the reach of human inquiry. The proper office, and highest boast of true philosophy, is, to bring us 
still nearer to the Deity, by leading us upwards, step by step, in the mighty scale of Nature.— 
Tuomas MELVvILL, “Observations on Light and Colours,’’ Essays and Observations, Physical 


and Literary (Philosophical Society, Edinburgh, 1756), Vol. II, pp. 12-90. (Contributed by 
E. C. Watson.) 





Computation of the Moment of Inertia of Various Bodies by the 
Approximate Sum of Numerical Series 


EppiE Ortiz Muwfiz 
College of Agriculture and Mechanic Arts, Mayaguez, Puerto Rico 


(Received March 12, 1952) 


The paper presents a simple method of computing moments of inertia of various bodies by 
approximating the sum of numerical series. The method lends itself to be used in a course of 
college physics where most students are not familiar with calculus. 


N the analysis of problems concerning the 

mechanics of rigid bodies, the expression 
I= mr’, where m is the mass of a particle of 
the body at a distance r from the axis of rotation, 
occurs frequently, notably in determining the 
behavior of a rotating body. We might appro- 
priately name >-mr? the rotational inertia. 

In determining the moment of inertia of a 
body with respect to a line it is always possible 
to select elements of mass so that: 


1. All points in each element of mass are 
equally distant from the axis with respect to 
which the moment of inertia is to be found. 

2. The moment of inertia of the element, with 
respect to the axis about which the moment of 
inertia of the whole body is to be found, is known; 
the moment of inertia of the body is found by 
summing up the moments of inertia of the 
elements. 

3. The centroid of the elements is known and 
also the moment of inertia of the element with 
respect to an axis which passes through the 
centroid of the element and is parallel to the 
given axis; the moment of inertia of the element 
may then be expressed by means of the parallel 
axes theorem, that is 


IT=I+méd, 


where I, I, m, and d are, respectively, the mo- 
ment of inertia of the element about the axis 
under consideration, the moment of inertia with 
respect to the centroid, the mass of the element, 
and the distance between the two parallel axes. 
The moment of inertia is then found by sum- 


ming the moments of inertia of the different 
elements. 


The purpose of this paper is to present a 
simple method of computing the moments of 
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inertia of various bodies by approximating the 
sum of numerical series. It is hoped that the 
method will be found useful to students not 
familiar with calculus. 


NUMERICAL SERIES INVOLVED 


Throughout this paper the body under con- 
sideration will be divided into a large finite 
number of units. As the size of the unit is arbi- 
trary, as small as one pleases, and hence the num- 
ber of units very large, the sum of the necessary 
arithmetical series can be approximated as shown 
below: 


Approxi- 
mate 


Series value 


Exact value 
N(N+1) 
2 
N(N+1)(2N+1) 


6 


i r 


n? 


(2) 


+ 
™ 
3 


N 

N?(N+1)? Nt 
4 4 

: N(N+1)(2N+1)[3N(N+1)—1] N® 
30 5. 


. nt (3) 


N 


n 
1 


(4) 


MOMENT OF INERTIA OF VARIOUS BODIES 


1. Moment of inertia of a triangular plate about 
its base. Suppose the height h of the triangle 
(Fig. 1) is divided into h units; then the tri- 
angle is built up by laying side by side A rec- 
rangles having unit height. 

The moment of inertia of the shaded portion 
n units from the base is AJ=can?, where a 
=b(h—n)/h, and o is the mass per unit area. 
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a 


Fic. 1. Moment of inertia of a triangular plate about 
its base. The mass element is selected according to the 
first rule. 


The moment of inertia of the triangular plate 
about AA is 


ab(h—n)n? ob_ kh h 


—h> wd n* |. 


h 
I=> 
1 h h 1 


Using the approximate value from series (2) 


and (3), 
, ob (- -) obh? 
“RNZ 47° 12° 


As the mass of the triangle is M=obh/2, then 
I= Mh?/6. 

2. Moment of inertia of a rod about an axis 
normal to the length at the center. Suppose the rod 
is divided into LZ units each of mass m (Fig. 2), 
then the moment of inertia of the mass units 
from the axis is AJ =mn?. The moment of inertia 
of the rod about AA is 


L/2 
I=2m > ns’. 
1 


Using the approximate value for series (2), 


Fic. 2. Moment of inertia of a rod about an axis normal to 
the length at the center. 


I=mL*/12. Since the mass of the rod is M=mL, 
then J= ML?/12. 

3. Moment of inertia of a rectangular plate 
about an axis normal to the plate at the center. 
Suppose the length 6 of the rectangular plate is 
divided into 6 units (Fig. 3); then the rectangle 
is built up by laying side by side b rods having 
unit width. As the moment of inertia of a rod 
about an axis normal to the length at the center 
is known, applying the parallel axes theorem, 
the moment of inertia of the shaded rod n units 
from the axis under consideration is Al =oan? 
+oa?/12, where o is the mass per unit area. 


Fic. 3. Moment of inertia of a rectangular plate about 
an axis normal to the plate at the center. The mass element 
is selected according to the third rule. 


The moment of inertia of the rectangular plate 
about AA is 


b/2 a ba* »/2 
I=20a >> n+ —) =200(— + n? }. 
1 12 24 ~ 
Using the approximate value for series (2), 


a? 0b? 


me Aone Bam 
24 24 


oba 


I=2ca — 
12 


(a?+-b?). 


Since the mass of the rectangle is M=cab, then 


M 
I=—(@+2). 
12 


4. Moment of inertia of a disk about an axis 
normal to the plane of the disk at the center. A 
body that has all its mass disposed at equal dis- 





MOMENT OF 


tances from its axis of rotation (for example, 
a circular ring) has a rotational inertia that 
equals its total mass times the square of its 
radius. 

Suppose the disk is divided into R concentric 
rings each having unit thickness (Fig. 4); then 
the moment of inertia of the shaded ring having 
a radius of » units about an axis normal to the 
plane of the ring is AJ=2z0n', where o is the 
mass per unit area. The moment of inertia of 
the disk about AA is 


R 
I=2ne > n'. 
1 


Using the approximate value to series (3), 


I=r7cR*/2. 


Fic. 4. Moment 
of inertia of a disk 
about an axis nor- 
mal to the plane of 
the disk at the cen- 
ter. The element of 
mass is selected ac- 
cording to the first 
rule. 


As the mass of the disk is M=2cR?, therefore, 
I= MR?/2. 


5. Moment of inertia of a right circular cone 
about its geometrical axis. Suppose the height / 
of the cone is divided into h units (Fig. 5); then 
the cone is built up by laying one over the other 
h thin disks having unit thickness. 

The moment of inertia of the shaded thin 
disk about the axis is AJ=pa‘/2, where p is 
the mass per unit volume and a=R(h—n)/h. 
The moment of inertia of the cone about AA is 


wpR* h 
_ ( 
2 1 


4n 6n 4n® n‘ 
1 -=+—-— +"). 
ROE OR aes 


INERTIA OF VARIOUS BODIES 
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_ Fic. 5. Moment of inertia of a right circular cone about 
its geometrical axis. The element of mass is selected accord- 
ing to the second rule. 


Using the approximate values of series (1), (2), 
(3), and (4), 


aphR' 


apR! h 
r= (i—2h+-2h—h4-) = 
2 5 10 


Since the mass of the cone is M=-zphR?/3, then 
I=3MR?/10. 

6. Moment of inertia of a sphere about any 
diameter. Suppose the radius R of the sphere is 
divided into R units (Fig. 6); then the sphere is 
built up by laying one over the other 2R thin 
disks having unit thickness. 

The moment of inertia of the shaded thin 
disk about the axis is AJ = rpr*/2, where p is the 


A 


MO 


Fic. 6. Moment 
of inertia of a sphere 
about any diam- 
eter. The element 
of mass is selected 
according to the 
second rule. 
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mass per unit volume and r=4/(R?—n?). The 
moment of inertia of the sphere about AA is 


R R 
I=np (R?—n?)?= xp (R!—2R'n?-+n!). 
1 1 


Using the approximate values of series (2) 
and (4), 


2R5 
tame R-+ 


R5 


5 


8rpR5 
15 


Since the mass of the sphere is M=42pR?/3, 
then J=2MR?/5. 

It has thus been shown that for a given body 
the moment of inertia can be expressed in terms 
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of the sum of a numerical series. The approxi- 
mate sum of the series being known, the moment 
of inertia can be found at once. It is obvious that 
the method could very well be applied to com- 
putation of center of mass. Because of its sim- 
plicity the method lends itself to be used in a 
course of college physics and in the light of the 
author’s experience it may prove to be a valu- 
able teaching aid. 

It is a pleasure to express my thanks to Pro- 
fessor Braulio Duejfio Leary, Professor Jose 
Medina Hernandez, and Professor Jorge Rod- 
riguez Arias for having given valuable sugges- 
tions and for reading the paper. 


Achievement of Physics Students With and Without Laboratory Work* 


HayM KRuGLAK 
University of Minnesota, Minneapolis 14, Minnesota 


(Received June 9, 1952) 


Two general physics laboratory groups were compared in achievement with a group enrolled 
in the same course without laboratory. The achievement was measured by two pencil-paper 
examinations and two laboratory performance tests. The analysis of variance and covariance 
method was used to hold constant three measures of initial differences among students. There 
were no statistically significant differences between the means of the groups on the mechanics 
theory test. The laboratory groups were significantly better than the no-laboratory group on 


all tests dealing with laboratory work. 


THE PROBLEM 


N experiment on the relative effectiveness of 
the conventional and demonstration labora- 
tory methods was described in a previous paper.! 
The analysis has been extended so as to include 
the data on a group of students who were not 
exposed to laboratory experiences. The present 
investigation was designed to test two null 
hypotheses: (1) There are no differences in the 
learning outcomes of general physics students 
who do laboratory work by the individual method 
and students who take the same course without 
laboratory. (2) There are no differences in the 
learning outcomes of general physics students to 
whom laboratory experiments are demonstrated 
and students who take the same course without 
laboratory. 


* This work was supported n part by the joint program 
of the AEC and ONR 


1H, Kruglak, Am. J. Phys. 20, 136 (1952). 


EXPERIMENTAL PROCEDURE 


The 28 students in the individual or conven- 
tional laboratories and the 28 students in the 
demonstration laboratories were enrolled in 
Physics 1a, during the fall quarter 1950, at the 
University of Minnesota. The characteristics of 
the course population, the method of sampling, 
and the contrasting treatment in the laboratory 
were described in the above-mentioned paper. 

The no-laboratory group consisted of 21 male 
volunteers, enrolled in Physics 1 during the same 
quarter as the two laboratory groups. Physics 1 
and Physics 1a have the same prerequisites, 
cover the same subject matter, and have the 
same lecturer. However, students in Physics la 
take two hours of laboratory per week for one 
additional quarter-credit, while students in 
Physics 1 take no laboratory at all. 

Since randomization is one of the basic as- 
sumptions in the analysis of variance and co- 
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TABLE I. Summary of the instructional outcomes. Physics 1, 1a, fall quarter 1950, University of Minnesota, 


Laboratory 


Max. possible 
method* score 


Criterion 
Mechanics theory C 
D 
N 


132 


Laboratory 
Written 


Short-item laboratory 
performance 


Long-item laboratory 
performance 


® C =conventional or individual, 
D =demonstration, 
N =no laboratory. 


Adjusted 


means> 


58.05 
54.48 
50.80 


Significance of differences between means 
(C—D) Not significant 
(C—N) Not significant 
(D—N) Not significant 
(C+D) —N Not significant 


42.07 
44.20 
30.49 


Not significant 

Significant at the 5 percent level 

Significant at the 1 percent level 
(C-+D)— N Sienificant at the 1 percent level 


90.51 
i234 
56.81 


(C—D) 
(C—N) 
(D—N) Significant at the 
(C+D) —N Significant at the 


Significant at the 1 
Significant at the 1 
1 
1 


percent level 
percent level 
percent level 
percent level 


78.68 
73.62 
52.03 


(C—D) Not significant 
(C—N) Significant at the 
(D—N) Significant at the 


(C+D) —N Significant at the 


percent level 
percent level 
percent level 





_ > Adjusted by analysis of variance and covariance for initial differences in scores on the American Council on Education Psychological Examina- 


tion, mathematics pre-test, and criterion pre-test. 


variance, the sample of volunteers was compared 
to the remainder of the male population in 
Physics 1. No significant differences between the 
means of the volunteers and nonvolunteers were 
found on the ACE, the mathematics pre-test, 
the mechanics theory pre-test, and the labora- 
tory written pre-test. The volunteer group was 
adjudged to be representative of the course 
population. 

The three groups in this study were homo- 
geneous in that every student was a graduate of 
an accredited Minnesota high school, was en- 
rolled in the course for credit, had a score on a 
scholastic aptitude test (ACE), took all the 
special pre-tests and post-tests, and completed 
the quarter’s work in physics. 

The control factors for the three groups were 
(1) identical lectures, (2) identical textbook and 
assignments, (3) identical course quizzes, and 
(4) identical pre-tests and post-tests. 

The students in the conventional groups 
worked in pairs, assembling and operating the 
apparatus as specified in the laboratory manual. 
The laboratory teaching assistants performed all 
the experiments for the demonstration groups, 
without any manipulation on the part of the 
students. The no-laboratory group had no 
laboratory instruction during the quarter. | 


The achievement of each group was measured 
by means of a mechanics theory test, a labora- 
tory paper-pencil test, and two laboratory per- 
formance examinations.! 

The volunteers were given a ten percent in- 
crease in their final grades as compensation for 
taking part in the experiment. The volunteer 
group appeared to be well-motivated and serious- 
minded during the test periods. 


ANALYSIS OF DATA 


The analysis of variance and covariance tech- 
niques were used to test the hypotheses of no 
difference between the means of the groups on 
the four criterion tests. The analysis of covari- 
ance made it possible to hold constafht the stu- 
dent’s scores on the ACE, the mathematics pre- 
test, and the pre-test corresponding to the 
criterion under consideration. The summary of 
the analysis is shown in Table I. 

Testing the null hypotheses by the F test 
gave the following results: 


(a) The hypothesis of no difference between 
the individual laboratory group and the no- 
laboratory group for the 


(1) mechanics theory—accepted, 
(2) laboratory written—rejected, 
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(3) short-item laboratory performance— 
rejected, 

(4) long-item laboratory performance— 
rejected. 


(b) The hypothesis of no difference between 
the demonstration laboratory group and the 
no-laboratory group for the 


(1) mechanics theory—accepted, 

(2) laboratory written—rejected, 

(3) short-item laboratory performance— 
rejected, 

(4) long-item laboratory performance— 
rejected. 


Further inspection of the table shows that on 
three of the tests the adjusted means are in de- 
scending order, with the conventional group at 
the top. On all of the criteria the no-laboratory 
group had the lowest means. 


CONCLUSIONS 


The experimental evidence of this investiga- 
tion supports the conclusion that students who 
get laboratory instruction, by the individual or 
demonstration method, are superior to students 
without such instruction, on tests designed to 
measure laboratory outcomes. 

Since no significant differences between the 
means of any two groups were found for the 
mechanics theory test, it is reasonable to con- 
clude that laboratory instruction does not sig- 
nificantly influence scores on pencil-paper tests 
designed to measure a knowledge of facts, prin- 
ciples, and applications of elementary mechanics. 

The difference between the means of the in- 
dividual and the demonstration groups on the 
short-item performance test is about equal to 
the difference between the means of the demon- 
stration and the no-laboratory groups on the 
same test. This fact might be used as evidence in 
support of the “learning by doing’’ pedagogical 


concept, since the short-item performance test 
attempts to evaluate simple manipulatory skills, 
measuring techniques, and apparatus knowledge. 


RECOMMENDATIONS 


On the basis of this investigation it is recom- 
mended that the evaluation of laboratory in- 
struction be made in terms of the objectives of 
laboratory work. 

It is recommended that pencil-paper and per- 
formance tests be used to measure the general 
and specific outcomes of laboratory instruction. 

Since there were no significant differences be- 
tween the means of the individual and the 
demonstration groups on the problem-solving 
(long-item) performance test, it is recommended 
that more time be spent in the laboratory on 
problematic situations. 

The evaluation of laboratory teaching is of 
considerable interest and importance. Although 
a great deal of staff time, apparatus, and building 
space is given over to instructional laboratories, 
there have been relatively few investigations 
dealing with measurable outcomes of laboratory 
instruction at the freshman and sophomore 
levels. It is recommended and hoped that similar 
studies will be carried out in other schools, with 
nonscience majors as well as technical students, 
with other methods and tests, and for longer in- 
structional periods. 

The writer recognizes that all the conclusions 
and recommendations of this study are limited 
by the population used in the three groups, by 
the design of the experiment, and by the four 
measures of instructional outcomes. 

It is a pleasure to acknowledge the invaluable 
help of Mr. W. J. Moonan in planning the sta- 
tistical treatment. Thanks are due to Mrs. N. 
Horwitz for carrying out most of the computa- 
tions. The encouragement and advice of Pro- 
fessor C. N. Wall are greatly appreciated. 


There have been lean seasons for adventurers before, the eighteenth century notably, when every- 
thing seemed owned, done, mapped. In such times the new is to be sought inwards, not in immutable 
Nature, but in the ever renewed flux of human life-—W1LL1AM BoLiTHO, Twelve Against the Gods. 





Determination of Molecular Structure from Microwave Spectroscopic Data 


J. KraltcHMAN* 
Westinghouse Research Laboratories, East Pittsburgh, Pennsylvania 


(Received April 17, 1952) 


A method is described for determining the position of an atom in a molecule from spectro- 
scopic measurements on two isotopic species of the molecule. The method is applied to various 
types of molecules; explicit expressions are derived for linear, symmetric top, planar, and non- 
planar asymmetric top molecules. The number of isotopic species on which measurements 
must be made to complete the structural determination, i.e., determine the position of every 
atom in the molecule, is discussed for various types of molecules. An application of the method 
to the determination of mass difference ratios is also considered. 


I. INTRODUCTION 
ETERMINATION of the size and shape of 


molecules from spectroscopic measurements 
in the microwave region has now become very 
common. Once the major work of identifying 
the spectrum of the molecule has been com- 
pleted, it is possible to calculate rotational con- 
stants for the molecule, and from these the 
molecular moments of inertia. Except for di- 
atomic molecules and a few asymmetric top 
molecules, e.g., NH2D, SOs, etc., the moments 
of inertia of one isotopic species of a molecule 
are not sufficient to determine the molecular 
configuration. However, when sufficient in- 
formation is available on the moments of in- 
ertia of other isotopic species of the molecule, 
it becomes possible to determine the molecular 
configuration completely. In this paper a general 
and quite straightforward method of using this 
information to determine structural parameters 
is discussed. In particular, it is found that simple 
but exact expressions can be derived for deter- 
mining the position of any atom in a molecule. 
Usually the spectroscopist is satisfied if he can 
determine the equilibrium structural parameters 
of the molecule, i.e., the positions the atoms 
would occupy if they were at rest in the mole- 
cule. However, all spectroscopic measurements 
must be made on vibrating molecules, since even 
in the lowest or. ground state the atoms have a 
certain amount of vibrational energy—the so- 
called zero-point energy. Consider, for example, 
a diatomic molecule in which the distance be- 
tween the atoms is 7. If the vibrations were 


* Present address: Co Nee York of Chemistry, Columbia 
University, New York 27, New York. 
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harmonic, the average distance between the 
atoms during a vibration rx, would be equal to 
the equilibrium distance r,. Thus, the equi- 
librium moment of inertia would be proportional 
to r2 or to (r)?. However, the average moment 
of inertia for the vibration is proportional to 
the average square of the distance between the 
atoms, (r?). Since (77), is not equal to (rw), the 
average moment of inertia is not equal to the 
equilibrium moment of inertia. Hence, even if it 
were possible to determine the average moment 
of inertia, the positions computed from it would 
not correspond to the equilibrium positions. 

Actually, however, the analysis of a rotational 
spectrum gives not the average moments of 
inertia, but rather rotational constants, which 
are essentially reciprocals of the moments of 
inertia averaged over a vibrational state. The 
moments of inertia computed from these rota- 
tional constants are called effective moments of 
inertia. Thus, in a diatomic molecule the rota- 
tional constant is proportional to the average 
value of 1/r?, i.e., to (1/r?)w. Since (1/r?) is not 
equal to 1/(rw)*, he effective moment of inertia 
will not be equal to the equilibrium ‘moment of 
inertia. Consequently, positions computed from 
it will not correspond to equilibrium positions. 
(Since (1/r?)« is not equal to 1/(r?), the effective 
moment of inertia will also differ from the aver- 
age moment of inertia.) 

In some cases the average distance between 
the atoms during a vibration r~ would be of 
greater interest than the equilibrium distance r.. 
As mentioned above, if the vibrations were 
harmonic, these two quantities would be equal. 
However, to complicate matters further the 
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vibrations are actually anharmonic, so that ra, 
is not equal to r,. Since it usually has been im- 
possible even to calculate equilibrium param- 
eters, little or no attention has been devoted to 
the problem of determining average distances. 

In the theory it is always possible to obtain 
equilibrium parameters. By determining the ro- 
tational constants of the molecule in a number 
of different vibrational states, it is possible to 
compute the equilibrium rotational constants.! 
From these the equilibrium moments of inertia 
and the equilibrium positions of the atoms can 
be determined. However, except for a few simple 
molecules, it has been impossible to measure the 
rotational constants in a sufficient number of 
different vibrational states to make this possible. 
In fact, in the vast majority of cases the micro- 
wave spectroscopist has been only able to deter- 
mine the rotational constants in the ground state. 

Of course, the spectroscopist could resign 
himself to using just effective moments of in- 
ertia. However, when it becomes necessary to 
use isotopic species of a molecule to complete 
the structural determination, the implicit as- 
sumption is made that the structural parameters 
are the same in all isotopic species of the mole- 
cule. Actually this assumption is only valid for 
the equilibrium structural parameters. Hence, if 
only to obtain a set of self-consistent parameters, 
equilibrium moments must be used. If for lack 
of sufficient information the spectroscopist is 
forced to use effective moments, not only will 
the calculated structural parameters differ from 
the equilibrium parameters, but they may not 
even exactly fit all the isotopic species measured. 

These factors, which are responsible for the 
differences between calculated and equilibrium 
parameters, are known as the zero-point vibra- 
tion effects; it is these zero-point effects that are 
responsible for the largest uncertainties in struc- 
tural parameters determined from microwave 
data. 


Il. PRELIMINARY REMARKS 


The inertia dyadic can be defined as? 


I= > m,(r?1—rzi), (1) 


1For details, see G. Herzberg, Infrared and Raman 
Spectra of Polyatomic Molecules (D. Van Nostrand Com- 
pany, Inc., New York, 1945). 

2See, for example, H. Goldstein, Classical Mechanics 
(Addison-Wesley Press Inc., Cambridge, 1950), Chap. 5. 
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where m; is the mass of the ith atom, r; the dis- 
tance of the ith atom from the origin, and 1 the 
unit dyadic. The diagonal elements of I are 
known as the axial moments of inertia or, 
simply, the moments of inertia, a typical one 
having the form 


I22= DLim(y2?+z2?). (2a) 


The nondiagonal elements are designated prod- 
ucts of inertiat and have the form illustrated by 


(2b) 


In matrix notation the inertia dyadic can be 
written 


Iy=Iyz2= — DV imxyi. 


Tes Iz: 
[I= Ta Ey, 
ae } 


It is always possible to rotate the coordinate 
axes so that this matrix is diagonal. In the 
language of matrix algebra such a diagonaliza- 
tion can be obtained by means of a similarity 
transformation. The rotated axes are called the 
principal axes, and the corresponding diagonal 
elements, Ix, Iy, Iz, are known as the principal 
moments of inertia. In practice the principal 
moments of inertia, being eigenvalues of the 
inertia matrix, are obtained as the three roots 
of the secular equation 


I.,—I lw as 
les Iy—I Tye 
| Tus I,,—I 


=0. 


Spectroscopic measurements only give in- 
formation as to the values of the principal mo- 
ments of inertia with respect to axes passing 
through the center of mass of the molecule. If 
the origin of the coordinate axes is chosen at the 
center of mass of the molecule and Eq. (2) is 
used to compute the elements of the inertia 
dyadic, the diagonalization procedure will auto- 
matically give the principal moments of inertia 
with respect to the center of mass. However, it 
is often desirable to start with an arbitrary 
origin for the coordinate system and still make 

3 Certain very small errors are introduced by assuming 
that the nucleus of an atom and its associated electrons 
are located at a point. However, these errors are negligible 


compared to the precision to which the moments of inertia 
can be determined. 


4 Often their negatives are called the products of inertia. 
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the computations with respect to the center of 
mass. This can easily be done by making use of 
the parallel axis theorem. Let r; represent the 
distance of a particular atom from the arbitrary 
origin, r;° the distance of this atom from the 
center of mass, and r the distance of the center 
of mass from the origin. Then 


r;=r,°+r, 
and by definition of the center of mass 
== dYmazi/ > mi, Lmr?= 0. 


Upon substitution of these relations into Eq. 
(1), the following expression is obtained : 


Ie= Db m(r?71—rx;) — (1/2 m,) [(Lam as) 
‘(Smazi1—(CLamaz) (Cams), (3) 


where I? is the inertia dyadic with respect to 


the center of mass. Typical elements have the 
form$ 


Te2°= Dimi(y?+2,*) 
, (Leamiyi)? (Li mzi)? 


, 4 
Lim; Lim; ”" 


(Limixi) (Limi) 
Lim; 


Iy°= —Ymxyit (4b) 


With an arbitrary origin for the coordinate axes 
the diagonalization procedure will automatically 
give the principal moments of inertia with re- 
spect to the center of mass if the elements of 
the inertia dyadic are computed from Eq. (4). 
Note, that if the arbitrary origin is the center of 
mass, Eq. (4) reduces to Eq. (2). In the future, 
since all moments of inertia will be computed 
with respect to the center of mass, Eq. (4) will 
always be used; for simplicity the superscript 0 
will be omitted. 

In dealing with asymmetric top molecules it 


will be more convenient to use another symmetric 
dyadic defined by 


(maxi) (Lima) 
Lm; 


where P is the planar dyadic with respect to the 


P= Lmrrti— ’ (5) 


5 The usefulness of these relations in computing moments 
and products of inertia has previously been pointed out 
by J. O. Hirschfelder, J. Chem. Phys. 8, 431 (1940). 


center of mass. If the origin of the coordinate 
axes is at the center of mass, P reduces to 


P=) mray. 
Typical elements of P have the form 
Pre= DLimax?—[(Limxi)?/Lmi), 
P,=Limxy: 
—C(Lamx,)(Lemyi)/Limi], (6b) 


which, when the origin of the coordinates axes 
is at the center of mass, reduce to 


Pi2= Lmx?, 


(6a) 


(7a) 
(7b) 


The diagonal elements of P are often called 
planar moments of inertia. The nondiagonal ele- 
ments, except in choice of sign, are equivalent to 
the products of inertia. 

Relations between the elements of P and I can 
readily be derived from the expressions given 
above. In this paper only the relations between 
the principal moments will be needed; the fol- 
lowing are typical of these: 


Px=}([y+Iz—Ix), 
Ixy=Py+Pz. 


Pay= DL mxiyi. 


(8) 
(9) 
Ill. GENERAL METHOD 


Consider a molecule of mass M for which the 
principal moments of inertia have been found to 
be Ix, Iy, Iz. Let a coordinate system be set up 
with the origin at the center of mass of the mole- 
cule and with axes X, Y, Z oriented in the direc- 
tions of the principal moments. Since 


Lmx;=0, Kmy:=0, YLwmez;=0, (10a) 


the elements of I given by Eq. (4) reduce to 
I= —> mxyi=0, 
| —Ymxz2;=0, 
Iy.= —> myz2:=0, 
Ixy =I22= Do mi(y2+2?), 
Ty =Iy= YL am(x2+22), 
Iz=I1.2=Dm(x?+y2). 


The nine Eqs. (10) will be referred to as the 
moment equations of the original molecule. 


(10b) 


(10c) 
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Let an atom of mass m in the original molecule 
be replaced by an isotope of mass m-+Am; de- 
note the coordinates of this atom with respect 
to the center of mass of the original molecule 
by x, y, z. It is then convenient to compute the 
elements of I for the isotopically substituted 
molecule using the coordinate system set up for 
the original molecule. Even if this is done, by 
using Eq. (4), the elements of I are given with 
respect to the center of mass of the substituted 
molecule. It is then found that for the substi- 
tuted molecule 


(Amy)? (Amz)? 
M+Am M-+Am 
= Ix+p(y?+2"), 


Taz! = Ix+Am(y?+2") — 


(11) 


Izy’ = — yxy, etc., 


where p= MAm/M-+Am, and M is the mass of 
the original molecule. From spectroscopic meas- 
urements the principal moments of inertia of 
the substituted molecule, Jy’, Iy’, Iz’, can be 
determined. The problem, then, is to compute 
x, y, in terms of Ix, Iy, Iz and Ix’, Iy’, Iz’. 

In general, the principal axes X’, Y’, Z’ of the 
substituted molecule will be rotated with re- 
spect to the principal axes X, Y, Z of the original 
molecule. However, for many simpler molecules 
at least one of the principal axes of the substi- 
tuted molecule will be parallel to a principal 
axis of the original molecule. This special con- 
dition will be met whenever one of the principal 
axes of the original molecule is perpendicular to 
the line joining the center of mass and the atom 
to be isotopically substituted. If this perpendicu- 
lar principal axis is the Z axis, then the Z’ axis 
will be parallel to the Z axis and the principal 
moments of the original and substituted mole- 
cules will obey the relation 


luo Tg ly! ~ Bet Ty — Te. (12) 


Since to date most investigations in microwave 
spectroscopy have involved molecules which 
meet this condition, such molecules will be con- 
sidered first. The more general case will be 
treated in Sec. VII. 


IV. LINEAR MOLECULES 


Let the Z axis be oriented along the molecular 
axis. Then Jxy=Ty, and Jz=0. With this choice 
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of axes the X’ axis can be chosen parallel to the 
X axis and the Y’ axis parallel to the Y axis; 
the Z’ axis will coincide with the Z axis. Since 
the x and y coordinates of all atoms in the 
molecule are zero, Eq. (11) reduces to 


Tey’ = Tez! =Iy2/ =0, 
Ty’ = Ix+p2’, 

Ty’ =Iy+u2*, 
Iz'=Iz=0. 


Thus, the absolute value of the distance of the 
substituted atom from the center of mass of the 
original molecule is 


|2| =[u-(Ly’ — Ty) }#=[u'(Ix’ — Ix) J}. 


(13) 


V. SYMMETRIC TOP MOLECULES 


Let the Z axis be oriented along the symmetry 
axis of the molecule. Then Jx = Jy. Consider the 
following types of substitutions: 


(A) Substitution on an Atom Located 
on the Symmetry Axis 


With the choice of axes given above the X’ 
axis can be chosen parallel to the X axis and the 
Y’ axis parallel to the Y axis; the Z’ axis will 
coincide with the Z axis. The x and y coordinates 
of any atom located on the symmetry axis are 
zero. In all the mathematical details this case is 
exactly similar to the linear molecule treated 
above, except that Jz and Jz’, while equal, are 
not zero. However, since the microwave spectrum 
of a symmetric top gives no information as to 
the magnitude of Iz or Iz’, this distinction is 
immaterial to the treatment being given here. 
Thus, the distance of a substituted atom, located 
on the symmetry axis, from the center of mass 


of the original molecule can be determined from 
Eq. (13). 


(B) Substitution on an Atom not 
Located on the Symmetry Axis® 


In this case it is always possible to orient 
either the X or Y axis so that one is perpendicular 
to the line joining the center of mass and the 
atom to be substituted. Let the axes be so ori- 
ented that the X axis is this perpendicular axis. 
With this orientation the x coordinate of the 


® While the original molecule is a symmetric top, the 
isotopically substituted molecule is not. 
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substituted atom is zero. Equation (11) becomes 
Iw’ =Iz,'=0, 
Tye’ = — py2, 
Ix’ =Ix+u(y*?+2’), 
Ly’ = Ty+ 2’, 
I,,'= Iz+ny”. 


The third of these equations immediately gives 
y+2?=p-(Ix’ — Ix). (14a) 


Furthermore, the inertia matrix of the substi- 
tuted molecule can be reduced to the form 


-— — pys ) 
Iz+uy 


— yz 
Since the trace of a matrix is invariant under a 
similarity transformation, it is possible to write 


Ty+u2?+Izt+uy?=TIy’+Iz’, 


py+2=p Ty’ —Iyt+Iz'’—Iz). (14b) 


Comparison of this expression with Eq. (14a) 
gives 


Tx’ —Ix =Ty’ —Tyt+ Iz’ —Tz. 


Except for a slight change in notation (since the 
X axis is the perpendicular axis in this case) this 
expression is equivalent to Eq. (12). Further- 
more, since Ixy = Ty, this relation reduces to 

Iz=Iz'+TIy'—Ix’. (15) 
It is also a well-known theorem of matrix algebra 
that the determinant of a matrix is invariant 


under a similarity transformation. This gives the 
equation 


(Iy+u2?)(Izt+ uy?) —w’y’2? =Ty’Tz'. (16) 


Equations (14b) and (16) can now be solved ex- 
plicitly for y and z. The results are 


| =< 
= (Iy —Iz) 
1 
=|-0e'-1)(14 
Le 
4 |- ae) 
13) =| — —_———____——_- 
5 Lu (Iz—Ty) 


1 I7z'—Iz\}3 
-|-v’-10)(14 )| ; 
M Iz—TIy , 


ao) 
Se~fat a 


The moment of inertia about the symmetry axis 
of the original molecule Jz can not be evaluated 
from the rotational spectrum. However, Jz can 
be eliminated from these expressions using Eq. 
(15), giving 
1 (Ix’ — Ty’) Ix’ —Iz2’)}3 
ly| -|- need » (18a) 
we (Ix’ —Iy’ —Iz'+Ty) 
1 a a } 
's| -|- (Iy’ —Iy) (Iz’ —Iy) |. (18b) 
u (Ir’+Iz'—Ix'—Iy) 
Thus, the y and z coordinates of the substituted 
atom can be determined from Eq. (18) ;’ the 
distance of the substituted atom from the center 
of mass can be obtained from Eq. (14). It is also 
possible to calculate the moment of inertia of 
the original molecule about its symmetry axis 
from Eq. (15). 


VI. PLANAR ASYMMETRIC TOP MOLECULES 


Let the Z axis be chosen perpendicular to the 
plane of the molecule. Then Jx+Jy=Jz. With 
this choice of axes the Z’ axis will be parallel to 
the Z axis, and the z coordinates of all atoms 
zero. Equation (11) reduces to 


I;'= Iz+u(x?+y"). 


In its mathematical details the treatment here 


is similar to that given for symmetric tops in 
Sec. V(B). Equation (14) becomes 


y= uz! Iz), 
while Eq. (17) becomes 
x --- | 

M (Ix —TIy) 


1 Ix’ -—Ix\}3 
-|-s'=1)(14 )| » (20a) 
b Ix —TIy 
- ee) 
Les aa 
KB (Iy — Ix) 
1 Ty’ —_ Ty ? 
-|-ux’-10(14 )| _ (20b) 
M Iy—Ix 
7 If the isotopically substituted molecule is only a very 
slightly asymmetric top, often Jz’ can not be evaluated. 


In this case, since Eq. (18) is inapplicable, an explicit 
general solution can not be given. 


(19) 
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VII. NONPLANAR ASYMMETRIC TOPS 


In this case it is much more convenient to 
use the planar dyadic P defined in Eq. (5). In 
terms of this dyadic the moment equations 
[ Eqs. (10b) and (10c) ] for the original molecule 
become 


Pu®™ Lmxyi=0, 
Pa™ “—mnwz2;=0, 


Py.= LU imye2i=0, 
Px - Lmex?, 
Py= Lmiy?, 
Pz= Lime?, 


while for the substituted molecule the elements 
of P are of the form 


P,z/ =Pxt+ px’, 
P.,' = pxy, etc. 


(21a) 


(21b) 


(22) 


The secular equation for the substituted mole- 
cule, of which the roots are the planar principal 
moments Px’, Py’, Pz’, is 


Pxyt+ px? — 
pry 
The] 


wXY 
Py+ py? = 
LYS 


pxe 
wyS 
P2+ys?— 


This secular equation can be expanded to give a 
cubic equation in P’. It is well known from the 
theory of equations that there are definite rela- 
tions between the coefficients of a cubic equation 
and the roots of the equation. In particular, if 
Px’, Py’, Pz’ are the roots of a cubic equation 
in P’, the equation may be written 


— (Px'+Py'+P2')P?+ (Px'Py'+Px'Pz’ 
+Py'Pz')P’ —Px'Py'Pz'=0. 


By equating corresponding coefficients of the ex- 
panded secular equation to those in this expres- 
sion, three equations in x?, y?, 2? can be obtained. 
Solution of these three equations yields 


| -|- —e—e 
ls) Pe 


Y ~—) 
Pr—Pyr 


Pz'—Pz\} 
«(14+ —)]. (23) 
Pz—Px 


-|- (Px’ -Px)(14+5 


Analogous expressions for y and z can be ob- 
tained by cyclic permutation of the subscripts.*® 
The distance of the substituted atom from the 
center of mass of the original molecule can be 
determined from the expression 


|r| =[w-\(Px’ —Px+Py'—Py 
+Pz'—Pz) }}. 


Use of the relations given in Eq. (8) makes it 
possible to rewrite these expressions in terms of 
the axial moments of inertia 


(24) 


1 
in { —C(le!—Iv)-+ (a! —12) — (Ix! —Ix)] 


ae 
xX{| 1+-—-——_--_------—- --— 
2(Ix —Iy) 
(Ix’ — Ix) + Uy’ — Iv) — (I2’-—Iz) 1)! 
pS 


2(Ix—Iz) 
(25) 


|r| =((2u)* Ux’ — Ix+ Ir’ —Iy+Iz' —Iz)}}. (26) 


VIII. DETERMINATION OF MOLECULAR 
STRUCTURE 


The parameter equations derived above enable 
the position of an atom to be determined with 
respect to the center of mass of the molecule. 
If the Cartesian coordinates of two atoms have 
been determined, it is a relatively simple prob- 
lem in geometry to calculate the distance be- 
tween them. Hence, when sufficient data are 
available to determine the positions of all atoms 
in the molecule, the bond distances and bond 
angles of the molecule can easily be calculated. 

One isotopic substitution made on a particular 
atorn is all that is required to determine the posi- 
tion of that atom. An additional substitution 
made on the same atom should not give any new 
information on the position of the atom. How- 
ever, when effective moments of inertia are used 
instead of equilibrium values, the calculated 
structural parameters will be in error due to the 
neglect of zero-point vibrations. This will cause 
the parameters calculated from the additional 

8D. K. Coles, Advances in Electronics (Academic Press, 
New York, 1950), Vol. II, p. 317, has given approximate 
expressions accurate to order (Am/M)* for x, y, and z. 


Except for a numerical factor, his results can be shown to be 
in agreement with the exact expressions given here. 
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substitution to differ from those calculated from 
the first substitution. These differences are some- 
times useful in estimating the magnitude of the 
zero-point vibration effects. Some attempts have 
been made to correct for the zero-point vibra- 
tions by correlating these differences with the 
change in mass with isotopic substitution.® As 
yet no completely satisfactory method for mak- 
ing such corrections has been developed. 

Substitution on an atom near the center of 
mass of the molecule, which produces only a 
small change in the moment of inertia, may still 
cause a marked change in the vibrational fre- 
quencies and, hence, in the effects due to zero- 
point vibrations. It has even been found that 
in some cases the zero-point vibration effects are 
so large that, if neglected, imaginary values 
would be obtained for certain of the bond dis- 
tances.!° Consequently, it is always desirable to 
avoid substituting on an atom near the center 
of mass. 

It is apparent from the above that if the 
molecule is assumed to be rigid, i.e., zero-point 
vibrations can be neglected, a single substitution 
on a particular atom suffices to determine the 
position of that atom. The problem of deter- 
mining the minimum number of isotopic sub- 
stitutions (on different atoms) necessary to com- 
pletely determine the positions of all the atoms 
in the molecule is somewhat complicated. One 
complication arises from the fact that the pa- 
rameter equations only give information as to 
the magnitude of the coordinates of the substi- 
tuted atom, i.e., the signs of the coordinates are 
not determined. These ambiguities arise directly 
from the physical nature of the problem, viz., 
moments of inertia depend on length squared. 
In general, by making additional substitutions it 
would be possible to resolve these ambiguities. 
However, a much more profitable approach is to 
make use of structural information available from 
other sources, e.g., infrared spectroscopy, elec- 
tron diffraction measurements, theoretical chem- 
istry, etc., to resolve them. In the ensuing dis- 
cussion the attitude will be adopted that any 
ambiguities arising can be resolved by recourse 
to other sources of data. When this is not pos- 

® Townes, Holden, and Merritt, Phys. Rev. 74, 1113 
(1948). 


10D. K. Coles and R. H. Hughes, Phys. Rev. 76, 178 
(1949). 


sible, it will be necessary to make additional 
substitutions above the minimum number. 

Consider, then, an 2 atomic linear molecule 
in which n—2 isotopic substitutions have been 
made. If the molecular axis is taken to be along 
the Z axis, the magnitude of the z coordinates of 
the »—2 substituted atoms can be determined 
from the parameter equation [Eq. (13) ]. The signs 
of these coordinates are determined by recourse 
to other data. The magnitude and sign of the 
coordinates of two atoms are, then, still unde- 
termined. Substitution of the 2—2 known values 
into the moment equations for the original 
molecule, 


Y~me2;=0, and Ymez2=Ty, 


yields two equations in the two unknown co- 
ordinates. These two equations give two sets of 
solutions for the coordinates of the remaining 
two atoms: 2)’, 22’ and 2,’’, ze’, where, in general, 
2,’ 2," and 22’2z.2’’. By recourse to other data 
it should be possible to choose the correct one 
of these two sets. If only 7—3 substitutions were 
made, the moment equations of the original 
molecule would yield two equations in three 
unknowns. Solution would then be impossible. 
Hence, the minimum number of isotopic sub- 
stitutions required to determine completely the 
molecular configuration of a linear molecule is 
n—2. 

For a planar asymmetric top molecule, n—2 
isotopic substitutions would allow the magnitude 
of the x and y coordinates of these »—2 substi- 
tuted atoms to be determined from the parameter 
equation [Eq. (20)]. The signs of these co- 
ordinates would again be determined from other 
sources. The coordinates of the remaining two 
atoms can be determined by solving the five 
moment equations of the original molecule: 


Lmx:=0, Limx2=ly, YLmyi=0, 
Lmy?=Ix, Limxyi=I,=0. 


These five equations again yield two possible 
sets of solutions: x;’, x2’, yi’, ye’ and x1’, x2”, 
yi’, y2’’. By recourse to other data it should be 
possible to choose the correct one of the two sets. 
Consideration of the data obtained from »—3 
substitutions again indicates that the available 
information is not sufficient to fix the six co- 
ordinates of the remaining three unsubstituted 
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atoms. Hence, as for linear molecules, at least 
n—2 substitutions are required to completely 
determine the structural configuration. 

A similar analysis for a nonplanar asymmetric 
top molecule leads to quite different conclusions. 
In this case all nine moment equations in Eq. 
(10) are available—a somewhat simpler set can 
be obtained from the three in Eq. (10a) and the 
six in Eq. (21). If »—3 substitutions are made, 
the three Cartesian coordinates of each of the 
substituted atoms can be determined from pa- 
rameter Eqs: (23) or (25). The signs of these co- 
ordinates would again be determined from other 
sources. Substitution of these data into the nine 
moment equations yields nine equations in nine 
unknowns—the unknowns being the coordinates 
of the remaining three unsubstituted atoms. 
Since six of these equations are of the second 
degree and three of the first degree, an exact 
solution of these equations would be quite diffi- 
cult. It would probably be more expedient to use 
a numerical technique in solving these equations; 
the initial approximate values for such a process 
could be obtained from other sources of data. 
Hence, it should be possible to completely deter- 
mine the structural configuration from only n—3 
substitutions. Obviously if measurements from 
n—2 substitutions were available, the deter- 
mination would become vastly simpler and the 
need for accurate data from other sources greatly 
reduced. 

For all types of molecules, when information 
from n—1 substitutions is available, the deter- 
mination becomes still simpler. In this case 
data from other sources are helpful aids but are 
not absolutely necessary. Molecules possessing 
certain elements of symmetry will require an 
even smaller number of substitutions than the 
minima indicated above. Finally, it should be 
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mentioned that when effective moments of in- 
ertia are used and the calculations are first made 
using the minimum number of substitutions and 
then repeated using a larger number of substitu- 
tions, exact agreement should not be expected 
unless the effects of zero-point vibrations are 
taken into account. 


IX. MASS DIFFERENCE RATIOS 


If more than one isotopic substitution has been 
made on the same atom, the additional informa- 
tion can be used to determine mass difference 
ratios. Thus, if Jx, Jy, Iz represent the mo- 
ments of inertia of an original molecule of mass 
M, and Ix’, Iy’, Iz’ and Ix”, Iy’’, Iz’’ represent 
the moments of inertia of isotopically substituted 
species of this molecule in which the mass of the 
ith atom has been changed by Am’ and Am”, 
respectively, then from Eq. (26), 


r? = (2u’)'(Ix’ —Ixtly’ —Iy+Iz'’—TIz), 
a (2y"")—"( Ix” —Ix+ Ty" —Iy+ Is"—Is), 


where p’=MAm'/M+Am’ and yp” = MAm"/M 
+Am"’. Hence, the mass difference ratio is 
given by 


Am! = M'(Ix' —Ix+Ty'—Iy+Iz' —Iz) 


(27) 


where M’=M-+Am’ and M”=M+Am”. Analo- 
gous formulas can be obtained, for linear, sym- 
metric top, and planar ‘asymmetric top mole- 
cules from Eqs. (13) and (19). . 


ACKNOWLEDGMENT 


The author is very grateful to Dr. D. K. Coles 
for many helpful discussions. 


In the modern theory there is no longer any attempt to make an imaginative picture of the atom. 
An atom only gives evidence of its existence when it emits energy, and therefore experimental evidence 
can only be of changes of energy. The new theory takes over from Bohr the doctrine that the energy 
in an atom must have one of a discrete series of values involving h; each of these is called an ‘‘energy 
level.’,—BERTRAND RUSSELL, Human Knowledge, 1948. 
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The mean energy per excitation wavelength of e.m. radiation propagated along a lossless 
uniform wave guide is found to be w=Q», Q being independent of frequency for a given ampli- 
tude of the longitudinal field intensity, mode of the field, and structure of the guide. The fre- 
quency »v and the group velocity v, of the radiant energy are related by the equation v/v, 
=[1—(v,/c)?]-?, where »e=c/d-. denotes the frequency at cutoff, and , the excitation wave- 
length at cutoff. This equation permits one to put 


hv=hy.+hv-{ [1 — (v,/c)* }-#—-1}. 


The amplitudes of the field intensities are chosen such that Q equals Planck’s constant h. 
Applying the mass-energy relation, the term hvy,=hc/A:= moc? may be interpreted as the 
‘potential’ (=rest) energy of the quantum fv and moas the rest mass. The momentum  asso- 
ciated with hy is found to be p=h/)g, A, being the wavelength along the guide. When radiant 
energy passes from one guide into another one of different cross-sectional structure via a 
tapered piece of wave guide, a partial conversion of “potential” to “kinetic” energy, and vice 
versa, occurs along the tapered section. In the case of TEM waves, which, however, can only 


be supported in a guide of infinite transverse dimensions, the group velocity becomes c and the 
rest mass zero. 





N an interesting paper! on particle and quan- 

tum properties of electromagnetic wave fields 
in free space, H. W. KGnig attempts a classical 
interpretation of the quantum nature of light 
which allows the possibility of a rest mass for 
the photon. His conclusions, however, are based 
on several arbitrary assumptions and do not 
appear to be consistent with results of non- 
classical theory (see reference 1, page 420). 

It is shown in the following that such particle 
properties of electromagnetic wave fields can be 
derived classically without any arbitrary as- 
sumptions and in accordance with established 
theories, if a geometrical structure, such as a 
wave guide, imposing periodicity requirements 
on the field, is present. 

For the sake of simplicity, a uniform non- 
dissipative wave guide of rectangular cross sec- 


equations (MKS units) ? 
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' 


marx nry 
H,= —Ak*® cos— - cos(wt — Bz). 
a 


In Eqs. (1) A, in amp-m, is an arbitrary con- 


tion may be considered. The perfectly conducting 
walls are the planes x=0 and x=a, y=0 and 
y=b (Fig. 1). The z axis is the axis of propaga- 
tion along which H,, , waves may be propagated. 
The field distribution is given by the following 


* Now at Sir George Williams College, Montreal, Canada. 
1H. W. Kénig, Acta Physica Austriaca 4, 405 (1951). 
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stant determining the amplitude, and m and n 
are integers defining the mode of transmission 
by indicating the number of half-period varia- 
tions in transverse field intensity along the x 
and y dimensions of the guide, respectively. The 

2L. G. Huxley, The Principles and Practice of Wave 


Guides (Cambridge University Press, Cambridge, 1947), 
p. 52. 
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Fic. 1. Illustrating the notation used with 
rectangular guides. 


other symbols used in these equations and in 
the further discussion are as follows: E, H are 
electric and magnetic field intensity, respec- 
tively, v=w/2m frequency, v, the so-called ‘‘cut- 
off frequency” (see below), A. wavelength at 
cut-off frequency, \ wavelength of field excita- 
tion or free space wavelength, \, wavelength 
along the guide, v, phase velocity, v, group 
velocity, €, u dielectric constant and permeability 
in free space, respectively, and c=(eu)~? ve- 
locity of light in free space. The quantity k is 
given by k?=(mz/a)?+(n2/b)?=(27/X.)?, and 
B=2x/d,. For any mode of transmission, the 
cut-off frequency is the lowest possible fre- 
quency at which such mode is possible; the 
cut-off wavelength \, is then given in terms of 
the guide dimensions a and b by 


(1/d.)? = (m/2a)?+ (n/2b)?. (2) 


The quantities A, \,, and A, are related by the 
equation 


(1/d.)? = (1/A)?— (1/A,)?, 
or, multiplied by (1/7), 
(1/vd.)? = (1/c)?—(1/vp)?, (4) 


where v»=A,v and c=)v. Frequency and wave- 
length at cutoff are linked by the equation 


(3) 


C=AXPe. 


(S) 

It is important to note that Eqs. (3), (4), and 
(5) hold quite generally for any kind of uniform 
guide and mode.* The choice of the cross-sec- 
tional structure (boundary conditions) and mode 


3N. Marcuwitz, Wave Guide Handbook, M.1.T. (Mc- 
Graw-Hill Book Company, Inc., New York, 1951), p. 8. 


of field only affect the value of \,. As c and XA, 
are independent of », it follows that the phase 
velocity v, depends on frequency, and the propa- 
gation is accompanied by dispersion. The 
group velocity is defined by 1/v,=(d/dv)(v/v,), 
which yields, with Eq. (4), 


UpV,=C*. 


(6) 


The mean energy W per unit length in the guide 
is found in the usual way to be 


W = A*kabwupr,/8e?, (7) 


and the total average power flow P along the 
guide in the z direction is 


P=A*kabwyB/8, (8) 


or, with Eqs. (6) and (7), P= Wz,, which states 
that radiant energy in the guide is transmitted 
with group velocity. The mean energy w per ex- 
citation wavelength X is given by 


w=)A?RabwyBv,/8c’, 
or rearranged, 


w= (Ak*abur?/2c) v= Qy, 
where 
Q=A*kabyur?/2c. (9) 


It is obvious that Q is a constant independent of 
v for a given value of A, mode of field and ge- 
ometry of guide, and that Q has the dimensions 
of Planck’s constant h. The radiant energy w, 
therefore, may be expressed in terms of hv, giving 


(10) 


n being defined by »=Q/h. As w represents the 
mean energy per excitation wavelength along 
the guide, the numeral m is not necessarily an 
integer. 

Since a signal, i.e., finite amount of energy, 
cannot be transmitted with a strictly mono- 
chromatic wave, the equation w=Qv may be 
regarded as the definition of a mean fre- 
quency of a group of waves, if in Q an appropri- 
ate average value of A is substituted. (Also, vice 
versa, using an appropriate average of v, w=Qyv 
may be regarded as the definition of Q and, 
thence, by Eq. (9) of A). Instead of the total 
cross-sectional area of the guide, one could 
choose equally well the smaller area between two 
pairs of planes within the guide, which are 


w=nhvp, 
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parallel to the actual walls of the guide and 
along which the necessary boundary conditions 
are fulfilled, e.g., x=a/m and x=2a/m, y=b/n 
and y=2b/n. Apart from a reduction of the value 
of Q, the decisive equations would not be 
affected. 

In order to avoid the factor m in the following 
equations, let the special case now be considered 
where the arbitrary constant A in the expression 
for Q has such value that Q becomes equal to 
Planck’s constant 4; hence m=1 and w=hp. 
The value of A is determined by the power 
output of the generator at the sending end of 
the guide. Since the output of the generator can 
be attentuated at will, the integration constant 
A may be regarded as arbitrary. It is obvious 
that a variation of frequency will not affect the 
value of A, if during such variation not the 
power output and not the amplitude of the 
electric or transverse magnetic field is kept 
constant, but rather the amplitude of the longi- 
tudinal magnetic field H,. 

According to Maxwell’s theory, the z com- 
ponent of the electromagnetic momentum ?, is 
originally defined by the equation 


a@ pb avr 
ep. f f [[ Savayas, 
0 0 0 


S, being the mean z component of Poynting’s 
vector. Employing Eqs. (8), (9), and (10) with 
n=1, it is readily found that 


On the other hand, the z component of the 
momentum of the mass m which is associated 
with the quantum of energy w=hy is equal to 


(12) 


since w, as an element of W, moves with group 
velocity v,. Equalizing Eqs. (11) and (12) and 
substituting c?/v, for v,, one obtains 


m=hv/c?. (13) 


Equations (11) and (13) serve to determine mass 
and momentum associated with the electro- 
magnetic energy hy streaming along the guide. 
The next step consists in referring the wave 
field, represented by Eqs. (1), to a coordinate 
system Ko(X, Y, Zo) that moves with group ve- 


pP.= M2, 


locity v, in the direction of propagation relative 
to the coordinate system K(X, Y, Z) of the guide. 
With respect to Ky the group velocity of the 
wave field becomes zero and the phase velocity 
(vp)o, because of Eq. (6), infinitely large. The 
frequency vo observed in Ko is defined by Eq. 
(4). Putting (v,)o= © and c=X,.»,-, one obtains 
Yo= ve, in accordance with the definition of the 
cut-off frequency. For the wavelength along the 
guide one gets (A,)o= ©, because of (vp)o= vo(Ag)o 
= ©. The mode of the wave field observed in Ko 
is obtained by putting B=0 and w=ao=«, in 
the field equations [Eqs. (1) ]. Thus, 


ne mmx MTy | 
(Ez)o= sno cos—— sin— sinwfo, 
a 


s mr mmx nry 
(E,)o=Aw-u— sin—— cos— sinw¢lo, 


a 
(E.)o=0, 
(Hz)o= (Hy)o=0, 


max mnry 
(H.)9= —Ak*® cos—— cos—— cosw fo. 
a 


Obviously there is no phase dependence on the 
axial distance; hence, there is no progressive 
wave in Ko. A given field component oscillates 
everywhere in phase with itself. The components 
of E vibrate in quadrature with the z component 
of H. There is no power flow along the 2» axis. 
This, incidentally, is the mode of field which 
De Broglie associated with a particle at rest.‘ 
It is worth noting that the Lorentz transforma- 
tion applied to the field equations [Eqs. (1) ] for 
the transition from K to Ko leads exactly to 
Eqs. (14). The excitation wavelength observed 
in Ky becomes \,, and the energy per excitation 
wavelength becomes wo=hy-. 

Mass and momentum associated with wo ap- 
pear to be, respectively, 


mo=(hve)/c?=h/(cd-), 


po=0. (15) 


Using Eqs. (5) and (6) one readily finds from 
Eq. (4) that 
v/ve=[1—(v,/c)? 1+. (16) 


4A. March, Quantum Mechanics of Particles and Wave 
Fields (Chapman and Hall, Ltd., London, 1951), p. 12. 
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Dividing Eq. (13) by Eq. (15) gives the rela- 
tionship between mass in motion and mass at 
rest, m=m)[1—(v,/c)?}-', a result which may 
be obtained also by subjecting energy and mo- 
mentum to the Lorentz transformation. 

Product hv represents the total quantity of 
energy w in K, relative to which it moves with 
group velocity, whereas hv, represents the energy 
of this energy “quantum” observed in a refer- 
ence system relative to which it is at rest. There- 
fore, hv. may be interpreted as the potential5 
energy w, of this ‘‘quantum.”’ The kinetic energy 
WwW, appears to be 


W,=W—Wyp=hv—hv.=hv[ (v/v) —1], 
or with Eqs. (16) and (15), 
wi =hve{[1 —(v,/c)?}-*—1} 
=moc*{(1 — (v,/c)? }-#— 1}. 
Equation (15) states that the rest mass, and 
thus the potential energy moc*, varies with the 
cut-off wavelength, which depends on the struc- 


ture of the guide and the mode of the field. 
This is understandable, in view of the fact that 


the group velocity of the electromagnetic energy 
in a guide depends on X,, according to Eq. (4), 
which by Eq. (6) may be rewritten as 


vg=cL1—(A/d-)?}}. (17) 


When radiant energy passes from one guide into 
another one of larger cross-sectional area via a 
tapered piece of guide, so that the mode of trans- 
mission does not change, the group velocity of 
the radiation and hence its kinetic energy in 


5 The term “potential energy’’ w, is used here for what 
in the theory of relativity is conventionally called ‘rest 
energy” Wo; Wp is the maximum energy the quantum can 
possess when at rest. 


the larger guide, is greater. An increase in 
kinetic energy, however, can only occur at the 
expense of potential energy, as the total energy 
must be preserved. The tapered guide section 
apparently acts as a transformer, converting 
potential energy of the radiation partially into 
kinetic energy and vice versa, according to the 
direction of propagation. An expansion of the 
electromagnetic field, however, reduces the field 
intensities. Therefore, the quantity Q=d,P/v 
=)?P/v, in the larger guide will also be smaller, 
and hence <1. This means that the section of 
the larger guide which contains the mean energy 
w=hy will be longer by a factor 1/n. In other 
words, the mean axial distance between two suc- 
cessive quanta is larger in the wider guide as a 
consequence of the acceleration occurring along 
the tapered section. 

The dominant mode in free space, which may 
be regarded as a uniform wave guide of infinite 


transverse dimensions? is the ‘‘TEM mode” . 


(transverse electromagnetic waves), which a hol- 
low pipe wave guide is not capable of supporting. 
In this case the cut-off wavelength becomes 
infinitely large, and therefore, 


Vp=V,=C, v-=0 and w,=0, mo=0. 


All the potential energy changes into kinetic 
energy as the radiation passes from confined 
space into free space. A quantum traveling at 
speed c has attained its maximum speed and 
possesses neither potential energy nor rest mass. 

In conclusion it may be pointed out that the 
group velocity is almost equal to c, if A>d. 
Unless the cross-sectional dimensions of the 
guide are microscopical, it is immaterial whether 
the Compton effect, for example, is studied with 
or without the employment of wave guides. 


Functional Writing 


At Colgate University, English department consultants take an opaque projector into 
science laboratories and other classrooms and show on a screen samples of student reports 
written for these courses. The papers are discussed by the class and criticized both by the 
English and the subject-matter teachers. This is part of the Colgate plan of functional writing, 
in which emphasis is placed on accurate and concise expression in all courses. “If it isn’t good 
English it isn’t good physics either’’ is the slogan. 
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In the light of private answers suggested for the public questions: What is physics? What is 
education?, a critical review is made of fluid dynamics in physics teaching. First, general 
physics is examined for what is being done in teaching fluid dynamics in physics. A critical 
analysis is made in particular of the important Bernoulli theorem as to what it is, what it says, 
what it means, and how it is used. Second, theoretical physics is examined for what is being 
done. In this case, the motion of a sphere has been selected as illustrative of present teaching. 
The typical course in theoretical physics is then examined with respect to its material about 
the motion of a sphere, both in a perfectly inviscid fluid and in a perfectly viscous fluid. This 
academic material is compared finally with experimental observations, particularly those in 
wind tunnels and in ballistics ranges. In conclusion, questions are noted on the foundations of 
fluid dynamics, particularly the role of viscosity, turbulence, the concept of a steady state, and 
the role of heat. The review ends with a moral question as to the students’ stake. 


PHYSICS AND EDUCATION 


HE subject of this discussion revolves about 

two fundamental questions. What is phys- 

ics? What is education? In general, these ques- 
tions lead to a third: how to teach physics, or 
more specifically, how to teach fluid dynamics in 
physics? In this connection, there is an interest- 
ing comment by F. K. Richtmyer' in the first 
article in the American Physics Teacher, entitled 
“Physics is Physics,” with reference to a remark 
by G. W. Jones, Professor of Mathematics at 
Cornell University: ‘‘to become a_ successful 


teacher of mathematics one must acquire a 
thorough knowledge of mathematics—and then 
by constant practice acquire the art of teaching 
t.”’ (The latter part of this statement is due to 


Richtmyer.) He states further ‘‘teaching is an 
art and not a science.’’ On the other hand, H. G. 
Wells has emphasized, ‘‘No man can be a good 
teacher when his subject becomes inexplicable.” 
Here we are faced with the dual implication of 
our question, namely, that science signifies es- 
sentially knowledge derived from nature, and 
art in this case is the method of its expression. 
To the above public questions every teacher 
has rightfully his own private answers, which 
* Address given at a symposium sponsored jointly by 
the American Association of Physics Teachers, the Fluid 
Dynamics Division of the American Physical Society, 


and the Physics Division of the American Society for 
Engineering Education, at Michigan State College, June, 
1951. 


7 At present with the National Science Foundation, 
Washington aap ©. 


IF. K K. Richtmyer, Am. Phys. Teacher 1, 1 (1933). - 


form the axiomatic foundation for all his in- 
struction. Accordingly, I must state at once my 
own views as to the primary goal of physics 
teaching and the methods of its attainment. 

Someone once asked the artist H. Matisse how 
he saw an orange. ‘‘That depends,” he answered, 
“Yon whether I am buying it, eating it, or paint- 
ing it.”” Likewise, our outlook on physics depends 
upon our point of view. Evidently any matter, 
including physics, can be regarded from the view- 
points of aesthetic enjoyment, of rational in- 
quiry, and of technological use. Any integrated 
education must involve all these aspects even 
though it may stress only one. Thus, A. N. 
Whitehead in his The Aims of Education says, 
“There is not one course of study which merely 
gives general knowledge and another which gives 
specific knowledge.’’? Even so-called pure science 
and applied science are intimately interrelated, 
as F. Klein emphasized long ago. W. F. G. 
Swann in his recent article ‘“Teaching of Phys- 
ics’ observes further that, ‘‘Physics has com- 
bined two phases, knowledge of facts and knowl- 
edge of ideas.’”" He proceeds to emphasize the 
ideas aspect. G. F. B. Riemann, too, has stressed 
the latter in the statement, ‘‘Science is the at- 
tempt to comprehend nature by means of con- 
cepts.”” In the same spirit physics, to me, is pri- 
marily the understanding of physical phenomena 
through concepts. 


2A. N. Whitehead, The Aims of Education (Macmillan 
Company, New York, 1929). 
3W. F. G. Swann, Am. J. Phys. 19, 182 (1951). 
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TABLE I. Distribution”offfluid dynamics topics in 46 
general physics textbooks. 


Topics Number of texts 
Bernoulli theorem 
Applications 
Torricelli theorem 
Vena contracta 

Viscosity 
Friction head 
Streamlining 
Hydraulic ram 
Poiseuille law 
Turbulence 
Water wheel 
Venturi meter 
Stokes law 
Vortex motion 
Reynolds number 


Whitehead suggests further a method of at- 
taining such a goal. He mentions ‘“‘two educa- 
tional commandments: do not teach too many 
subjects, and again, what you teach, teach 
thoroughly.’ Thus, our first imperative is to 
select the elemental concepts and the associated 
phenomena. Physics, therefore, must be re- 
garded not as a dogmatic presentation of a self- 
sufficient body of learning, but rather as a 


creative search for a predictive description of 
nature. It should strive continually for libera- 
tion from superstitions, both old and new. Thus, 
physics teaching of fluid dynamics should stimu- 
late students to think critically about its prob- 


lems. Breadth, moreover, is to be achieved 
through coherent ties rather than through wide- 
spread spottiness, as e pluribus unum should 
signify an orchestration rather than a melting 
pot. Above all, may I suggest that quality in 
teaching, as in acoustics, depends upon over- 
tones and not just upon fundamentals. 

Education, in short, is not so much like a static 
structure, which is fixed according to a blue 
print; it is more like an organic plant, which 
grows up and down and spreads out. From this 
point of view let us now regard fluid dynamics 
teaching at two levels: that of undergraduate 
general physics and that of graduate theoretical 
physics. 


FLUID DYNAMICS IN GENERAL PHYSICS 


First of all, let us review what is done in 
general physics courses by glancing at some ele- 
mentary textbooks that have appeared between 


SEEGER 


1930 and 1951. It is significant that in fifteen 
textbooks (522 pages average) there is abso- 
lutely no mention at all of any topics in fluid 
dynamics, although the textbooks themselves 
have inviting titles like Everyday Physics; Man 
and his World; The Universe; Atoms, Rocks, and 
Galaxies; From Galileo to Cosmic Rays, etc. In 
my personal opinion the quality itself in some of 
these cases is poor because of the very overtones 
that are present or absent. For example, in one 
textbook there is a mention of a so-called 
Aristotelian law of falling bodies, but no word at 
all about Archimedes principle. Such a selection, 
I believe, gives a perverted idea as to the status 
of Greek science. In another textbook there is a 
description of the measurement of atmospheric 
pressure but no indication at all as to any physi- 
cal meaning of the term pressure—not even a 
guiding definition. 

In those 46 textbooks (634 pages average) 
which devote some part to fluid dynamics (in- 
cluding viscosity) there is an average of nine 
pages (1.42 percent) used for this subject. The 
impression gained from this cursory review of 
college physics textbooks is reflected in the 
American Journal of Physics, the organ of physics 
teachers. In its entire history since 1933 there 
have been only 27 articles and 21 letters (152 
pages out of 6955 pages total, approximately 
2.19 percent) that can be construed in any sense 
to be concerned with fluid dynamics. (As a 
matter of fact, three of the articles (57 pages) 
are actually on rockets.) Evidently, this Journal, 
as well as the textbooks, reflects merely the in- 
terest of physics teachers themselves, as empha- 
sized also in the popular histories of physics. 
For example, F. Cajori‘* fails to mention the 
fluid dynamics contributions of L. Euler and 
D. Bernoulli in the 18th century, and of O. 
Reynolds, M. H. Navier, and G. G. Stokes in 
the 19th century. As a matter of fact, there is no 
reference to turbulence including Reynolds’ in- 
vestigations either in this book or in his com- 
plementary history of mathematics.® All in all, 
fluid dynamics has not been particularly inter- 
esting to physicists themselves. 

4F. Cajori, A History of Physics (Macmillan Company, 
New York, 1929), rev. edition. 


5 F. Cajori, A History of Mathematics (Macmillan Com- 
pany, New York, 1919), second edition. 
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Let us now turn to the topics of fluid dynamics 
that are actually selected in general physics 
textbooks. The findings are listed in Table I. 

It is not surprising to find that flows in pipes 
and in channels, which are everyday experi- 
ences, are among the popular topics. The initial 
scientific achievement along this line was the 
series of investigations of J. L. M. Poiseuille 
(see simultaneous work of G. Hagen) begun in 
1835, announced in 1840, and published in 1846.°® 
Specifically, Poiseuille observed that the quan- 
tity of fluid flowing from a tube is directly pro- 
portional to the excess pressure across its ends 
and the fourth power of its diameter, and in- 
versely proportional to its length. It is interest- 
ing that the motivation for this work was con- 
cerned originally with the problem of blood 
circulation in capillaries, but that later studies’ 
show Poiseuille’s results not to be valid for fine 
capillaries (about a 3-mm diameter is the lower 
limit for blood flow). 

It is difficult to understand completely the 
pedagogical reason for the inclusion of the 
Poiseuille law in an elementary textbook. Per- 
haps the simplified theoretical deduction that 
can be made of the law is the primary objective. 
Surprisingly enough, few authors note that 
whereas Stokes undoubtedly derived the law 
(to date, one of the few existing solutions of the 
Navier-Stokes equations with boundary condi- 
tions), he himself never published his proof, 
apparently because of its disagreement with the 
available experimental data (primarily for large 
pipes). Accordingly, it is the later published 
derivation of G. Wiedemann and of E. Hagen- 
bach® that is the one commonly used in text- 
books today. What is truly disappointing, how- 
ever, in an elementary textbook is the failure to 
stress—or even to note—that Poiseuille’s original 
work is a classical example of the observational 
phase of the inductive method. In unparalleled 
fashion he isolated the important physical vari- 
ables and analyzed the sources of error. His 
resulting empirical law for water flowing through 


6 J. L. M. Poiseuille, Mémoires des Savants Etrangers 
ti 339) (1846). G. Hagen, Pogg. Ann. Phys. Chem. 46, 423 
1 ; 
7G. Barr, Viscomeiry (Oxford University Press, New 
York and London, 1931). 
8G. Wiedemann, Pogg. Ann. Phys. Chem. 99, 217 
teeny E. Hagenbach, Pogg. Ann. Phys. Chem. 109, 385 
1 k 


fine capillaries (0.01 to 0.65 mm) had an accuracy 
between 0.1 and 0.2 percent. Students should be 
encouraged to read for themselves a translation 
of the original article of Poiseuille’—not to men- 
tion the equally stimulating presentations of B. 
Pascal.!° 

As early as 1854 and 1869 Poiseuille pointed 
out also that there are definite limitations to the 
validity of his formula. The classical work in this 
direction, of course, is that of O. Reynolds,"! who 
began experimenting in 1880 with pipes } in., 
3 in., and 1 in. in diameter and 4 ft 6 in. long. 
The work of Reynolds is concerned primarily 
with the limitation due to the initiation of tur- 
bulence (see Fig. 1). It is now known that 
Poiseuille’s law for laminar flow is strictly valid 
only up to a certain “‘critical’’ Reynolds’ number 
R, of about 1000, although for smooth rounded 
pipes with special entrances the law may be 
valid up to a Reynolds number of 75 000. Two 
other physical factors may be important also. 
If, for instance, there is slipping along the walls, 
the quantity of fluid is classically proportional to 
the factor” (1+4y/AR), where yp is the coefficient 
of viscosity, A a slip coefficient, and R the 


Reynolds number. Furthermore, as the pressure 
of a gas decreases so that molecular structure 
becomes significant, one has to multiply by the 
so-called Knudsen™ factor (ratio of molecular 
mean free path to a body length, such as the 
diameter of the tube). All the above remarks 


Fic. 1. Reynolds’ data for the nonlinear, variation of 
the quantity of liquid issuing from a tube with the speed 
(see reference 11). 


® Rheological Memoirs: J. L. M. Poiseuille, Experi- 
mental Investigations upon the Flow of Liquids in Tubes of 
Very Small Diameters, translation by W. H. Herschel 
(Easton, Pennsylvania 1, V940). 

10 The Physical Treatises of Pascal, translation by I. H. B. 
and A. G. H. Spiers (Columbia University Press, New 
York, 1937). 

110, Reynolds, Trans. Roy. Soc. (London) 174, 935 
(1883); Papers on Mechanics and Physics Subjects I, 51 
(Cambridge, 1900-1903). 

12H. A. Lorentz, Lectures on Theoretical ~— (Mac- 
millan Company, New York, 1927), Vol. 1, p. 81. 

18 M. Knudsen, Ann. Physik 28, 75 (1909) ; ; The Kinetic 
Theory of Gases (Methuen and Company, London, 1946), 
second edition. 
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point up to a severe criticism of so many modern 
textbooks in that experimental information is 
rarely given completely or even referenced by a 
caution sign. In the case of the Poiseuille law, 
moreover, the historical significance of some of 
the great investigations in physics has also been 
neglected. 

Even more popular than the Poiseuille law, 
it seems, is the vena contracta. Perhaps, the 
appeal here is the gobbledygook one of “‘far 
away places with strange sounding names.’’ One 
is reminded of H. C. Andersen’s story of the 
‘Nightingale’: ‘“The music master wrote a 
twenty-five volume about the artificial bird. It 
was learned, long-winded, and full of hard 
Chinese words. Yet everyone said they had read 
and understood it, lest they should show them- 
selves stupid.’’ Sometimes the phrase vena con- 
tracta is used just as a label, occasionally it is 
defined, but hardly ever is one told how the 
concept originated or why it is now important. 

The practical Romans, for example, were par- 
ticularly concerned with the vena contracta in 
connection with their problem of increasing the 
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Fic. 2. Title page of Bernoulli’s Hydrodynamica 
(see reference 17). 
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rate of flow from orifices connected to their 
water supply. The inclusion of such a concept in 
an elementary course illustrates Whitehead’s 
“inert ideas’: ‘‘Above all things, we must be- 
ware of what I will call ‘inert ideas’—that is to 
say, ideas that are merely received into the 
mind without being utilized, or tested, or thrown 
into fresh combinations.’”? 

Finally, the selection of phenomena and of 
ideas for the fluid dynamics part of a general 
physics course is not obvious and certainly not 
modern. At best, it may be said to be classical. 
First class in education, however, may not ne- 
cessarily mean the highest. 


BERNOULLI’S THEOREM IN GENERAL PHYSICS 


In view of the widespread consideration given 
to Bernoulli’s so-called theorem, as well as its 
basic engineering importance, we shall consider 
this topic in some detail. A. Sommerfeld" insists : 
“Bernoulli's equation is the most important 
theorem in elementary fluid dynamics,” but 
F. Magie does not bother even to mention it in 
his source book in physics.!® An educational argu- 
ment for its inclusion in an elementary course is 
that it is one of the few simple concepts con- 
cerned with the mechanics of a continuum as 
distinct from that of physical bodies. In this 
respect, it is strange that J. C. Maxwell in his 
Matter and Motion (1877) does not discuss fluids 
at all. 

The first question that arises is whether the 
Bernoulli statement itself is essentially a theo- 
rem, a principle, or a law. First of all, what isa 
theorem? Is it ‘‘a general statement that has 
been proved’’'* particularly with respect to 
mathematics? On this basis it would appear 
historically, at least, that the relation is a theo- 
rem; for D. Bernoulli in his classical work!’ on 
hydrodynamics (see Fig. 2) utilized the principle 
of vis viva, first employed by Huygens, in order 
to show that the change in kinetic energy for a 
fluid element displaced during a constant dis- 
charge is equal to the work done on the element. 

144A. Sommerfeld, Mechanics of Deformable Bodies (Lec- 


a Theoretical Physics) (Acadernic Press, Inc., 1950), 
ol. 


16 W.F. Magie, A Source Book in Physics (McGraw-Hill 
Book Company, Inc., New York, 1935). 

16 Webster's New International Dictionary (G. & C. 
Merriam Company, Springfield, 1949), second edition. 

17D—D. Bernoulli, Hydrodynamica (Argentorati, 1738). 
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Fic. 3a. Text from Bernoulli’s Hydrodynamica (see reference 17). 


Inasmuch as Bernoulli used this vis viva principle 
throughout his entire book, it is difficult to locate 
the specific place where the so-called Bernoulli 
theorem can be said to be first introduced. My 
colleague, G. Kuerti, believes that page 2587 of 
the Hydrodynamica is a likely location (see 
Fig. 3). 

Secondly, what is a principle? Is it ‘‘a compre- 
hensive law from which others are derived, or 
upon which others are founded’’?!* On this basis 
the mere fact that many textbooks derive the 
equally famous Torricelli theorem from the 
Bernoulli theorem would justify calling the latter 
a principle. Historically, however, this interpre- 
tation is not correct. E. Torricelli!* himself ob- 
tained his law directly from observations of 
liquids flowing from orifices (here the concept of 
vena contracta is absolutely essential for the 
correct statement of the theorem). Indeed, the 
historical development of fluid dynamics con- 
cepts is from Torricelli through Bernoulli to 


18 E. Torricelli, Lezioni Accademiche (Opere di Evangelista 
Torricelli) (Faenza, 1919), Vol. 2, p. 185. 


Euler, whereas we usually trace the logical de- 
velopment from Euler through Bernoulli to 
Torricelli. It is for this reason that authors are 
at variance even as to the so-called ‘Father of 
Hydrodynamics.”” Thus, E. Mach’® so desig- 


Fic. 3b. Diagram from Bernoulli’s Hydrodynamica (see 
reference 17). 


19. Mach, The Science of Mechanics (Open Court, La 
Salle, 1942), fifth edition. 





as. 


Fic. 4, Farwell’s “Two dia rams appearing in different 
current textbooks”’ (see reference 22). 


nates Torricelli, whereas E. T. Bell?® points out 
that Bernoulli is often called the ‘Father of 
Mathematical Physics’; L. Prandtl and O. G. 
Tietjens#! choose Euler. 

Finally, what is a scientific law? Is it ‘‘the 
observed regularity in nature’’?!® If so, the 
Bernoulli relation is wanting in precise con- 
firmation owing to the existence of frictional 
effects, which are usually ignored in the theo- 
retical treatment of an idealized fluid. Inci- 
dentally, even when external friction is actually 
emphasized in an elementary course, another 
difficulty often arises from failure to specify the 
magnitude of the velocity head (see Fig. 4).” 

The confusion indicated above illustrates, I 
believe, P. Frank’s criticism of present-day teach- 
ing in his discussion of the ‘‘Philosophy of Sci- 
ence in the Physics Curriculum” 2% “The result 
of conventional science teaching has not been a 
critically minded type of scientist, but just the 
opposite---main source is the failure to teach 
students a clear distinction between the observa- 
tion of facts and the drawing of logical conclu- 
sions.” Above all, it is important for the teacher 
to discriminate between the roles of induction 
and of deduction in scientific research.** A valu- 
able educational overtone is their interplay in 
the development of a useful theory. 

There is an apparent uncertainty in an intro- 
ductory course not only as to how to classify 
the Bernoulli relation, but even as to how to 
state it. Typical differences of opinion are shown 
by the following excerpts from three well-known 
texts: 


20E. T. Bell, Men of Mathematics (Simon & Schuster, 
New York, 1937). 

27. Prandtl and O. G. Tietjens, Fundamentals of 
Hydro- and Aerodynamics (McGraw-Hill Book Company, 
Inc., New York, 1934). 

2H. W. Farwell, Am. J. Phys. 12, 307 (1944). 

23 P, Frank, Am. J. Phys. 15, 202 (1947). 

2 R. von Mises, Positivism (Harvard University Press, 
Cambridge, 1952). 
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Bernoulli’s theorem, which applies only to an incom- 
pressible and nonviscous fluid. j 


Because liquids have so little compressibility, Ber- 


noulli’s principle holds more precisely for them than 
for gases. 


The incompressible fluid Bernoulli equation can be 
interpreted as an energy equation, but the compres- 
sible fluid equation is only a momentum equation. 


The point at issue here is not merely the se- 
mantics of a label, but actually the physics of its 
content. There is certainly need for reviewing 
elementary concepts of physics, too, from an ad-: 
vanced point of view.?5 In this connection, let 
us consider various pressure-speed relations per- 
tinent to common fluid flows. For simplicity we 
shall confine our attention to the one-dimensional 
case. 

We start with the basic dynamic relations for 
changes in mass (continuity equation), in mo- 
mentum (force equation), and in energy (con- 
servation of energy). One must distinguish be- 
tween the equations of motion for two different 
ideal fluids, namely, for the perfectly inviscid 
incompressible fluid of L. Euler®* and for the per- 
fectly viscous incompressible fluid of C. L. M. H. 
Navier and of G. G. Stokes.?’ In both cases the 
equations can be readily modified for compres- 
sible fluids. In general, inclusion of viscous effects 
(for example, the coefficient of viscosity) necessi- 
tates consideration also of thermal conduction, 
inasmuch as the so-called Prandtl number (re- 
lating viscous and thermal conductive effects) 
is usually of the order of unity. Only if thermal 
conduction is negligible, moreover, does the en- 
ergy equation reduce to the mechanical energy 
equation, which originates from a space integral 
of the force equation, provided all forces are 
included. For example, the time-rate of change 
of the mechanical energy (T kinetic, V potential) 
of a given quantity of a perfectly inviscid com- 

23F, Klein, Elementary Mathematics from an Advanced 
Point of View, translation by E. R. Hedrick and C. A. 


Noble. Arithmetic, Algebra, Analysis (Macmillan Com- 
any, New York, 1932); Geometry (Macmillan Company, 
ew York, 1939). 
2%6L. Euler, aol Berlin (1755). J. L. de LaGrange, 

1788). Acad., Berlin (1781); Mécanique Analytique (Paris, 
27C. L..M. H. Navier, Memoire de l’Académie Royal 

des Sciences 6, 389 (1823). G. G. Stokes, Trans. Cambridge 

Phil. Soc. 8, 287 (1845). 
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pressible fluid is given by 
d/dt(T+ V)= foa-do+ f pv-adr, 


where > is the pressure and q the fluid velocity 
(vector) in a region of volume v and of external 
surfaces s. The energy equation can be modified 
to include thermodynamic changes. 

Pressure-speed relations along a streamline for 
one-dimensional steady (no time variation at a 
fixed point in space) motion are given below for 
different kinds of fluids at a height z above the 
level of zero potential energy: 

(a) Inviscid, incompressible fluid, 


p 
3¢’+g2+-=A. 
p 
(b) Inviscid, compressible fluid, 
3g°+es+ f dp/p=B, 


where the pressure is a function of the density 
only (so-called polytropic changes). If p=kp?, 
then 


hott g2+—— a8. 
¥~t% 
For a perfect gas with negligible viscosity, y is 
the adiabatic exponent (ratio of the specific heat 
at constant pressure to the specific heat at con- 
stant volume). 
(c) Viscous, incompressible fluid, 


_ 3P+g2+p/p=A+rF, 


where » is the kinematic coefficient of viscosity 
(ratio of coefficient of viscosity to density) and 
F is a function of both the fluid rotation and the 
streamline curvature. 

The parameter A (and B) is constant only 
along a particular streamline. Its gradient, how- 
ever, can be shown to be equal to the sum of a 
term involving rotation and of another term in- 
-volving the entropy. Thus, for irrotational, isen- 
tropic flow the value of such a parameter is uni- 
form for all streamlines (this fact is sometimes 
called the strong form of Bernoulli’s theorem). 

For unsteady irrotational motion of an in- 
-viscid compressible fluid the following relation 
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holds: 
0¢ 
at’ 


dp 
ag +gs+ f —=C+ 
p 


where the value of C is uniform throughout the 
fluid and where ¢ is the velocity potential. 

The uncertainty about the Bernoulli theorem 
in a general physics course extends not only to its 
mathematical formulation, but also to its physi- 
cal interpretation. It revolves largely about the 
phrase ‘“‘pressure energy,’’ which is sometimes 
used to designate the term involving the pres- 
sure. (Three of the 21 letters in the American 
Journal of Physics deal with this subject.) The 
most famous public argument along this line 
took place in 1925 between E. H. Kennard and 
W. S. Franklin.2® The former argues that the 
pressure term is certainly not the potential en- 
ergy of compression and cannot strictly be said 
to belong to the moving fluid; rather it is ‘‘trans- 
mitted through the fluid.” Franklin, however, 
insisted that potential energy is never localized 
in space, but ‘‘an idea which makes up for things 
ignored”’ (see H. Reichenbach’s”® introduction of 
the concept of interphenomenon to provide for 
just such abstractions as localization of potential 
energy). In this connection one is reminded of 
E. J. Routh (1876), who pointed out that so- 
called ignorable or cyclic dynamical coordinates 
result in an expression for the kinetic energy of 
spin that is formally equivalent to an additional 
potential energy. On this basis J. J. Thomson 
(1886) proposed to explain all potential energy 
by concealed motions, thus stamping an igno- 
ramus on any hope of locating it. It would seem 
that the usage of the term pressure energy en- 
joys a certain slang value, but that it is physically 
confusing owing to the actual doing Of external 
work by the moving fluid particle. (Pressure actu- 
ally plays the mathematical role of a Lagrangian 
multiplier in the derivation of Euler’s equations 
of motion from Hamilton’s principle.) Indeed, 
from a purist point of view, pressure itself is a 
misleading name in connection with energy in- 
asmuch as it is only a change in stress which is 
associated with work. 


28 E. H. Kennard, Science 62, 243 (1925); W. S. Franklin, 
Science 62, 397 (1925). 


29H. Reichenbach, Philosophic Foundations of Quantum 
Mechanics (California, 1946). 
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I personally cannot leave this topic of Ber- 
noulli’s theorem without urging some moderniza- 
tion in the use of illustrative applications of it. 
In an age of supersonic flight why not call 
attention to certain important features of super- 
sonic flow inherent in the following form of 
Bernoulli’s theorem for air (negligible viscous 
effects) flowing at a constant gravitational level : 


ip 
37° +— =B, 
2p 


y= 1.4. 


First of all, we note that whereas there is no 
limit to the value of the Mach number (the ratio 
of the flow speed to the acoustic speed ; its square 
is proportional to the ratio of the fluid’s available 
kinetic energy to its random thermal energy), 
there is an absolute limit to the flow speed of a 
fluid initially at rest (about 751 m/sec for air 
originally at room temperature). (An amusing 
paradox in the use of a concept like Mach number 
is that the effective Mach number of a projectile, 
which would have a Mach number of two if fired 


ee em, 


Fic. 5. Maxwell’s picture of Rankine’s conception of the 
formation of a shock wave (see reference 32). 


under room conditions, has a Mach number that 
is neither 0 nor four when fired into a supersonic 
stream flowing with a Mach number two.) 

Secondly, a flow speed equal to the acoustic 
speed is in a real sense a “critical speed,’’ inas- 
much as flows with speeds higher than this speed 
are supersonic and behave quite differently from 
subsonic flows. For example, if one considers the 
equation of continuity in the hydraulic form 
pgs is constant (s is the cross-sectional area), one 
can show that for a convergent nozzle a super- 
sonic flow decreases monotonously in speed but 
a subsonic flow increases; on the other hand, for 
a divergent nozzle a supersonic flow increases 
but a subsonic flow decreases. 

Supersonic flow, moreover, is associated with 
shock waves, which being waves of finite ampli- 
tude are an actual occurrence rather than the 
idealized (mathematical) acoustic waves of in- 
finitesimal amplitude. Although the basic equa- 
tions®® describing shock-wave phenomena were 


30 H. Hugoniot, J. Ecole Polyt. 57-59 (1887-89). 


SEEGER 


questioned by so eminent an authority as H. 
Lamb*! as late as 1932, their physical formation 
was understood by W. J. M. Rankine as early 
as 1869. Figure 5° shows how the greater speed 
of the more displaced fluid results in a steeper 
wave front. More recently, the study of shock- 
wave interactions* has revealed new physical 
phenomena. For example, shock waves not only 
differ quantitatively with respect to the acoustic 
laws of regular refraction and reflection [Fig. 
6a], but actually exhibit qualitatively dis- 


‘tinctive patterns such as Mach intersections 


[Fig. 6b]. It is believed*that such modern 
illustrations, supplemented with experiments in 
a shock tube* (the physicists’ wind tunnel), will 
definitely serve to make fluid dynamics a ‘“‘live”’ 
subject in physics. 


FLUID DYNAMICS IN THEORETICAL PHYSICS 


Let us now turn our attention to the treat- 
ment of fluid dynamics in a typical first-year 
graduate course such as the one commonly 
called ‘‘Theoretical Physics.” First of all, what is 
being done? 

Most of the fluid dynamics of such a graduate 
course is customarily concerned with a per- 
fectly inviscid incompressible fluid. One starts 
with the basic equations of continuity and of 
motion. Laplace’s partial differential equation 
follows from these if the fluid is assumed to be 
incompressible and irrotational. The major por- 
tion of this part of an American course is then 
devoted to the elegant mathematics associated 
with two-dimensional flows satisfying Laplace’s 
equation, namely, the velocity potential and the 
stream function, conjugate functions, conformal 
mapping, the method of images, etc. Kelvin’s 
circulation®® and Helmholtz’s vortices*® form a 
natural mathematical continuation. Applica- 
tions are usually restricted to the motion of a 

31H. Lamb, Hydrodynamics (Cambridge University 
Press, Cambridge, 1932), sixth edition. 

3 J. C. Maxwell, Theory of Heat (Longmans, Green and 
Company, London, 1871). 

83H. Polachek and R. J. Seeger, Proc. Symposia in 
Applied Math. (American Mathematical Society, 1949), 
Vol. 1, p. 119; J. Appl. Phys. 22, 640 (1951); Phys. Rev. 
84, 922 (1951). 

3% W. Bleakney and A. H. Taub, Revs. Modern Phys. 
21, 584 (1949). Bleakney, Weimer, and Fletcher, Rev. Sci. 
Instr. 28, 807 (1949). 

35 W. Thomson (Kelvin), Trans. Roy. Soc. (Edinburgh) 


25 (1869). 
36 H. von Helmholtz, Crelle’s Journal 55, 25 (1858). 
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Fic. 6a. White-light interferogram of regular shock- 
wave intersection—Mach number 2.48, Reynolds number 
10° (NOL Aerophysics Tunnel No. 3). 


sphere, with particular, attention to the virtual 
mass associated with unsteady conditions. Occa- 
sionally, there is a discussion of wave motion 
and of tides. It is quite disconcerting, however, to 
find little reference to Helmholtz’s discontinuous 
motions*’ and to modern nonlinear problems, in- 
volving the possibility of specific principles of 
superposition (to replace the general principle of 
linear superposition first used by D. Bernoulli). 
Gas dynamics, dealing with compressible fluids, 
is usually ignored. 

In the case of a perfectly viscous incompres- 
sible fluid one usually derives the Navier-Stokes 
equations of motion and then uses them to de- 
duce the Poiseuille law for a tube and the Stokes 
law for a sphere. In neither application, however, 
is there usually any discussion of the important 
phenomenon of turbulence or of the significant 
role of the Reynolds number. Inhomogeneous 
fluids*® (non-Newtonian) such as sols and gels 
are not even mentioned. 

It is not surprising that such superficial regard 
for phenomena per se has evoked severe criticism 
as exemplified by the hyperbolic expressions be- 


37H. von Helmholtz, Ber. Akad. Wiss., Berlin 215 
(1868). 


38 E. G. Richardson, Dynamics of Real Fluids (Edward 
Arnold and Company, New York, 1950). 


low. E. T. Bell in his sketch of L. Euler under the 
heading ‘‘Analysis Incarnate’’ says: ‘‘The physi- 
cal universe was an occasion for mathematics to 
Euler, scarcely a thing of much interest in itself; 
and if the universe failed to fit his analysis, it 
was the universe which was in error.’’”®° H. Rouse 
states: ‘‘During the past 150 years mathemati- 
cians and mathematical physicists have at- 
tempted completely rational solutions of the 
problems of fluid motion without any recourse 
to experimental measurements. These rational 
methods were made possible, however, only by 
such simplifying assumptions that the resulting 
solutions often bore little or no resemblance to 
the phenomena in question.’’*® R. von Mises 
claims in 1945: “There exists so far no mathe- 
matical theory that would account satisfactorily 
for the actually observed types of fluid motion.’’*° 
L. Prandtl and O. Tietjens remark: ‘‘A great 
amount of literature exists on the theoretical 
solution around moving bodies in the ideal fluid. 
However, we shall not dwell on these theories 
here, since they do not lead to a solution which 
conforms in the least with actually occurring flow 


Fic. 6b. Monochromatic interferogram of Mach shock- 
wave intersection—Mach number 2.48, Reynolds number 
105 (NOL Aerophysics Tunnel No. 3). 


89H. Rouse, Fluid Mechanics for Hydraulic Engineers 
(McGraw-Hill Book Company, Inc., New York, 1938). 

40R. von Mises, Theory of Flight (McGraw-Hill Book 
Company, Inc., New York, 1945). 
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phenomena.’”*! S. Goldstein, too, concludes: “In 
this book we are concerned with the motion of 
fluids whose viscosity is small; the phenomena 
in such motions sometimes approximately agree, 
sometimes violently disagree with the predic- 
tions of ideal fluid theory.---In these circum- 
stances it is surprising not that the results of 
theory and observation of disagree violently, but 
rather that they agree fairly well, as for the mo- 
tion of airship bodies.’’*? (Experimental observa- 
tions for fluids with small viscosity are well 
described by theory up to a Mach number of 
approximately 0.25 except for boundary layer 
and wake phenomena.) 

Without apologizing for the above criticisms 
we shall review for ourselves briefly the theo- 
retical motion of a particular body in a per- 
fectly inviscid fluid and in a perfectly viscous 
fluid. We choose a sphere for this purpose be- 
cause its symmetry aids in laboratory control 
and calibration and because its simplicity affords 
little opportunity for misunderstanding. The 
sphere, moreover, is of specific interest alike to 
mathematician, physicist, aerodynamicist, me- 
teorologist, and astronomer. In each case we 
shall look also at some observational data. Fi- 
nally, we shall mention some well-known ques- 
tions pertinent to the foundations of fluid 
dynamics. 


SPHERE IN A PERFECTLY INVISCID FLUID 


A practical problem of everyday concern is 
drag. Let us consider first a sphere moving at 
constant speed in a perfectly inviscid fluid. Note 
that the motion of the fluid with respect to the 
body is strictly steady motion, but that the mo- 
tion of the fluid with respect to a fixed point 
infinitely distant (body moving) is not steady. 
(This distinction may be of considerable signifi- 
cance for turbulent phenomena.) 

First of all, let us recall that, in general, there 
are three possible irrotational motions: 

(1) For simple potential (Dirichlet) motion of 
an inviscid, incompressible fluid (probably theo- 
retically true also for a compressible fluid), 

“L. Prandtl and O. G. Tietjens, Applied Hydro- and 


Aeromechanics (McGraw-Hill Book Company, Inc., New 
York, 1934). 


4S. Goldstein, Modern Developments in Fluid Dynamics 
> 32) (Oxford University Press, London and New York, 
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without acceleration and without free or solid 
surfaces, there is no force (but possibly a couple, 
if the body is not a sphere). This so-called para- 
dox of J. LeR. d’Alembert (or P. L. Dirichlet) is 
evident for a sphere in that no change in mo- 
mentum exists across an enveloping surface at a 
great distance away, inasmuch as the fluid mo- 
tion here can be regarded as the superposition 
of a uniform motion and that of a dipole or 
doublet source, for which the disturbing effect 
decreases with the cube of the distance. 

(2) Two-dimensional potential (N. Joukowski) 
motion with circulation only about a singular 
point® is characterized by a lift proportional to 
the density, the speed, and the circulation (cir- 
cuit line-integral of the velocity) itself, provided 
only that the Kutta-Joukowski condition of 
finite speed at the trailing edge is fulfilled. 

(3) In the two-dimensional discontinuous mo- 
tion for a flat plate“ there is a vortex sheet 
(observed by H. Benard and explained by T. von 
K4rman for a cylinder)** which arises from sheets 
of velocity discontinuity, although the motion as 
a whole is irrotational. It is noteworthy that 
apparently symmetrical causes here do not neces- 
sarily produce symmetrical consequences. 

For these several types of steady irrotational 
motion the following questions persist: which 
solution(s) occurs in nature? which solution is 


Fic. 7. Bubblegram (second photographs taken after a 
brief time interval) showing the velocity field about a 
sphere entering water (G. Birkhoff, Harvard). 


4 W. M. Kutta, Illust. Aeron. Mitteilungen 133 (1902). 
N. o Joukowski, Z. Flugtechnik Motorluftschiffahrt 1 
(1910). 

“ G. Kirchhoff, Crelle’s Journal 70, 289 (1869). 

45H. Bénard, Comptes rend. 147, 970 (1908). T. von 
Karm4n, Gott. Nachr. 509 (1911); 547 (1912). T. von 
Karman and H. L. Rubach, Physik. Z. 13, 49 (1912). 
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Fic. 8. Goldstein's comparison of theoretical values and 
experimental data for Stokes’ law (see reference 42) for 
a sphere. 


most stable? how does the introduction of vis- 
cosity modify the stability? 

Potential theory can be well utilized for 
streamlined bodies also in unsteady motion. As 
G. Green and G. G. Stokes** showed inde- 
pendently for an inviscid incompressible fluid 
there is virtually an increase in mass owing to the 
setting of the fluid in motion; for a sphere this 
so-called ‘induced mass’’ is } the mass of the 
displaced fluid. (In general, the virtual mass 
depends upon the direction of motion, as well as 
upon the body shape.) Applications of this 
seemingly abstract concept are rarely given. One 
example‘? is the motion of a spherical balloon 
when large differences in fluid density are in- 
volved. More striking, however, is the unsteady 
motion of a weightless sphere in a liquid.4* The 
sphere will have an acceleration, relative to an 
observer, three times that of the liquid itself. 
Another interesting example occurs in water 
entry (see Fig. 7), where the variation of cavity 
drag (dependent upon wake underpressure) may 
be interpreted as an ‘‘impact”’ induced mass due 
to the setting of the water in motion.*® 


SPHERE IN A PERFECTLY VISCOUS FLUID 


Let us now consider a perfectly viscous fluid 
for which, by definition, there is a constant co- 
efficient of tangential stress during motion and 


46G. Green, Math. Papers 316 (1833) (London, 1871). 
G. G. Stokes, Camb. Trans. 8 (1843); Math. and Phys. 
Papers (Cambridge, 1880-1905), Vol. 1, p. 17. 

47G. Birkhoff, Hydrodynamics (Princeton University 
Press, Princeton, 1950). 

48G. Birkhoff and T. E. Caywood, J. Appl. Phys. 20, 
646 (1949). 

49 A, May and J. C. Woodhull, J. Appl. Phys. 21, 1285 
(1950). ; 
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no slipping along the boundary. The equations 
of motion were formulated first by Navier.?7 The 
establishment of the boundary condition, how- 
ever, was early in doubt with respect to the 
very question of slipping. Stokes pointed out 
three possible hypotheses: slipping of fluid along 
fluid, slipping of molecules, no slipping. He de- 
cided in favor of no slipping at all.?” 

The equations of motion for a perfectly vis- 
cous fluid were solved by Stokes®® for the slow 
steady motion of a rigid sphere, without slipping, 
in an infinite homogeneous fluid. In the case of 
such “‘creeping’”’ motion, the viscous effects are 
more important than the inertial ones so that 
“‘semiquadratic”’ terms like (pu)du/dx can be 
neglected (u, the x velocity component). The re- 
sisting force turns out to be proportional to the 
speed ; it can be written 


pn Qaneyirstieatt 
= a =Onryau, 
pU(2a) 


where a is the radius of the sphere and U the 
undisturbed stream velocity (u/pU2a=R). 
In this connection one recalls Newton’s com- 
ment: ‘‘That the resistance of bodies is in the 
ratio of the velocity is more a mathematical 
hypothesis than a physical one.’’™! 

Evidently the viscous effects become less sig- 
nificant as one gets farther away from the body 
so that the inertial terms, in turn, then become 
relatively more important. For this case C. W. 
Oseen® obtained a better approximate solution 


TaBLE II. Stokes’ law applied to raindrops. 


Atmospheric con- 7 
densations of Diameter 
water (cm) 


Terminal 
velocity 


Reynolds 
(cm /sec) 


number 


Fog 0.002 
Mist A 1.67 
Drizzle Y 10 
Light rain 
Moderate rain 
Heavy rain 
Excessive rain 
Cloudburst 
(Maximum size) 


60 
267 
500 
840 

1400 
2667 


0.50 


50 G. G. Stokes, Trans. Cambridge Phil. Soc. 8, (1843). 

51]. Newton, Philosophiae Naturalis Principia Mathe- 
matica (1687), translation by F. Cajori (California, 1946), 
Book II, Prop. 35, Prob. 7, Cor. V. and Book II, Scholium 
to Prop. IV. 

8 C, W. Oseen, Arkiv. Math. Astron. Fysik 6, 29 (1910) ; 
Neuere Methoden und Ergebnisse in der Hydrodynamik 
(Leipzig, 1927). 
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for a sphere. Goldstein’s®:** more complete ex- 
pression for the drag coefficient is 


24, 3 19 
Co=—114—2-—_—_- 3 


71 
commas ss 
R 16 1280 


20480 


The two theoretical solutions for a sphere differ 
by 0.05 percent at Reynolds number 0.2, by 1.0 
percent at 1.0, and by 10 percent at 5.0; experi- 
mentally, the Stokes law can be considered valid 
for Reynolds numbers less than 0.05, and the 
Oseen law for Reynolds numbers less than 0.8 
(see Fig. 8). 

The failure of average texts on theoretical 
physics to specify generally experimental limita- 
tions is probably reflected in statements made in 
more elementary texts about the applicability 
of Stokes law. For example, one book states: 
“‘In the case of a sphere this resistance is pro- 
portional to the radius and to the velocity at 
moderate speeds. For this reason falling rain- 
drops reach a constant terminal velocity.’’ One 
immediately wonders what is meant by moderate 
speeds and what actually is the Reynolds-number 
range. With respect to the applicability of the 
Stokes law for the terminal velocity here, Table 
II, based upon W. J. Humphrey’s®™ data, speaks 
for itself. Humphreys himself warns: ‘‘The equa- 
tions of fall in the case of raindrops of average 
and large size: - -are little more than empirical.” 

One more experimental fact is worth noting 
here, namely, that for Reynolds numbers greater 
than 50, small air bubbles in a liquid rise in 
spherical spirals.*> For smaller Reynolds num- 
bers D. Riabouchinsky and J. Hadamard** 
showed in 1911 that for creeping motions the 
Stokes formula has to be modified theoretically 
for gas bubbles, inasmuch as the boundary con- 
dition at the surface of such a sphere specified 
continuity of velocity across it. Observations,’ 
however, indicate that it is the formula for a 
rigid sphere (with zero fluid velocity on its sur- 
face) which actually holds. 
ssa3} Goldstein, Proc. Roy. Soc. (London) ¥A123, 225 

54W. J. Humphreys, Physics of the Air (McGraw-Hill 
Book Company, Inc., New York, 1940). 

56H. Nisi and A. W. Porter, Phil. Mag. 46, 754 (1928); 
A. W. Porter, The Method of Dimensions (Methuen and 
Company, London, 1933). 

56D. Riabouchinsky, Bull. Acad. Cracovie 40 (1911); 
J. Hadamard, Compt. rend. 152, 1735 (1911). 
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It was early recognized by I. Newton® that 
for bodies traveling at higher speeds drag is 
proportional to the square of the speed rather 
than to the speed itself. For example, in the case 
of the inelastic impact of fluid particles upon a 
flat plate transverse to the flow, both the mass 
impinging per unit time and the change in speed 
itself are proportional to the speed. Only the 
front of the plate, of course, is significant in this 
instance. In 1907, A. G. Eiffel showed that the 
resistance of two flat plates, together but in- 
clined to the flow, is only 60 percent as much as 
when they are separated. On the Newtonian 
picture separation of the plates should have 
made no difference; accordingly the rear of a 
plate must have some role to play. In other 
words, the Newtonian view is adequate for a 
rarefied medium where the particles are inde- 
pendent, but not for a continuum where the 
particles are interdependent. 

Another general problem of considerable prac- 
tical importance is the existence of similar flows, 
which permit the scaling of information obtained 
experimentally with models. Similar inviscid 
flows are determined for incompressible fluids by 
what may be called the Euler number (ratio of 
the pressure to the so-called dynamic pressure), 
and for compressible fluids by the Mach number. 
(The Euler number can be regarded as being 
inversely proportional to the square of the Mach 
number for comparative purposes.) An addi- 
tional parameter, however, is needed for specify- 
ing similar viscous flows, namely, the Reynolds 
number, which measures the inertial influence 
relative to the viscous effect. The kinematic 
viscosity is important here. For example, air is 
effectively 15 times more viscous than water 
(with respect to kinematic viscosity), although 
the coefficient of viscosity for air is only 1 per- 
cent that for water. Now large kinematic vis- 
cosity (low Reynolds number) itself may be said 
to favor laminar flow inasmuch as it smooths out 
transversal components of motion, whereas rela- 
tively large inertia (high Reynolds number) 
favors turbulence by continuing transverse com- 
ponents. Thus, an increase in the coefficient of 
viscosity will increase the speed at which the 
critical Reynolds number occurs, ceteris paribus, 
but an increase in density will decrease this 
speed. 
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Viscous effects, in general, are confined to a 
small region about the body. In this boundary 
layer®’ even a small coefficient of viscosity may 
result in large forces, inasmuch as the latter de- 
pend also upon the velocity gradient. If the 
boundary layer is sucked away with little ex- 
penditure of energy, the viscous flow closely 
approximates the potential flows with and with- 
out circulation. “In general, the agreement be- 
tween theory and experimental results becomes 
better when viscosity becomes smaller.’’4! Much 
more significant, however, is the indirect effect 
that even small viscosity may have in selecting 
stable flow from among the possible types of flow. 

Paradoxical descriptions of phenomena have 
colored the differing outlook of G. Birkhoff in his 
stimulating book, Hydrodynamics.” ‘“‘No treat- 
ment of slightly viscous fluids which is based on 
this theory [classical asymptotic theory of par- 
tial differential equations] can come close to 
describing the real behavior of such fluids. It 
will never predict a turbulent wake, a vortex 
sheet, or transition at R,.’’ A small change in 
the coefficient of viscosity, for example, from 
0.5X10-5 to 110-5, doubling the Reynolds 
number, may apparently change the type of 
flow. completely, as from laminar to turbulent 
flow. He follows Oseen® in stressing the fact 
that the limit of a solution of a partial differential 
equation for coefficients of the highest order 
terms approaching zero is not necessarily the 
same as the solution of the equation with these 
coefficients already zero. The truth of this state- 
ment is evident in that second-order physical 
effects would otherwise be always negligible. 
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Fic. 9. Sphere drag coefficient as a function of 
Reynolds number—subsonic speeds. 


57 L, Prandtl, Proc. Math. Congress, Heidelberg (1904). 


Fic. 10a. Shadowgram of flow about a 37-mm sphere 
with separation in a laminar boundary layer—Mach 
number 0.445, Reynolds number 3.68 X 105 (A. C. Charters, 
BRL free-flight range). 


Moreover, no matter how infinitesimal the co- 
efficient of viscosity may become, finite slip flow 
will certainly not develop in a fluid initially at 
rest along a boundary. In any event more rigor- 
ously thoughtful experimentation is needed be- 
fore the role of viscosity will be known precisely. 

In the light of the above complexity the fol- 
lowing naive conclusion from a popular (1934) 
textbook on theoretical physics exemplifies strik- 
ingly how inadequately critical and somewhat 
meaningless our teaching may be: “The dis- 
crepancy [with respect to a circular cylinder] is 
to be attributed to the fact that even for fluids 
of low viscosity the friction at the surface is not 
negligible.” 


SPHERE IN A REAL FLUID 


We shall now consider briefly some particular 
phenomena of fluids flowing about spheres, first 
of all, as observed in wind tunnels. G. Eiffel5* 
found experimentally in 1911 that the drag co- 
efficient of a sphere was 0.18 as measured in the 
French tunnel at Auteuil. A year later L. Féppl 
obtained a value of 0.44 in the German tunnel at 
Gdttingen. It is amusing to find Eiffel explaining 
away the discrepancy as being due to different 
critical velocities for the flows, only to have 
Rayleigh explaining away the latter in view of 


58G. Eiffel, Nouvelles Recherches sur la Resistance de 
l’Air et l’ Aviation (Dounod and Rivat, 1914). 
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Fic. 10b. Shadowgram of flow about a 37-mm sphere 
with a separation in a turbulent boundary layer—Mach 
number 0.519, Reynolds number 4.29 X 105 (A. C. Charters, 
BRL free flight-range). 


the compensating effect of different diameters so 
that the Reynolds number was unchanged. The 
real explanation was given by L. Prandtl®® in 
1914 in terms of differences in the turbulence 
levels of the tunnels. For example, the critical 
Reynolds number varies from 100,000 in a very 


turbulent tunnel to 368,000 in a low-turbulence 
tunnel; in the atmosphere it is about 385,000. 
The phenomenon associated with these early 
experiments is sometimes called the Eiffel para- 
dox, namely, that near the critical Reynolds 


Fic. 11a. Schlieren photograph of flow about a 10-cm 
sphere—Mach number 3.24, Reynolds number 710° 
(R. Lehnert, NOL Aeroballistics Tunnel No. 2). 


°L. Prandtl, Nachr. Ges. Wiss. Géttingen, Math.- 
physik KI. 177 (1914). 
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Fic. 11b. Schlieren photograph of flow about a 10-cm 
sphere—Mach number 4.24, Reynolds number 4X105 
(R. Lehnert, NOL Aeroballistics Tunnel No. 2). 


number the resistance of a sphere actually de- 
creases as its speed increases owing to a large 
decrease in the drag coefficient (see Fig. 9). The 
decrease is due to the separation of the fluid 
from the body farther back in the boundary 
layer for turbulent than for laminar flow there, 
with a consequently greater neutralizing (back) 
pressure at the base [see Figs. 10a and b]. 
(Turbulence mixing in the stream thickens the 
boundary layer; although the flow-transition 
point in the boundary layer may move forward, 
the outer fluid effectively pulls the boundary fluid 
farther along the surface to the separation point.) 

At supersonic speeds there is a shockwave at 
the head [see Figs. 11a, b, c], which ac- 


Fic. 11c. Schlieren photograph of flow about a 20-cm 
sphere—Mach number 5.18, Reynolds number 6X105 
(R. Lehnert, NOL Aeroballistics Tunnel No. 2). 
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Fic. 12a. Sphere drag coefficient as a function of 
Reynolds number—supersonic speeds. 
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counts for more than 90 percent of the total 
drag of a sphere, so that there is little influence 
resulting from changing Reynolds number [see 
Figs. 12a and bj]. An unusual headwave 
formed by the faster moving muzzle gas on the 
back of a projectile is seen in Fig. 13. An inter- 
ferogram of a sphere in a supersonic jet is shown 
in Fig. 14. The supersonic flow about a sphere 
in a low density fluid with slip flow is shown by 
the nitrogen after-glow technique. 

Ballistics ranges®® afford a complementary 
means (see Fig. 16) of investigating supersonic 
phenomena in that the physical conditions of 
pressure, temperature, and turbulence in free 
flight are, in general, quite different from those 
in a wind tunnel. Typical wave phenomena, first 
photographed with schlieren technique by E. 
Mach,® are evident in the shadowgrams of Figs. 
17a, b, c, for different pressures in the 
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Fic. 12b. Sphere drag coefficient as a function of Reynolds 
number R; computed behind a normal shock wave. 


60 A. C. Charters and R. N. Thomas, J. Aeronaut. Sci. 
12, 468 (1945). 


61 EF. Mach, Wien. Ber. 98, 1310 (1889). 
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Fic. 13. Headwave on rear of projectile due to faster 
moving muzzle gases (W. Bleakney, Palmer Physical 
Laboratory). 


Naval Ordnance Laboratory Pressurized Ballis- 
tics Range. The wavelets behind the headwave 
are an indication that supersonic flow exists 
there. A shock wave occurs also at the separa- 
tion point, where there is a change of flow direc- 
tion. At the rear of the body supersonic vortices 
produce wavelets which coalesce in a shock tail- 
wave, where the vortices merge into large ir- 
regular ones. A comparison photograph of a 
sphere at Mach number 16.7 in monatomic kryp- 
ton is given in Fig. 18. For diatomic gases a new 
phenomenon becomes increasingly important as 
the temperature increases, namely, molecular 
dissociation. For example, bromine dissociation 
is found to have occurred already at a Mach 
number of 13.5 (observed at the Naval Ordnance 
Laboratory). Such preliminary results in the new 


Fic. 14. Interferogram of the flow of an axially sym- 
metric jet about a -in. diameter sphere—Mach number 
1.7 (R. Ladenburg, Palmer Physical Laboratory). 
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Fic. 15. _ Nitrogen afterglow photograph of a shock wave 
for a 0.7 in. sphere in a low density, supersonic flow— 
Mach number 2.6, Reynolds number 760/inch, P.:=110 
microns (S. A. Schaaf, University of California Low 
Density Wind Tunnel). 


field of aeroballistic research have been described 
elsewhere by the author in an article entitled, 
“On Aerophysics Research.” 

All the above complex phenomena observed 
in free flight and in wind tunnels give credence 
to the statements in 1934: “It is a deplorable 
fact that no theory of drag yet exists which even 
approximately does justice to the experimental 
results.” Although there have been subsequent 
improvements, the explanation of the drag of a 
sphere is still an enigma. 
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Fic. 16. Comparison of drag-coefficient data from free 
flight and from wind tunnel for sphere. 


® R. J. Seeger, Am. J. Phys. 9, 459 (1951). 
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Fic. 17a. Shadowgrams of 3-in. (}-in.) sphere at Mach 
number 2.8, Reynolds number 6150 (4200), and pressure 


0.0097 atmos (0.0094) (A. R. Witt, NOL Pressurized 
Ballistics Range No: 3). 


FOUNDATIONS OF FLUID DYNAMICS 


What is more, many questions persist as to 
the very foundations of fluid dynamics, for ex- 
ample, with respect to even the equations of 
motion. In continuum mechanics the conserva- 
tion equations (five in number) contain time 
derivatives of two thermodynamic variables (say, 
pressure and density) and of three kinematic 
variables (velocity components). In addition, the 
equations contain components of a stress tensor 


Fic. 17b. Shadowgram of 4-in. sphere at Mach number 
2.0, Reynolds number 1X 108, and pressure 2 atmos (A. R 
Witt, NOL Pressurized Ballistics Range No. 3) 
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Fic. 17c. Shadowgram of 3-in. sphere at Mach number 
1.2, Reynolds number 0.9 X 108, and pressure 4 atmos (A. R. 
Witt, NOL Pressurized Ballistics Range No. 3). 


and of a heat-flow vector. Hence, special assump- 
tions about these quantities have to be added if 
the equations are to form a determined system. 
A simple theoretical assumption is the vanishing 
of both the stress tensor and the heat-flow vector. 
In this case, the equations of motion become 
those of Euler for a perfectly inviscid fluid, and 
the energy equation specifies isentropy along a 
streamline. A simple phenomenological assump- 
tion, on the other hand, is that the stress tensor 
should be proportional to the velocity gradient, 
and the heat-flow vector to the temperature 
gradient. Under the further stipulation of tensor 
invariance these conditions lead to the Navier- 
Stokes equations for a perfectly viscous fluid. 
The viscosity appearing here may be formally 
classified into the usual shear viscosity and a so- 
called dilatational viscosity. The second viscosity 
coefficient is found to be equal to — (2/3) if the 
average principal stress terms are chosen to 
agree with the mean static pressure for com- 
pressible fluids (this stipulation need not be 
made). As early as 1851 Stokes recognized the 
lack of a rigorous justification for this hypothesis. 

The only theoretical argument is given by the 
kinetic theory of a dilute monatomic gas.“ Yet 
the derivation of the macroscopic behavior of a 
gas from its molecular properties is still lacking 
in rigorous treatment. For example, specifying 
the initial values of the above five variables 

6 PD. Enskog, Dissertation (Uppsala, 1917). Y. Rocard, 


L’hydrodynamique et la théorie cinétique des gaz (Gauthier- 
Villars, 1932). : 


of the hydrodynamic equations yields only an 
approximate solution (not unique) to the funda- 
merital integro-differential equation for the mo- 
lecular distribution function. Other approxima- 
tions can generally be obtained by assigning 
initial values to various moments of the distribu- 
tion function with respect to the molecular ve- 
locity (e.g., in the neighborhood of equilibrium 
the five-moment approximation of D. Hilbert, 
and the thirteen-moment extension by H. Grad). 

Experimental data are sparse. N. Lieberman’s® 
recent measurements of the so-called second co- 
efficient of viscosity apparently explain the 
acoustic absorption of high frequency waves in 
liquids. Nevertheless, the fundamental question 
remains: does the mean pressure depend upon 
the rate of expansion? Perhaps, the expansion- 
rate-measuring-apparatus developed for the Na- 
val Ordnance Laboratory by A. M. F. J. Michels 
will be useful for investigating such phenomena.® 
(It has already been used to show experimentally 
the constancy of the entropy of a rapidly ex- 
panding real gas.) 

Turbulence, of course, is still a perennial 
puzzle. Most success has been achieved to date 
by regarding it phenomenologically, with as- 
sumed ignorance as to the underlying equations 
of motion. A basic question here is whether or 
not quasi-steady turbulence is a second solution 
of the Navier-Stokes equation.®” 


Fic. 18. Shadowgram of }-in. sphere shot into krypton 
at a pressure of 30.5 mm of mercury and at a temperature 
of 103°K—Mach number 16.7 (C. R. Donaldson, NACA 
LAL free-flight range). 


64H. Grad, Lecture Notes on Kinetic Theory and Sta- 
tistical Mechanics (New York University, 1950). 

65 L. N. Lieberman, Phys. Rev. 75, 1415 (1949). 

66 A. Siegel, Doctoral dissertation, University of Amster- 
dam, 1952. 

67 M. Planck, ‘‘The Mechanics of Deformable Bodies’”’ 
(Introduction to Theoretical Physics 11), (Macmillan Com- 
pany, New York, 1932). 
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Even steady-flow phenomena, indeed, are sub- 
ject to critical analysis. In view of the importance 
of time variations in determining stability of 
flow, the following uncomfortable question lurks: 
does strictly steady flow actually occur in nature? 

Finally, the precise role of heat in fluid dy- 
namics opens up a whole new field of questions. 
In this connection one is reminded of Mach’s 
remarks: ‘‘Purely mechanical phenomena do not 
exist.---Questions concerning the motion of 
gases cannot be dealt with as questions of pure 
mechanics, but always involve questions of heat. 
[Ed. note: Nor can even a thermodynamical 
treatment always suffice; it is sometimes neces- 
sary to go back to the consideration of molecular 
motions. ]’"!® Rapid changes in gases, of course, 
are associated with temperature changes and 
expansive forces. Much more startling is P. 
Duhem’s challenging observation: ‘The sup- 
position of an incompressible viscous fluid is 
thermodynamically incompatible with the con- 
cept of a fluid.’’47:®8 As flow speeds about Mach 
number 2.5 are attained, a heat similarity pa- 
rameter, the Péclet number (measures the rela- 
tion of convective heat to conductive heat), has 
to be considered in addition to the Mach number 
for compressible flow and to the Reynolds num- 
ber for viscous flow. An unanswered question 
here is whether these parameters can be regarded 
as being independent or whether they become 
involved in a functional relation to describe phe- 
nomena at hypersonic (Mach number greater 
than five) speeds.® 


TRADITION AND BETRAYAL 


In spite of all the above interesting physics, 
what is the traditional attitude of physicists to 
fluid dynamics? Perhaps, two statements will 
suffice. P. Lenard says with respect to the hydro- 
dynamic equations: ‘“‘They nevertheless contain 
only old knowledge, nothing fundamentally new 
is added, or made evident by their treatment.’’® 

68 P. Duhem, Traite d’energetique 2 (Paris, 1930). 


69 P. Lenard, Great Men of Science (Macmillan Company, 
New York, 1933). 
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In M. von Laue’s recent History of Physics,” 
in which the Navier-Stokes equations are not 
even mentioned, one reads: ‘“‘Nobody expects 
these results to go beyond the Newtonian founda- 
tions.”” True—no poem contains new letters, 
possibly not even unorthodox words—but fresh 
ideas or possibly even fresh expressions! Pity 
the poor student who is not brought into con- 
tact with the existence of new phenomena in fluid 
dynamics and the need for new concepts in 
its physical description. Lest we teachers forget: 
Tradition does mean basically a handing over; 
but it also comes from the same etymological 
source as betrayal. The statement of the mathe- 
matician G. Birkhoff must be provocative to us 
physicists: ‘‘It seems surprising that this subject 
[fluid dynamics], to which so much was con- 
tributed by Stokes, Helmholtz, Rayleigh, and 
other physicists of the nineteenth century, should 
owe all its more recent progress to mathemati- 
cians and engineers. I feel that a renewed interest 
by physicists in the accurate correlation of theory 
with experiment, would help greatly in providing 
it with a more solid scientific basis.’’*” 

In summary, if the teaching of fluid dynamics 
in physics involves essentially the understanding 
of concepts on the part of the student so that he 
grows intellectually, and if the teaching of fluid 
dynamics in physics is primarily a matter of 
selecting the phenomena and their related con- 
cepts so that culture develops out of the quality 
of the overtones, then both in the general physics 
course of the beginning college student and in the 
theoretical physics course of the beginning gradu- 
ate student Shakespeare’s words (with respect 
to the ghost’s appearance at the end of Act I) 
are pertinent: 


Horatio: ‘‘O day and night, but this is wondrous 
strange. 

Hamlet: ‘‘And therefore asa strarger give it welcome. 

There are more things in heaven and earth, 
Horatio, 


Than are dreamt of in your philosophy.” 


70 M. von Laue, History of Physics (Academic Press, Inc., 
New York, 1950). 
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The possible motions of a sphere suspended on a string are found to consist of combinations 


of (1) a pendulum type of motion and (2) a rolling motion, in which the ball rolls in angular 
harmonic motion about an axis very near its center. For (1) it is shown (a) the usual formula 
for the length of the equivalent simple pendulum should be replaced by L(1+272/5L?+-41r3/25L8 
+--+-), where r is the radius of the ball and L the distance of point of suspension to center of 
the sphere, and (b) the best way of starting this motion is to pull the ball aside by a horizontal 
force through its center. For (2) it is shown that the period of oscillation is very close to 


2av (27/5g). 


Photographs are given showing experimental verifications of these effects. 





I. INTRODUCTION 


T is frequently stated that for a sphere of 
radius r on the end of a string of length 
(L—r), swinging as a pendulum, the length of 
the equivalent simple pendulum is L(1+3r?/L?). 
It is not generally realized that this is an approxi- 
mation. It is accurate for a sphere on the end of 
a massless rigid rod, where the couple to cause 
rotation of the sphere is provided by the rigidity 
of the rod; but for the sphere on a string this 
couple must be provided by the sphere’s weight 
and the string tension, and consequently the 
line of the string does not, in general, pass 
through the center of the sphere. 

As the sphere is free to rotate about the point 
of contact with the string, the system has a de- 
gree of freedom more than is assumed for the 
simple pendulum, and this permits motions other 
than the simple pendulum motion. The various 
possible motions are discussed below. 






II. BASIC THEORY 


The problem of the double pendulum has been 
treated by Timoshenko,! as an example of the 
Lagrange method using generalized coordinates. 
In the following the problem is treated by the 
more direct force methods, as a preliminary to 
obtaining further results not given explicitly by 
Timoshenko. The usual assumption is made that 
all angles are kept sufficiently small for their 
sines to be replaced by the angles, and their co- 
sines equated to unity. 

* Now at Southwest Research Institute, San Antonio 6, 


Texas. 


1S. Timoshenko, Vibration Problems in ataeniniliits 


(Van Nostrand Company, Inc., 1941), p. 203. 
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Consider the system, when swinging, in the 
position shown in Fig. 1. The coordinates of Q are 


x=(L—r)6+r(6+¢) =L6+r¢, 


(1) 
y=L. 
Resolving vertically and horizontally, and using 
Eq. (1), one obtains 


) 


ae dod 


L—+r—+ g0=0. (2) 
de dt 


Taking moments about Q and writing mbr? for 
the moment of inertia, where } is a numerical 
constant (=? for a sphere), gives 


a6 





: (3) 
dt? 


Subtracting Eq. (3) from Eq. (2) leaves 


2 


g 
(L—r)—+g6—-o=0. 
dt? b 


(4) 


Eliminating ¢@ between Eqs. (4) and (2) yields a 
fourth-order differential equation in 6. As we 
expect harmonic oscillations, we look for solu- 
tions of the form 0=A sin(wi+y). Substitution 
of this in the differential equation gives 


br(L—r)w*—g(L+br)w?+ g?=0 (5) 


This is a quadratic in w, with two positive real 
roots, w; and we. It is convenient to treat the two 
modes, involving w; and we respectively, sepa- 
rately. Defining a dimensionless constant a by 
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Fic. 1. Coordinate 
system for the freely 
swinging pendulum. 


a=r/L, the solution of Eq. (5) can be written 
g 1+bea 4ba(1i—a)\? 
we {14(1-—-—=) |. 6) 
2baL 1—a (1+)ba)? 
Since r<L, a<1, Eq. (6) can therefore be expan- 


ded as a power series in a, giving 
positive sign: w;?=(g/br){1+a+(1+))a? 
+(1+6+0?)+ eee a 
negative sign: we? =(g/L)(1—ba? 
—Bai+---), 


(7) 


(8) 


To obtain the relation between @ and , Eq. (4) 
can be rewritten 


d’6 
dt? 


bL 
= b6+—(1—a) 
g 


=| 1 ~~ —a)et|, 


Ill. THE LOW FREQUENCY SOLUTION 


Equation (8) gives the period of the pendulum 
type of motion as T = 27/we. Since the equivalent 
simple pendulum has a length L, given by T 
=2n/(L./g), inverting and expanding Eq. (8) 
gives 


L,=L(1+be?+b%at+ - ++). (10) 


This equation shows that the usually quoted 
formula for the length of the simple pendulum 
equivalent to a ball on a string is only an approxi- 
mation, being in fact the first two terms of Eq. 


(10). The approximation is of course a very good 
one if L>r, as is usually the case. 
Equation (9), with we from Eq. (8), yields 


$= bab{1+ba—b(1—b)a?+---}. (11) 


Thus, if L>r, a1 and ¢<86, the radius PQ 
does not move much from the line of the string. 
This justifies the opening sentence of this section, 
in which it is stated that this is “‘pendulum”’ 
type of motion. 

If the swinging is started at random, it is 
probable that a component involving w, will be 
introduced. This is avoided if the angles satisfy 
Eq. (11) at the moment of release. Suppose the 
pendulum is drawn aside by a horizontal force, 
as in Fig. 2. 

Resolving vertically and horizontally and tak- 
ing moments round P, one obtains 


f=ro/6. (12) 


To satisfy Eq. (11), retaining only the first 
power of a, therefore requires f=rba. Thus, the 
horizontal force should be applied below the 
center of the ball by an amount 2r?/5L. For 
L>r, a very good approximation to this will be 
produced by pulling the ball aside by a hori- 
zontal force through its center. If the force were 
applied at P, or at the bottom of the ball, then 
f=+r, and from Eq. (12) this would give 
= +06. Consequently both these methods would 
introduce a large amplitude of the unwanted 
component. 


Fic. 2. System of 
forces for pendulum 
pulled aside, before 
releasing. 
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IV. THE HIGH FREQUENCY SOLUTION 
Equation (7) gives the period of oscillation as 
br 3 2r\3 
T=2r|—<1 —a+:: | ~2r(=-) j 
g Sg 
Equation (9) gives 
¢=b0(1—L/br) = —0(1—ba)/a. (13) 


It is seen that ¢ is opposite in sign to 0, and 
larger. The ball ‘‘rolls,”’ with relatively little 
angular motion of the string. To determine the 
center of rotation, consider a point whose rest 
coordinates on Fig. 1 are (0, L+ Y), so that the 
point is Y below the center of the sphere. When 
swinging, its x coordinate will be 


x=(L—r)0+(r+ Y)(6+¢) =L0+aL¢+ Y(6+¢). 
Thus, x=0, if 
Y= —L(1+a@/6)/(1+¢/6). 


Substitution from Eq. (13), retaining only the 


Fic. 3. Double exposure photograph of 
simple pendulum motion. 


Fic. 4. Triple exposure photograph of bob undergoing 
rotational oscillations. 


lowest power of a, yields 
Y=oLb=br(r/L). 


Thus, the ball will roll about a point only slightly 
below the center of the sphere. 


V. EXPERIMENTAL VERIFICATION 


For experimental work, a flat disk (b=0.5) 
was used, since the effects should be more pro- 
nounced with the larger 6. The disk was painted 
white, and the diameter PQ was marked in 
black. The disk was 9 inches in diameter. 

Results of Sec. III (pendulum motion). For 
this system a string 3 inches long was used. 
Thus, r = 43 inches; L=7} inches; a=0.6; b=0.5. 

The system was set swinging and photo- 
graphed, two exposures of 1/100 sec being made 
on the same film. The photograph is reproduced 
in Fig. 3. Application of a straight edge will 
show that PQ is not in line with the string and 
that the angle ¢ has the same sign as 6. The 
angle between the two string positions was 
measured and called 6’; and the angle between 
the two diameters was called (6’+¢’). The meas- 
urements gave 6’ = 13.7° and 6’+¢’ = 18.1°; thus, 
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experimentally, ¢/8=0.32. From Eq. (11), theo- 
retically ¢/6=0.36. The agreement between ex- 
periment and theory is therefore reasonably good. 

An attempt was made to check the period. 
The usual formula for the equivalent simple 
pendulum gives T=0.948 sec, while Eq. (10) 
gives 0.975 sec. Owing to difficulty in keeping 
the pendulum swinging in one plane, the experi- 
mentally determined periods varied from 0.960 
to 0.970 sec. The result is therefore not con- 
clusive, but tends to favor the theory given here. 

Results of Sec. IV (rolling motion). For this a 
string 9 inches long was used. Thus, r=43 
inches; L = 134 inches; a=}; b=0.5. 

Three exposures, at different instants of the 
motion, were taken on the same film, and the 
photograph is reproduced in Fig. 4. It was 
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difficult to start the rolling motion without some 
component of the pendulum motion also occur- 
ring, which accounts for the fact that the three 
diameters do not intersect at a point. Never- 
theless, it can be seen that the center of rotation 
is on the disk, slightly below the center of the 
disk. 

The period was measured as 0.36 to 0.38 sec 
(counting such fast oscillations was difficult). 
Equation (7) gives a period of 0.384 sec. The 
agreement is as good as can be expected. 
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The Tape Recordings of Important Speeches 


E. U. Connon, K. K. Darrow, E. Fermi, J. C. SLATER 


Free tape recordings are available on loan to members 
of the Association who wish to use them for presentation 
to classes, seminars, or science’ clubs. The subjects are 
four of the six lectures presented as a Symposium on 
Physics Today at the Twentieth Anniversary meeting of 
the American Institute of Physics, Chicago, Illinois, on 
October 25, 1951, as follows: 


The atom. E. U. Conpon. 

Physics as science and art. K. K. Darrow. 
The nucleus. E. FERMI. 

The solid state. J. C. SLATER. 


Professor Slater’s speech has been added to the list 
since our last announcement. 


Recordings are now available in two forms: 


1. Double track. Each speech is complete on one reel, 
which must be played at 3} inches per sec. No 
rewinding is necessary. Track No. 1 is played to 
the end, then the reels are reversed (and turned 
over) to play track No. 2. This procedure rewinds 
the tape into its original condition. 


2. Single track. Two 7-inch reels are required for 
each address. They must be played at 7} inches 
per sec, and, of course, rewound to restore them 
to their original condition. 


Requests for recordings should be made to THomas H. 
Oscoop, Editor, American Journal of Physics, Michigan 
State College, East Lansing, Michigan. To allow some 
flexibility in scheduling, a request should specify a pre- 
ferred date and two alternative acceptable dates, as well 
as the type of recording (single or double track, or either). 

No charge is made for the loan of recordings. The bor- 
rower is expected, however, to pay the return postage to 
the Editor’s office or to any other address the Editor 
specifies. Transportation is to be by first class mail, special 
delivery, which will cost less than one dollar per reel at 
present rates. Reels will be mailed in convenient reusable 
boxes. 

Prospective borrowers may consult the published papers 
in Physics Today, November, 1951, January, 1952, and 
March, 1952. 





Reproductions of Prints, Drawings, and Paintings of 
Interest in the History of Physics 


50. Fluents and Fluxions 


E. C. Watson 
California Institute of Technology, Pasadena 4, California 


(Received April 21, 1952) 


The frontispiece of John Colson’s The Method of Fluxions and Infinite Series together with 
explanatory discussion is reproduced. It is an interesting attempt to explain the method of the 


calculus pictorially. 





HE method of fluxions (which we now call 

the infinitesimal calculus) apparently oc- 
curred to Newton while he was still an under- 
graduate at Cambridge (he took his B.A. degree 
in January, 1665, shortly after his twenty-second 
birthday). The Methodus Fluxionum was begun 
in 1664 and finished in 1671, but remained un- 
published until 9 years after his death, when in 
1736 it appeared in an English translation by 
Joun Cotson as The Method of Fluxions and 
Infinite Series. The Latin original did not appear 
until 1779, when it was published in the first 
volume of the Opera Omnia. In the meantime, 
however, the fundamental principles of the 
method had been published in Lemma II in the 
second book of the Principia (1687), in the 
Opticks of 1704, and elsewhere. 

Cotson’s translation of 1736 is of especial 
interest from the point of view of this series of 
reproductions because of the frontispiece. This 
engraving is here reproduced and shows, as an 
illustration of the method of fluxions, a hunter 
attempting to shoot down two birds at once. 

The full title of this important work reads as 
follows: The Method of Fluxions and Infinite 
Series; with its Application to the Geometry of 
Curve-Lines. By the Inventor, Sir Isaac Newton, 
Kt, Late President of the Royal Society. Trans- 
lated from the Author's Latin Original not yet 
made publick. To which is subjoin’d, A Perpetual 
Comment upon the whole Work, Consisting of 
Annotations, Illustrations, and Supplements, In 
order to make this Treatise A compleat Institution 
for the use of Learners. By John Colson, M.A. 
and F.R.S., Master of Sir Joseph Williamson’s 
free Mathematical-School at Rochester. London: 
Printed by Henry Woodfall; and Sold by John 
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Nourse, at the Lamb without Temple-Bar. 
MDCCXXXVI. 

CoLson’s translation of Newton’s definitions 
of the terms ‘‘Fluents” and ‘‘Fluxions”’ reads as 
follows: 

“‘Now those Quantities which I consider as 
gradually and indefinitely increasing, I shall 
hereafter call Fluents, or Flowing Quantities, and 
shall represent them by the final Letters of the 
Alphabet v, x, y, and z; that I may distinguish 
them from other Quantities, which in Equations 
are to be consider’d as known and determinate, 
and which therefore are represented by the initial 
Letters a, b, c, etc. And the Velocities by which 
every Fluent is increased by its generating 
Motion, (which I may call Fluxions, or simply 
Velocities or Celerities,) I shall represent by the 
same Letters pointed thus é, #, y, and z. That is, 
for the Celerity of the Quantity v I shall put 3, 
and so for the Celerities of the other Quantities 
x, y, and z, I shall put 4, y, and z respectively.” 


The frontispiece is explained as follows: 


“But lastly, I believe it may not be difficult 
to give a pretty good notion of Fluents and 
Fluxions, even to such Persons as are not much 
versed in Mathematical Speculations, if they are 
willing to be inform’d, and have but a tolerable 
readiness of apprehension. This I shall here 
attempt to perform, in a familiar way, by the 
instance of a Fowler, who is aiming to shoot 
two Birds at once, as is represented in the 
Frontispiece. Let us suppose the right Line AB 
to be parallel to the Horizon, or level with the 
Ground, in which a Bird is now flying at G, 
which was lately at F, and a little before at E. 
And let this Bird be conceived to fly, not with 
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Ta yore Jaras, 7H awa yowss. 


Fic. 1. A hunter attempting to shoot down two birds 
at once, as an illustration of the method of fluxions. 
(From John Colson’s English translation of Newton’s 
Methodus Fluxionum, 1736.) 


an equable or uniform swiftness, but with a 
swiftness that always increases, (or with a Ve- 
locity that is continually accelerated,) according 
to some known rate. Let there also be another 
right Line CD, parallel to the former, at the 
same or any other convenient distance from the 
Ground, in which another Bird is now flying at 
K, which was lately at J, and a little before at 
H; just at the same points of time as the first 
Bird was at G, F, E, respectively. But to fix 
our Ideas, and to make our Conceptions the 
more simple and easy, let us imagine this second 
Bird to fly equably, or always to describe equal 
parts of the Line CD in equal times. Then may 
the equable Velocity of this Bird be used as a 
known measure, or standard, to which we may 
always compare the inequable Velocity of the 
first Bird. Let us now suppose the right Line EH 
to be drawn, and continued to the point L, so 
that the proportion (or ratio) of the two Lines 
EL and HL may be the same as that of the 


Velocities of the two Birds, when they were at 
E and G respectively. And let us farther sup- 
pose, that the Eye of a Fowler was at the same 
time at the point L, and that he directed his 
Gun, or Fowling-piece, according to the right 
Line LHE, in hopes to shoot both the Birds at 
once. But not thinking himself then to be 
sufficiently near, he forbears to discharge his 
Piece, but still pointing it at the two Birds, he 
continually advances towards them according to 
the direction of his Piece, till his Eye is presently 
at M, and the Birds at the same time in F and J, 
in the same right Line FJM. And not being yet 
near enough, we may suppose him to advance 
farther in the same manner, his Piece being 
always directed or level’d at the two Birds, 
while he himself walks forward according to the 
direction of his Piece, till his Eye is now at N, 
and the Birds in the same right Line with his 
Eye, at K and G. The Path of his Eye, described 
by this double motion, (or compounded of a pro- 
gressive and angular motion,) will be some 
Curve-line LMN, in the same Plain as the rest 
of the figure, which will have this property, that 
the proportion of the distances of his Eye from 
each Bird, will be the same every where as that 
of their respective Velocities. That is, when his 
Eye was at L, and the Birds at E and H, their 
Velocities were then as EL and HL, by the 
Construction. And when his Eye was at M, and 
the Birds at F and J, their Velocities were in 
the same proportion as the Lines FM and JM, 
by the nature of the Curve LMN. And when his 
Eye is at N, and the Birds at G and K, their 
Velocities are in the proportion of GN to KN, 
by the nature of the same Curve. And so uni- 
versally, of all other situations. So that the 
Ratio of those two Lines will always be the 
sensible measure of the ratio of those two sen- 
sible Velocities. Now if these Velocities, or the 
swiftnesses of the flight of the two Birds in this 
instance, are call’d Fluxions; then the Lines 
described by the Birds in the same time, may be 
call’d their contemporaneous Fluents; and all 
instances whatever of Fluents and Fluxions, may 
be reduced to this Example, and may be illus- 
trated by it. 

“And thus I would endeavour to give some 
notion of.Fluents and Fluxions, to Persons not 
much conversant in the Mathematicks; but such 
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as had acquired some skill in these Sciences, I 
would thus proceed farther to instruct, and to 
apply what has been now deliver’d. The con- 
temporaneous Fluents being EF=y, and HI=x, 
and their rate of flowing or increasing being 
suppos’d to be given or known; their relation 
may always be express’d by an Equation, which 
will be compos’d of the variable quantities x 
and y, together with any known quantities. 
And that Equation will have this property, be- 
cause of those variable quantities, that as FG 
and IK, EG and HK, and infinite others, are 
also contemporaneous Fluents; it will indiffer- 
ently exhibit the relation of those Lines also, 
as well as of EF and HI; or they may be sub- 
stituted in the Equation, instead of x and y. 
And hence we may derive a Method for deter- 
mining the Velocities themselves, or for finding 
Lines proportional to them. For making FG= Y, 
and IK=X;; in the given Equation I may sub- 
stitute y+ Y instead of y, and x+X instead of x, 
by which I shall obtain an Equation, which in 
all circumstances will exhibit the relation of those 
Quantities or Increments. Now it may be plainly 
perceived, that if the Line MIF is conceived con- 
tinually to approach nearer and nearer to the 
Line NKG, (as just now, in the instance of the 
Fowler,) till it finally coincides with it; the 
Lines FG=Y, and IK=X, will continually de- 
crease, and by decreasing will approach nearer 
and nearer to the Ratio of the Velocities at G 
and K, and will finally vanish at the same time, 


and in the proportion of those Velocities, that is, 
in the Ratio of GN to KN. Consequently in the 
Equation now form’d, if we suppose Y and X to 
decrease continually, and at last to vanish, that 
we may obtain their ultimate Ratio; we shall 
thereby obtain the Ratio of GN to KN. But 
when Y and X vanish, or when the point F 
coincides with G, and I with H, then it will be 
EG=y, and HK=x; so that we shall have 
y:%::GM:KN. And hence we shall obtain a 
Fluxional Equation, which will always exhibit 
the relation of the Fluxions or Velocities, be- 
longing to the given Algebraical or Fluential 
Equation. 

“Thus, for Example, if EF=y, and HI =x, and 
the indefinite Lines y and x are supposed to in- 
crease at such a rate, as that their relation may 
always be express’d by this Equation x*—ax? 
+axy—y*'=0; then making FG= Y, and IK=X, 
by substituting y+ Y for y, and x+X for x, and 
reducing the Equation that will arise, (see before, 
pag. 255).) we shall have 3x°X+3xX?+X° 
—2axX —aX*+axY+aXy+axX Y—3y°Y —3yY? 
— Y*=0, which may be thus express’d in an 
Analogy, Y:X ::3x?—2ax+ay+3xX+X?—aXx: 
3y*—ax—aX+3yY+ Y*. This Analogy, when Y 
and X are vanishing quantities, or their ultimate 
Ratio, will become Y:X :: 3x?—2ax+-ay:3y?—ax. 
And because it is then Y:X::GN:KN::4:4, it 
will be y:%:: 3x?—2ax+ay:3y?—ax. Which gives 
the proportion of the Fluxions. And the like in 
all other cases. Q. E. I.” 


An Electronic Differential Analyzer 


KeEItH R. SYMON AND ROBERT P. PoPLAwskKy* 
Wayne University, Detroit, Michigan 
(Received April 10, 1952) 


An electronic differential analyzer designed primarily for classroom teaching is described. 
The value of building such a device as a project for advanced students is emphasized. The 
method of setting up differential equations on the differential analyzer is explained, and a few 


sample solutions are illustrated. 


HE present article describes briefly a small 
electronic differential analyzer built at 
Wayne University. The project had two prin- 
cipal objectives: first, to provide an instructive 
experience in the design and construction of 


* Now with University of Michigan Research Institute. 


electronic apparatus; and second, to produce at 
a minimum cost in time and money a machine 
which would be a useful teaching aid in courses 
in intermediate and advanced physics, and which 
might also be useful in solving certain kinds of 
research problems. 











Fic. 1. Com- 
puting amplifier 
connected to 
multiply by a 
constant. 





As to the first objective, the project has been 
eminently successful. There is no substitute for 
first-hand experience in working out a design, 
building a machine, and then spending the weeks 
or months required to get the flaws out of it 
until it runs as planned. As to the second ob- 
jective, we can claim moderate success. The 
machine was designed, built, and put in working 
order by one of us (R.P.P.) as a master’s thesis 
project, representing half-time work for one 
semester, under the guidance and supervision of 
one or two staff members. At present, it is capable 
of solving linear differential equations with con- 
stant coefficients or sets of such equations up to 
about the seventh order, with a design accuracy 
of about 1 percent for each mathematical opera- 
tion performed. Its performance is continuously 
being improved, and is already entirely satis- 
factory for classroom demonstration. Additional 
components which will permit the solution of 
nonlinear differential equations and equations 
with variable coefficients are under development 
as a second master’s thesis project. The total 
cost of materials, including parts for additional 
components to be built, has been $244.00. 

It is the purpose of this article to call to the 
attention of physics teachers the desirability of 
initiating projects of this type. We hold no brief 
for the superiority of the particular design at 
which we have arrived. On all questions of de- 
sign, we made those decisions which seemed to 
us to offer the greatest likelihood of leading to 
something workable in the shortest time with the 
least expense, rather than those decisions which 
might have lead ultimately to the highest degree 
of perfection. 


PRINCIPLE OF THE COMPUTING AMPLIFIER 
The principle of operation of the computing 
amplifier has been explained by many authors.! 


1E.g., Ragazzini, Randall, and Russell, Proc. Inst. 
Radio Engrs. 35, 444 (1947). 
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P. POPLAWSKY 
In Fig. 1, the computing amplifier has a large 
negative gain —A, so that 


Ey= —AE’, (1) 


where Ey and £’ are the voltages across the 
indicated points in the circuit. Assuming that 
the amplifier has an infinite input impedance, 
the currents through resistors R; and Rp» are 
equal : 

(Ei—E’)/Ri=(E’ —Eo)/Ro. (2) 


Substituting E’ from Eq. (1) in Eq. (2), and 
solving for Eo, we obtain 


Eo= —(Ro/RYEwWL1+(1/A)(1+Ro/Ry Tt". (3) 


If A is greater than 1000, and if Ro/R, is less 
than 10, then within an accuracy of about 1 
percent, 


Eo= —(Ro/R:) Ei. (4) 


Thus the computing amplifier may be used as 
in Fig. 1 to multiply an input voltage EZ; by a 
constant negative factor. 

Figure 2 shows a computing amplifier used as 
an adder. Applying a similar analysis to that 
given in the preceding paragraph, we obtain 


Eo= — {(Ro/ Ri) E1+(Ro/R2)E2+(Ro/ Rs) Es} 
X {1+ (1/A)[1+ (Ro/R1) 
+(Ro/R2)+(Ro/Rs) ]}-, (5) 


or, to a close approximation, 


Ry Ro Ry 
Eo= ee a (6) 


1 R,. Rs 


Similarly, for Fig. 3, to a close approximation if 
A is large, we have 





1 t 
aa f Exdt, (7) 
RiCoro 


Fic. 2. Com- 
puting amplifier 
connected as ad- 
der. 
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ELECTRONIC DIFFERENTIAL ANALYZER 


where it is assumed that at t=0 the condenser 
Co is uncharged. It is clear that one can combine 
the operations of addition and integration in one 
computing amplifier circuit. In principle, a com- 
puting amplifier can be made to differentiate by 
using a feed-in condenser and feedback resistor, 
but differentiators are to be avoided, as they 
tend to amplify any hum or stray noise present. 

With computing amplifiers performing the 
basic operations of addition, multiplication by a 
constant, sign reversal, and integration, any 
homogeneous linear ordinary differential equa- 
tion with constant coefficients can be solved. 
Examples will be given in a subsequent section. 


DESIGN OF THE COMPUTING AMPLIFIER 


Several important and often conflicting con- 
siderations enter into the design of the computing 
amplifier. Perhaps the first decision to be made 
is the time required to complete a solution. A 
distinction is usually made between high speed 
computers which complete a solution in a few 
milliseconds, and low speed computers which 
complete a solution in a few seconds to several 
minutes.? With the former, the solution is re- 
peated many times per second and presented on 
the face of a cathode-ray tube, while with the 
latter, the solution is usually drawn by a re- 
cording voltmeter. Partly from perversity, and 
because we could see a number of possible ad- 
vantages, we chose a solution time of one-tenth 
second, midway between the high and low com- 
puting speeds ordinarily used. Since our com- 
puter is to be used for classroom demonstration, 
it seemed desirable to use a cathode-ray tube 
presentation. A solution requiring one-tenth 
second repeated five times per second on a 
cathode-ray tube with a long persistence screen 
gives a satisfactory presentation. By using a rela- 
tively long solution time, we hoped to avoid 
some of the high frequency response difficulties 
which arise with high speed computers. Also, at 
this speed, it is possible to use relays to control 
the computing cycle. 

The greater the gain of the computing ampli- 
fier, the smaller will be the error in each com- 
puting operation, as shown in the preceding 
section. However, increasing the gain makes 
other design problems more severe. In our de- 


2J. L. Ryerson, Am. J. Phys. 19, 90 (1951). 


Fic. 3. Com- : 
puting amplifier £, 
connected as in- ; 
tegrator. 


sign, we decided on a gain of 1000, corresponding 
to an error of about 1 percent in each computing 
operation with allowable multiplication factors 
no greater than 10 in any one operation. The 
actual amplifier gain as designed turned out to 
be about 1200. 

The frequency response of a computing ampli- 
fier must be very carefully controlled. At low 
frequencies, the gain must not fall appreciably 
at least down to the solution repetition fre- 
quency, or 1/7, where T is the solution time. 
A dc amplifier would meet this requirement, but 
dc amplifiers are subject to zero drift difficulties. 
We decided to use resistance-coupled amplifiers 
with an RC time constant of 2.5 seconds. The 
high frequency response is determined by the 
maximum time rate of change of the functions 
to be handled by the computer. We selected 1 
percent of the solution time, or 1 millisecond, as 
the minimum time within which a significant 
change in any function would be allowed. A 
frequency response up to 1000 cycles per second 
would seem to allow the computing amplifiers 
to handle functions varying at the maximum 
allowed rate. However, detailed analysis shows 
that when the computer is generating oscillating 
solutions of a differential equation, any falling 
off in high frequency response of integrators, 
adders, and constant multipliers results in a 
progressive increase in amplitude of the oscillat- 
ing voltage.* For example, in order to generate 
a sine wave of 1000 cycles per second with the 
circuit of Fig. 9, with an increase in amplitude of 
no more than 1 percent in 0.1 second, the upper 
frequency limit of the amplifiers at the half- 
power point must be at least 150 kilocycles per 
second. We based our amplifier design on a half- 
power point of about 20 kilocycles per second, 
which should allow generation of a sine wave at 


3A. B. MacNee, Proc. Inst. Radio Engrs. 37, 1315 
(1949). 
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Fic. 4. Computing amplifier circuit. 


300 cycles per second with a 1 percent increase in 
amplitude in 0.1 second. 

Since a computing amplifier is used with a 
feedback resistor or condenser connected directly 
from output to input, its phase shift characteris- 
tics are very important. Within the band of 
frequencies required for computation the phase 
shift is nearly 180°, or an odd multiple thereof, 
and the feedback is negative. The precise condi- 
tions for stability of a feedback amplifier are 
rather involved, but roughly we can say that if 
the phase shift through the amplifier and feed- 
back path becomes 0° or 360° at any frequency 
at which the amplifier gain is greater than or 
equal to one, the amplifier will oscillate. This 
implies that the phase shift must be controlled 
not only within the band of frequencies required 
for computation, but throughout the entire 
band where the gain is more than one. It can be 
shown‘ that if the phase shift is not to differ by 
more than 180° from the phase shift at mid- 
frequencies, the gain of the amplifier must not 
fall off at high frequencies more rapidly than 12 
db per octave, (that is, when the frequency is 
doubled, the voltage gain must not decrease by 
more than a factor of four). This means that if 
the gain is 1000 at 20 kilocycles, the gain can 
not drop below one up to at least 600 kilocycles. 
Since a resistance-coupled amplifier with one 
RC network per stage will have an additional 
phase shift at high frequencies of 90° per stage, 
a feedback amplifier of more than two stages is 
inherently unstable unless special steps are taken 
to keep the phase shift small. After a few experi- 


4F. E. Terman, Radio Engineering (McGraw-Hill Book 
Company, Inc., New York, 1947), third edition, p. 316. 


ments with three-stage amplifiers, we decided to 
try to eliminate the phase shift problem by using 
only one stage of amplification. _ 

Our present amplifier circuit is shown in Fig. 4. 
In it Ji, J2, and J; are input phone jacks leading 
into the amplifier through double banana jacks 
into which can be plugged the computing re- 
sistors Ri, Re, and R;, (Fig. 2). Also, Js, Js, Je, 
and J; are output phone jacks. A double banana 
jack is connected between input and output for 
insertion of the feedback resistor or condenser 
Zo. The amplifying stage is.a single 6SH7 tube 
which gives a voltage gain of about 1400 in the 
circuit shown. The two cathode follower stages 
provide a high impedance load for the 6SH7, a 
low impedance output, and a means for auto- 
matically zeroing the amplifier before each com- 
puting cycle. The over-all gain of the amplifier 
is about 1200, and the amplifier will deliver over 
100 volts positive or negative output to the load of 
25,000 ohms. Resistance values of 100,000 ohms 
to 10 megohms are used in computing. In opera- 
tion, the relay contacts K, and K, are closed for 
0.1 second before each computation begins. Con- 
tact K, grounds the amplifier input, and contact 
Kz connects the grid of the last cathode follower 
to a tap on Ry which is adjusted so that the 
amplifier output voltage is then exactly zero. 
Condenser C3; is charged up during this 0.1 
second to the proper voltage so that when K, 
and Ky, are opened, the amplifier output will be 
zero whenever the input is zero, during the 0.1 
second required to complete a solution. Resistor 
Rg is adjusted so that the voltage at the cathode 
of the first cathode follower is 150 volts. Re- 
sistor Rip is adjusted to compensate for the 
leakage through C; so that with the amplifier 
input grounded and with Kz open, the output 
remains at zero without drifting for a period of 
several seconds. Unfortunately, the amplifiers 
showed a slight tendency to oscillate when feed- 
back was applied through Zo, and this has been 
corrected temporarily by adding a 500 ywyf con- 
denser between the plate of the 6SH7 tube and 
ground. This has no apparent effect for most 
applications we have made of the computer, but 
it will of course reduce its accuracy in handling 
rapidly oscillating functions. It is hoped, when 
time permits, that a less drastic expedient can 
be worked out to eliminate the oscillations, which, 
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incidentally, were not present in the original 
breadboard model—a kind of experience shared 
by everyone who has tried to design electronic 
circuits. 

Our differential analyzer makes use of ten 
computing amplifiers like that shown in Fig. 4, 
built on a 13 in. X17 in. chassis behind a stand- 
ard 82 in. high relay rack panel. The panel 
arrangement is shown in Fig. 5. 


CONTROL CIRCUIT 


Figure 6 is a diagram of the control circuit 
which regulates the computing cycle. In auto- 
matic operation, a multivibrator alternately 
energizes and relaxes relay M. When relay M is 
energized, a 45-volt dc cycling signal is sent 
through contact M, to each amplifier, energizing 
each relay K (Fig. 4). Relay K returns the 
amplifier input and output to zero. When relay 
M relaxes, the cycling signal is removed, relay K 
in each amplifier opens, and the amplifiers are 
ready to compute. The solution time can be 
adjusted from 0.1 second to about 1 second by 
adjusting R;. The off time during which the 
amplifiers are returned to zero is about 0.1 
second. 

The cycling signal also controls relay N, which 
supplies a regulated voltage of 100 volts to phone 
jacks J;-Jio during the computing part of the 
cycle. This 100 volts is the standard dc reference 
voltage used as a basis for all computations. All 
other dc voltages required are generated by 
computing amplifiers (Fig. 1) as constant sub- 
multiples of this reference voltage. Relay N 
operates in such a way that jacks Ji-Ji are 
connected to +100 volts about 0.5 millisecond 
after contacts K, and Kez are opened on all 
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Fic. 5. Computing amplifier panel. 
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Fic. 6. Control circuit. 


amplifiers. Thus, for about 0.5 millisecond all 
amplifiers remain at zero output after K, and 
Ke are open, and then the 100-volt reference 
voltage is applied simultaneously at appropriate 
points in the computing circuit and computation 
begins. This arrangement insures that all com- 
putations shall begin simultaneously, independ- 
ently of slight differences in time of operation 
of relays K in the various amplifiers. 

Switch S; permits either manual or automatic 
control of the computing cycle. For manual con- 
trol, S: is placed in the ‘‘Stop’’ position; this 
closes relay M and allows the amplifiers to come 
to zero. Switch S; is then placed in the “‘Manual”’ 
position; this opens relay M and allows com- 
putation to proceed. 

Any desired variable in a problem can be 
plotted horizontally or vertically on the face of 
a cathode-ray tube by connecting the deflection 
amplifier to the computing amplifier output at 
which that variable is to be found. We use a 
5CP7 tube with a long persistence screen. The 
deflection amplifiers are dc amplifiers with a 
variable gain up to a maximum of about two. 
In most cases a variable is to be plotted against 
the time ¢ which can be obtained most conveni- 
ently from an integrating amplifier (Fig. 3) as 


(8) 


t 
t= 0.01 f 100dt. 
0 


By plotting ¢ from Eq. (8) horizontally, we 
achieve exact synchronization of the sweep with 
the computing cycle. By using a relay, operated 
by a trigger circuit, triggered by the cycling 
signal, it is possible to plot two variables during 
alternate cycles so that both are visible together 
on the screen. 











Ko 2 


Fic. 7. Circuit for computing U= —C+ f(aP+o)at. 
(RiC,=100/a, R2C.=}, R;C3= .. RiC3=1 b, 
R:/Rs= i and R:/Re= C). 


The precision of the computer during auto- 
matic cycling is quite adequate. For all ordinary 
problems, the solution repeats itself from one 
cycle to the next as nearly as can be seen on the 
(long persistence) oscilloscope screen. 


SETTING UP THE COMPUTER TO 
SOLVE PROBLEMS 


As a first example, let us take the differential 
equation 
dU/dt=at?+b, (9) 
U(t=0) = —C. 


In order to write Eq. (9) in a form suitable for 
the differential analyzer, we integrate to obtain 


t 
U= —c+ f (at?+5)dt. (10) 
0 
To generate the function a#*, we write 
t 
at=0.010 f 100dz?, (11) 
0 
; 
af? = 2f atdt. (12) 
0 





Fic. 8. Circuit for computing U= —17+0.02 
t 
x f'(0.0067#+-33)dt. 





S¥MON AND KR: FP. 








POPLAWSKY 





Fic. 9. Circuit for solving U+e*U=0. 


A circuit for solving Eq. (9) can then be arranged 
as in Fig. 7. 

We have so far omitted all reference to units 
and scale factors. Before Eq. (9) can be set up 
on the differential analyzer, all terms must be 
expressed in volts, and ¢ must be expressed in 
some convenient unit of time (seconds or milli- 
seconds) so that the interval of ¢ within which 
the solution is required is not more than 0.1 
second. This means that a scale factor must be 
chosen by which each quantity in Eq. (9) is to 
be multiplied in order that the resulting equation 
have the proper time scale, and in order that 
each term be expressed in volts with no term 
greater in magnitude than 100 volts during the 
solution interval. Suppose for example that we 
have the equation 


dy/dx =100x?+-200, O<x<2, 
y(x=0) = — 100. 


In order to achieve a suitable time scale, we 
introduce the variable ¢, expressed in milli- 
seconds, and write 


x =t/(50 msec), 


(13) 


0<t<100 msec. (14) 


At x=2, we have y=567, hence we introduce 
the variable U, measured in volts, 


y=6U/(1 volt). (15) 


Substituting in Eq. (13), we have (omitting 
units), “ 
50d U/dt = 0.00672+- 33, 
U(t=0) = —17. 


0<t<100. 
(16) 


A circuit with values appropriate to solving 
Eq. (16) is shown in Fig. 8. For simplicity, all 
considerations regarding units and scale factors 
will be omitted in further examples. 

A typical linear differential equation with 
constant coefficients is the harmonic oscillator 
equation 

U+wU=0, 


(17) 
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Fic. 10. Circuit for solving U+yU+o2U =f(t). (Uo =0.) 
where dots indicate time differentiation. Let us 


rewrite Eq. (17) in a form suitable for the 
differential analyzer: 


U=—, 


U= Ut f Udt, 
0 


(18) 


t 
U=Uot f Udt. 
0 


A circuit for solving this equation is shown in 
Fig. 9. 

The equation for a damped harmonic oscillator 
with forcing function f(t) is 


U+yU+U=f(t). (19) 


(a) Undamped. 


(c) Critically damped. 


Fic. 11. Circuit for solving U+(e—v(t)) U=0. 


Rewritten, this becomes 


U=—wU—yU+f(), 


t 
U= Ut f Udt, 
0 
t 
U=Uot f Udt. (20) 
0 


The corresponding circuit is shown in Fig. 10 
for the case Up=0. If the forcing function f(t) 


(b) Underdamped. 


(d) Overdamped. 


Fic. 12. Solution of equation for unforced harmonic oscillator. 
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(b) Strong coupling. 
Fic. 13. Two coupled harmonic oscillators (w:=»). 


is the solution of a linear differential equation 


with constant coefficients, it may be generated 
by an appropriate arrangement of computing 
amplifiers. If not, it must be obtained from a 
special function generating device. 


POPLAWSKY 


As a final example, consider the one-dimen- 
sional Schrédinger equation 


h? & 


—— —+ V(x) U=EU. 
2m dx? 


(21) 


Rewriting, replacing x by ¢, and introducing 
appropriate scale factors, we can put the equa- 
tion in the form 


U+(e—v(t)) U=0. (22) 


In order to solve this equation, we need a device 
to generate the function v(t), and a device to 
multiply (e—v(t)) by U. Assuming such devices, 
we can solve this equation with the circuit of 
Fig. 11. To find the eigenvalues e of Eq. (22), we 
apply appropriate initial values Up and Up» to 
satisfy the boundary conditions at t=0, and vary 
the resistor labeled 100 R/e until the boundary 
condition at the other end of the range of in- 
tegration is satisfied. The ratio of R to R/e then 
gives the eigenvalue e. Since the differential 
analyzer is limited inherently to equations in- 
volving only one independent variable (the time 
t), it cannot be used to solve partial differential 
equations unless they can be separated into a 
number of ordinary differential equations, each 
involving only one independent variable. The 
many dimensional Schrédinger equation cannot 
be solved directly on a differential analyzer, 


Fic. 14. Two coupled harmonic oscillators. 
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(c) U; alone initially excited. 


Fic. 15. Two coupled harmonic oscillators. 


although the ordinary differential equations 
which result when the variables can be separated 
are solvable. 


SAMPLE SOLUTIONS 


Figure 12 presents a series of tracings of solu- 
tions of Eq. (19) (with f(t) =0) for four different 
values of damping, showing undamped, under- 


damped, critically damped, and overdamped 
motion. 

Figure 13 shows solutions of the equations for 
two coupled harmonic oscillators: 


Ui+02U,;=kU2, 


: (23) 
U2t+w2U2=kU 4. 


Figure 13 shows U;,(t) for the case w;=we, with 
two different values of k, and with U, initially 
excited, U: at rest. Figure 14 is a study of the 
behavior of U; and U2 when U; alone is initially 
excited. Note that when w;=we, U; periodically 
exchanges all its energy with U2, whereas when 
#1 we, U; gives up only part of its energy to U2. 


Fic. 16. Wallpaper pattern. 


Figure 15 illustrates how the motion shown in 
Fig. 14a (wi=w2) can be regarded as a super- 
position of two normal modes of vibration. In 
the high frequency normal mode (curve (a)), U; 
and Uz are initially excited out of phase with 
equal amplitudes. In the low frequency normal 
mode (curve (b)), Ui and U? are initially excited 
in phase with equal amplitudes. In curve (c) of 
Fig. 15, U; alone is initially excited, and the 
motion of U2 can be seen as a beat ‘between the 
two normal modes of vibration, the beat fre- 
quency being the difference between the two 
normal mode frequencies. 

Figure 16 is an example of the use of the 
differential analyzer as a hoot-nanny. In this 
case, we have the solutions shown in Fig. 14(a), 
with U, plotted horizontally and U; vertically. 
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Shearing Stress in a Closely 
Coiled Helical Spring 


REUBEN BENUMOF, City College, New York, and Pace College, 
New York, New York 


AND 


MITCHEL BENUMOV, McKiernan-Terry Corporation, 
Harrison, New Jersey 


HE theory ef helical springs has important applica- 
tions in industry. Such springs are discussed in a 
number of books,! but unfortunately the derivation of the 
formulas is not indicated. The purpose of this article is to 
show how one of the formulas in everyday use may be 
obtained very simply. 

This brief but accurate formula for the stress in a closely 
coiled helical spring, widely used by engineers, is that given 
by A. M. Wahl.? Teachers, in setting up equipment for an 
experiment, will find it quite valuable in that it indicates 
when equipment is likely to be damaged through overload. 
To apply the formula to a given spring, a micrometer 
caliper for measuring wire and coil diameters, and a slide 
rule for calculation are the only tools needed. 

The following discussion is direct and shows how im- 
portant results may be obtained by the elementary methods 
taught in the classroom. The formulas derived in this paper 
reduce to those given by Wahl for the case of a large ratio 
of radius of helix to radius of rod (of which the spring is 
composed). The method used here is elementary in the 
sense that the usual conditions for the equilibrium of a 
section are imposed. The basic assumption made is that 
a plane section through the spring before load is applied 
remains plane after the load is applied. 

To simplify the analysis, it is considered that a single 
coil of the complete spring is isolated and two bars are 
attached as in Fig. 1(b). Obviously, the stresses in the 
isolated coil are still the same as though it were part of the 
complete spring. Since the‘coil or pitch radius is R, the 
moment acting on any section, such as P—Q, is FR, where 
F is the applied force. 

Let two sections, say P—Q and P’—Q’ in Fig. 1(b), be 
taken so as to subtend an infinitesimal central angle da. 
The length of any fiber between the sections is, then, 
(R+<x)da, where x is the abscissa of the fiber. 

Referring to Fig. 1(c), section P—Q rotates relative to 
P’—(Q’ through angle d¢ about some point (—A/, k) under 
the applied moment FR. For a straight cylindrical rod, the 
rotation would have taken place about O. Point (x, y) in 
section P—Q moves relative to point (x’, y’) in section 
P’—(Q’ through distance d¢[(x+h)?+(y—k)?]}. 

By definition, the modulus of elasticity in shear is 


G=stress/strain 
=1/(do/da)(1/R+x)[(x+h)?+(y—k)*}', (1) 
where r is the shearing stress due to torsion. Therefore, 


7 =G(do/da)[(x+h)?+(y—k)? (R+x). (2) 
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The section is in equilibrium. Consequently, 
are: _ fr FV a 
2x=0= f fr singdydx = {” awe k) 
X C(x +h)?+(y—k)? idydx 


0= —2rk[R—(R?2=--r?)*]. 
Therefore, 
k=0. 


r 6 am 
zv=0=f {> cospdydx =f" nate th) 


XC(e+h)?+(y—k)? idydx. 
Therefore, 


h=(R/2)[1—(1—r°/R)*] = (r/2)[ 1 -{1 -(=)}'] 


where m=R/r. If m becomes very large, then 
h=R/4m’. 


The load F causes a torque FR to be applied to the 
section. 


Therefore, 


(FR/G) (da/dg) = (wr?R/2)[3 — (1 —?/R?)4] 
— (2¢R*/3)(1 — (1 —r2/R?)8). 


Vv (7? —x*) 


_ yt — ay Le ThWP+ Ok) Hdyde. (8) 


F 
COMPLETE SPRING 
IN EQUILIBRIUM 
UNDER FORCES "F* 
( 


2 SINGLE COIL WITH 
BARS ATTACHED 


Fic. 1. Analysis of the action of a helical spring. (a) Plan and 
elevation views of a complete helical spring in equilibrium under forces 
F. (b) Plan and elevation views of a single coil in equilibrium. The 
short horizontal bars are imaginary, being introduced to clarify the 
physical situation. (c) Stresses and resultant strains in section P —Q of 
the single coil. 
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As before, if R/r=m becomes very large, then 
(FR/G) (da/d¢) = ar?R/2m?. (9) 
Stress 7 is a maximum at x= —r. Therefore, 
Tmax = G(dp/da)(r —h)/(R—r) 
= (F/ar*)(m)(4m—1)/(2m—2). (10) 


The above formula, which is the one given by Wahl, is 
an expression for the torsional shearing stress. The average 
direct shearing stress r’=F/zr* is usually negligible in 
comparison to tmax. However, in cases where 7’ is not 
sufficiently small, the maximum resultant shearing stress 
is, to a very good approximation, 


Tmax = (F/ar*)[(m(4m—1)/(2m—2)+1], 


where m=R/r. 


(11) 


1S. Timoshenko, Strength of Materials. Part 1. ‘‘Elementary Theory 
and Problems’ (D. Van Nostrand Company, Inc., New York, 1930), 
pp. 81-84 


.M. Wahl, Trans. Am. Soc. Mech. Engrs. 51, No. 16, 185-193, 
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An Entertaining Demonstration Using 
Window Glass as a Mirror 


ALBERT V. BAEz* 
UNESCO Technical Assistance Mission, University College, Baghdad, Iraq 


LMOST every elementary physics textbook describes 

the following illusion: The real image of a flower 
appears inside a beaker full of water. Unfortunately, the 
aberrations inherent in the large concave mirror which is 
required make the illusion visible only in a very narrow 
cone of space. 

The following demonstration has the same element of 
surprise and can be performed with an ordinary piece of 
window glass. Window glass is easily obtainable and in- 
expensive. A plane mirror is, in principle, free of aberra- 
tions, and the demonstration can be set up to be clearly 
visible to an entire seated audience. 

The principles are as old as optics. The new features of 
the following information center around the entertainment 
and instructional potential of the demonstration. 

In Fig. 1, a socket A; with a frosted bulb B, is placed 
in front of a large piece of window glass W. The bulb and 
socket are hidden from the class by a cover C, which is 
painted black inside. The virtual image of A, B, formed 
by the glass acting asa reflector is at A’, Be’. A real socket 
Az and a real bulb Bz are placed to coincide with the 
image A2’, Be’. 

The wires which go to the real socket at A2 are shorted 
there, although from a distance the student thinks that 
they are connected to Bz and controlled by a switch S. 
Actually S controls the current to bulb B; only. A foot- 
switch F, which the class cannot see, is connected as shown. 
After connecting plug P to the house current and starting 


with S open and F closed, the demonstration proceeds as 
follows: 


1. The students can view Ae, Bz directly through the 
glass. 


Fic. 1. Identical frosted bulbs Bi and Bz are placed in identical 
sockets Ai and Ag, respectively, and situated symmetrically with 
respect to the window glass W. When Bi lights, Bz appears lit. 


2. By closing S, the bulb B, lights up and its image By’ 
makes Bz appear ‘lit. The students think that S has 
turned on Bz. 

3. By unscrewing Bz and opening F it seems that Bz has 
gone out because it was unscrewed. Bz is put aside and no 
bulb appears in A». This occurs because in ordinary class- 


room illumination the image of B, is visible only when 
B, is lit. 


Fic. 2. This arrangement of the optical parts is possible when the 


lecture table has a hole in it. The electrical connections are identical 
with those of Fig. 1, 


4. Next, F is closed, whereupon B; lights up producing 
a lit image at B,’. To the students a lit bulb seems to have 
appeared in A». Time out for surprise. 

5. But just as suddenly, by opening F, it disappears. 

6. Apparently confused, the instructor puts Bz back 
into As. When he finishes screwing it in, it lights! (Or 
seems to, when F is closed, of course.) 

7. In disgust, the instructor takes a hammer and strikes 
Bz a blow which smashes it. Yet to the students who see a 
lit bulb at B,’, the bulb B seems intact and lit! (Obviously, 
B, can be a bulb with a burned out filament, and Az must 
be held firmly in its place). 


This finishes the entertainment portion of the demon- 
stration. 
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The instructional aspects concerning geometrical optics, 
transmission and reflection, virtual images, parallax, and 
the quality of the image formed by a plane mirror are 
left to the instructor. 

Figure 2 shows an arrangement suggested by Professor 
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The Behavior of a Flat Whisky 
Bottle When Heated 


N earlier notes? the mechanical effects of pressure upon 

a flat whisky bottle were reported. A conventional 
whisky bottle is equipped with a tight-fitting, cne-hole 
stopper into which is inserted a long glass tube of small 
bore. The bottle is filled with colored water; the initial 
level of the water is above the stopper in the tube. When 
the wide faces are squeezed by hand the liquid level rises; 
when the edges are squeezed the liquid level falls. Professor 
Dodd accounts for this nicely in terms of changing cross- 
sectional area. 

This note reports the very interesting scheme of things 
when the flask is heated. If a wide face is heated locally by 
directing a Bunsen flame upon it, the liquid level falls. This 
is to be expected, since the face expands, the bottle 
“rounds” up, and its volume is increased. 

Now if an edge is heated how will things go? The result 


quite astonished me. 


Jutius SUMNER MILLER 
Dillard University, 


New Orleans 22, Louisiana 


1J.S. Miller, Am. J. Phys. 18, 164 (1950). 
2?L. E. Dodd, Am. J. Phys. 18, 398 (1950). 


Kundt’s Tube Projection Demonstration 


HE usual Kundt’s tube, consisting of a cylindrical 

glass tube, was found unsatisfactory for classroom 
projection purposes because the refractions introduced by 
the curved glass surfaces rendered projection impossible. 

A rectangular box measuring about 2 in. X2 in. X30 in. 
was made of wood except that two of the long faces op- 
posite one another were made of flat glass. At one of the 
ends an ordinary earphone was attached permanently. 
Before closing the box a transparent celluloid ruler was 
cemented to one of the glass sides and some fine cork dust 
was introduced. The box was closed and screwed tight to 
make it airtight. 

The earphone was connected directly to the output of a 
variable audiofrequency generator and resonances were 
observed at several values of the exciting frequency. A 
standard-size lantern-slide projector, in which the bellows 
connecting the lantern-slide holder support and projection 
lens were easily removable, was used. By removing the 
bellows from the lantern-slide holder support it was found 
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Paul Kirkpatrick for use with a lecture demonstration 
table with a hole in it. The arrangement in Fig. 1 has been 
found quite suitable in a lecture hall with a sloping floor. 


* On leave from the University of Redlands, Redlands, California. 
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that the rectangular box could be rested on the long rods 
which supported the projection lens holder. The lens 
could then be moved out so that the celluloid scale and the 
cork dust simultaneously appeared focused in the image 
on the screen. 

The curious behavior of the cork dust particles as 
projected on the screen could be clearly seen by a large 
class. The enlarged image of the celluloid scale made a 
quantitative measure of the distance between nodes pos- 
sible. This, together with the frequency reading on the os- 
cillator, made possible a determination of the speed of 
sound in class. Since the long axis of the box was normal to 
the optic axis of the projection lens, it was possible to slide 
the box along its length and observe the pattern of the 
excited cork dust at different places within the box. 

As an aside it may be mentioned that during the same 
lecture a loud-speaker placed near a soap film whose image 
was projected on a screen produced interesting patterns at 
certain frequencies (eigenvalues) impressed on the loud- 
speaker by the audiofrequency oscillator. The fact that 
the films produce colors by interference makes them more 
attractive, although the colors are irrelevant to the sound 


vibration demonstration. 


ALBERT V. BAEz* 
UNESCO Technical Assistance Mission, 


University College, Baghdad, Iraq 


* On leave from the University of Redlands, Redlands, California. 


Remarks Concerning the Paper “Conditions for 
the Derivation of the Stress 
Deviator Tensor” 


BELIEVE that it may be appropriate to make the 

following comments (largely relating to terminology) 
concerning O. L. Anderson’s interesting paper! of the above 
title. 

The author seeks the ‘physical requirements, i.e., 
conditions, of such a decomposition”’ (of the stress tensor 
into a deviator and a spherical component). Now this 
cannot mean the physical conditions which must be 
satisfied in order that such a decomposition be possible, 
for these conditions are neither more nor less than the 
conditions which make it possible to represent stress by a 
a second-order symmetric tensor. For any physical quantity 
which is representable as a second-order symmetric tensor, 
e.g., the inertia tensor of a rigid body in dynamics, such 
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decomposition is possible. It appears, therefore, that the 
author is concerned not with the possibility of such de- 
composition but rather with its usefulness, and in par- 
ticular, its usefulness in the theory of plasticity. In this 
connection we remark that such decomposition had 
previously been used in other disciplines, such as elasticity 
theory? and hydrodynamics,? and was carried through in 
substantially the same way as described by the author, 


which he appears to attribute to Russian writers on 
plasticity. 


Since we are concerned with fundamental physical | 


meaning, it may be well to mention that the yield criteria 
listed are those for a perfectly plastic material, which some 
real materials approximate; they are false for work- 
hardening materials,‘ although the criterion can be ex- 
pressed in terms of the invariants of the stress deviator ; 
they are false for metal single crystals,5 and the criterion 
for yield is not expressible in terms of the invariants of the 
stress deviator in this case. 

In the third paragraph the author says: ‘“‘The question 
arises, are there other ways of decomposing S (the stress 
tensor) so that one of its constituents produces the desired 
yield criteria? If this decomposition is unique, what 
physical meaning is implied by this uniqueness itself? 
What is implied, in terms of elastic theory, by making this 
decomposition unique?”’ I am ata loss as to the meaning of 
the last statement. It seems to me that the decomposition 
either is, or is not, unique in the sense given by the first 
sentence. How.can we make it unique? 

We note that the results in the second section of the 
paper may be deduced® from (1) definition of stress, (2) 
definition of (infinitesimal) strain, (3) symmetry of the 
stress tensor, (4) generalized Hooke’s law, and (5) elastic 
isotropy. In particular, if item (5) is dropped, the results 
are quite different, there being in the least symmetric case 
21 independent elastic constants instead of two. 

In the third section the author assumes (in addition 
to the implicit assumption of the five items just mentioned) 
“two postulatory conditions: (1) the stress and strain 
tensors can be decomposed into two and only two funda- 
mentally different independent constituents, and (2) the 
total elastic strain energy is simply the sum of the strain 
energy of the individual constituents.” It would have been 
sufficient for the succeeding deductions to have assumed 
that the stress and strain tensors can be decomposed into 
two constituents such that the total strain energy is the 
sum of the strain energies associated with the two con- 
stituents taken separately. In particular, the author’s 
characterization of ‘fundamentally independent” appears 
to be irrelevant to the argument which he presents. 

In the final section of the paper the author states: “If a 
strain is proportional to a stress, i.e., if es; =aaiz, Hooke’s 
law is valid for all components of the stress and strain. 
It is well known that this proportionality is not observed 
in the general strain-stress relationships as shown by Eq. 
(12) which indicates that each of the «xz terms is a linear 
function of all the ox; terms. The lack of a causal relation- 
ship is more forcibly demonstrated for the case of simple 
axial extension in which there is lateral strain with no 
component of stress in the direction of that strain. The 


THE EDITOR 65 


question arises, what causes the lateral strain? The concept 
of Poisson’s ratio cannot explain away the above, for this 
number is experimentally determined and has no theo- 
retical significance.”’ The set of equations oi; = Ci;x1€%1 ex- 
pressing the stress components as linear functions of the 
strain components is known as the generalized Hooke’s 
law® and is usually regarded as expressing a causal connec- 
tion between stress and strain. The author’s implied re- 
striction on the concept of causality is an innovation of 
doubtful value, since it would deny causality in numerous 
situations where a functional relationship exists among 
three or more physical entities, while affirming causality 
where a functional relation involves only two physical 
entities. The answer to the question raised by the author 
is, in my opinion, that the lateral strain is caused by the 
longitudinal stress. There is nothing peculiar about this 
behavior if it is viewed in the light of the general aniso- 
tropic case, where a change in a single stress component 
produces an alteration of every strain component. That 
the Poisson ratio is determined by experiment is not 
unusual; all elastic constants are thus determined. I do 
not see that the Poisson ratio has any more or any less 
theoretical significance than any other elastic constant. 

In conclusion, the physical significance which the 
author attaches to the stress deviator and the spherical 
stress tensor cannot be of general validity. If the number 
of independent elastic constants exceeds two, the argument 
must necessarily fail. Perhaps in the general anisotropic 


case it may be possible to decompose the stress tensor 
21 
Sand the strain tensor E into 21 components, S= Z S; and 


t=1 
21 
E= 2 Ei, such that E;=a;S;(i=1, 2, 


, 21). However, 
t=1 


even if this could be done, and if the decomposition could 
be proved to be unique, I should still be unwilling to 
concede that the S; enjoyed any special physical signifi- 
cance not enjoyed by any other set of constituents. Such a 
resolution might be useful, but I do not see that it would 
be physically more meaningful. 


MILton O. PEACH 
University of Notre Dame, 


Notre Dame, Indiana 


- ~?~ Js Phys. 20, 236 (1952). 
H. Love, A Treatise on the Mathematical Theory of Elasticity 
(Cambridge University Press, Cambridge, 1927), fourth edition, p. 83. 
3H. Lamb, Hydrodynamics (Cambridge University Press, Cambridge, 
1932), sixth edition, p. 574. 
4A. M., Freudenthal, The Inelastic Behavior of Engimeering Materials 
and Structures (John Wiley and Sons, Inc., New York, 1950), p. 301. 
5 See reference 4, Pee 129. 
S. Sokolnikoff, Mathematical Theory of Elasticity (McGraw-Hill 
Book Company, Inc., New York, 1946), Chap. III. 


Extension of Remarks on the Derivation of the 
Stress Deviator Tensor 


HE remarks! of Mr. Peach concerning my paper,? 
“The Derivation of the Stress Deviator Tensor,” 
deserve further explanation. 

In his second paragraph, Mr. Peach felt that the paper 
had to be concerned with either the possibility or usefulness 
of the decomposition of the stress and strain tensor, 
implying that there were no other reasons for the paper. 
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Now the paper was not intended to show that the de- 
composition was possible, for it should be evident that any 
tensor quantity, symmetric or otherwise, can be decom- 
posed in an infinite number of ways. The paper was not 
intended to show that the decomposition was useful, a 
point I emphasized in the first paragraph of my paper. 
References 1, 2, 3, and 4 of my paper were listed to give 
ample evidence of usefulness of this decomposition. The 
paper was intended to show that the particular decomposi- 
tion was significant, a possibility Mr. Peach apparently 
overlooked. 

As Mr. Peach says, ‘‘this decomposition has previously 
been used in elasticity theory and hydrodynamics,” but 
the fact that the decomposition is well known is irrelevant. 
The point of the paper is (a) workers in plasticity have 
found this decomposition particularly useful, (b) this 
decomposition is significant because it requires that the 
strain energy have a simplified form, and (c) this form 
implies that the constituents of stress and strain are inde- 
pendent of each other, a fact that has not been noted 
previous to this paper. 

In Mr. Peach’s third paragraph, he states that the yield 
criteria are ‘false for work-hardened materials . . . and false 
for metal single crystals.’’ This should not be so surprising, 
since these materials are excluded by the assumptions made 
in the derivations of the yield criteria. In work-hardened 
materials, there is an internal distribution of strain which 
is unrelated to the external stress; hence, there is no well- 
defined strain-energy relation. For metal single crystals, 
there are more than two independent elastic constants, and 
the strain-energy expression has a form different from that 
used in the derivations in my paper. 

In Mr. Peach’s sixth paragraph, he makes the following 
statement, ‘‘It would have been sufficient for the succeeding 
deductions to have assumed that the stress and strain 
tensors can be decomposed into two constituents such 
that the total strain energy is the sum of the strain energies 
associated with the two constituents taken separately. In 
particular, the author’s characterization of ‘fundamentally 
independent’ appears to be irrelevant to the argument he 
presents.” Mr. Peach erroneously confuses an equivalent 
argument with an irrelevant argument. Let us follow Mr. 
Peach’s line of reasoning. 

Let A=Ai+A:2 and B=B,+B:; then let us assume that 
A-B=A,-B,+A:2-B: for all values of A and B. 

This means that A; acting on B: plus A: acting on B, 
is null. Now if Ai were not independent of Be, and likewise 
if A. were not independent of B;, we could take arbitrary 
values of Ai, Ae, Bi, Be such that the sum 


Ai-B2+A:2-Bi, 


would be different from zero, contrary to our assumption. 
I cannot agree that the characterization of “fundamentally 
independent” is irrelevant to the energy argument. It is 
rather equivalent. 

In Mr. Peach’s seventh paragraph, he objects to my 
restriction on causality. The concept of causality is a 
philosophical question and is used variously by different 
authorities. Let us avoid the work and discuss the meaning 
of Hooke’s law itself with regard to lateral strains. In G. 
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Joos’ Treatise on Theoretical Physics, we find: ‘‘the founda- 
tion of the so-called classical theory of elasticity is Hooke’s 
law, which states that the normal stress is proportional to 
the elongation per unit length (ut tensio, sic vis).”’ The 
Hookian relation does not include lateral strains. While 
the generalized Hooke’s law is a useful functional relation- 
ship, it is not a fundamental one since there is no principle 
of mechanics which provides for orthogonal responses to 
forces. The cause of a lateral strain can only indirectly be 
the longitudinal force. The direct cause must be a stress 
component in the direction of the lateral strain—other- 
wise Hooke’s law cannot explain the lateral strain. That 
such a stress component exists is shown by my paper. 
Consider an elemental domain under a simple tensile load 
along the x axis. According to the theory, we have for the 
dilatational stress system, 


a = ce = Pat = 4P, 
and for the distorting stress system, 
oz’ =3P, oy = =P, o;/= =—$F. 


These two stress systems are equivalent to a tensile load, 
since oz'+o,""=P, o,'+0,'’=0, and o,’+a,’’=0. Yet, they 
predict a net lateral strain since the elastic modulus for the 
distorting and dilatational stress systems are different. In 
this example, there is no reliance on a generalized Hooke’s 
law. Hooke’s law is used, as he gave it, to produce a lateral 
extension. It is necessary to assume, however, that there 
are two fundamentally different types of stress systems. 
To argue that the tensile load produces the lateral strain 
misses the point. Unless one relies cn a complicated argu- 
ment based on directed valence forces between atoms (of 
which we know too little) a normal force cannot be con- 
sidered as the prime cause of a lateral strain and be con- 
sistent with the principles of Newtonian mechanics. 

In the seventh paragraph, Mr. Peach states, “I do not 
see that the Poisson ratio has any more or any less theo- 
retical significance than any other elastic constant.’’ While 
the magnitude of the elastic constants are all determined 
by experiment, there is a difference in our knowledge 
concerning the derivation of these quantities. The bulk 
modulus is defined by a precise thermodynamic equation 
and can be calculated provided one knows the equation of 
state. The other elastic moduli (excluding Poisson’s ratio) 
are all related to the second derivatives of the potential of 
the chemical binding energies of the atoms themselves. 
While the evaluation of these elastic moduli is difficult 
in practice, they are nevertheless well known in principle 
and are as fundamental as our knowledge of atomic forces 
themselves. Poisson’s ratio is named after a French mathe- 
matician, who determined it using a postulated molecular 
theory of structure. Poisson’s model does not account for 
chemical binding energies in the form known to us today, 
however, and we cannot consider it as having theoretical 
significance. All we can say for the Poisson ratio today is 
that it is a ratio of-certain elastic moduli. I submit that 
the elastic moduli, excepting Poisson’s ratio, have a funda- 
mental theoretical significance. 

To dismiss the previous arguments on the premise that 
they are only valid for isotropic materials (as Mr. Peach 
would do in his last paragraph) is contrary to the spirit of 
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the development of the science of elasticity. I dare say 
over 95 percent of all elasticity papers are valid only for 
isotropic materials. Airy’s function itself is derived with the 
assumption of isotropy. Our knowledge of the mechanics 
of crystal systems is meager compared to our knowledge 
of the mechanics of isotropic systems. This does not prevent 
us from making progress, restrictive as it may be, when 
we can. 

While Mr. Peach does not agree that the decomposition 
of the tensors made in my paper is unique, he fails to offer 
another set. I should be very interested to be shown another 
set of constituents subjected to the conditions of Eqs. (11), 
(12), and (13) of my paper. 

There is yet one minor point of contention. If Mr. Peach 
will read the paper carefully, he will note that I do not 
attribute the original decomposition of stress to Russian 
authors of plasticity; I only give their method by way of a 


review of the problem. 


Orson L. ANDERSON* 
University of Utah, 
Salt Lake City, Utah 


1 Milton O. Peach, Am. J. Phys. 21, 66 (1953). 

20. L. Anderson, Am. J. Phys. 20, 236 (1952). 

* Presently a member of the Technical Staff, Bell Telephone Labora- 
tories, Murray Hill, New Jersey. 


Heat Transfer to a Boiling Liquid 


HE letter by Mr. Edwards! on the phenomena ob- 

served when a warm body is immersed in liquid air 
recalled the following measurements made several years 
ago on the same phenomena. 

When a warm body is placed in liquid air, boiling occurs 
at the warm surface. Vaporized liquid air is evolved at a 
moderate and steady rate until just before the body 
reaches the temperature of the liquid air. Then occurs very 
violent boiling that quickly subsides as the body attains 
liquid air temperature. 

This phenomenon is not a property of liquid air, but of 
boiling liquids. The same series of events occurs when a hot 
body is dropped into boiling water. The two types of heat 
transfer, the efficient ‘‘contact” transfer and the less 
efficient transfer through a gas film, are well known to heat 
engineers,” but not the mechanism suggested below for the 
former. 

A conical piece of copper was set base-upward into a 
sleeve of nickle alloy tube so that only a circular surface of 
copper was exposed. The sleeve was an alloy with low 
thermal conductivity. An electrical coil was wound on the 
cone. A constantan wire was soldered at the base of a small 
hole that reached almost to the exposed surface of the 
copper cone from inside the sleeve. This formed a thermo- 
couple junction that was connected to measure the tem- 
perature difference between the copper surface, heated by 
power metered to the coil, and surrounding liquid nitrogen 
when the apparatus was immersed therein. The sleeve kept 
liquid nitrogen off the unexposed parts of the cone. 

Stable temperature-power conditions were obtained over 
the range shown in Fig. 1 for “‘contact”’ transfer. The points 
on the two curves were not all obtained in sequence, as 
more and more power was applied; but the power succes- 
sively increased and decreased, each time retracing the 
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Fic. 1. Temperature difference between liquid nitrogen and sol'd sur- 
face vs rate of heat transfer. The two lower curves represent ‘‘contact”’ 
or thin-film transfer, the upper for a copper surface, the lower for a 
copper surface sheathed in ice. The upper curve represents heat trans- 
fer through a thick film of gas. 


same curve. For power inputs over 72 watts/cm? the 
“contact” type of transfer ceased and the temperature of 
the surface rose several hundred degrees. The lower curve 
with the reciprocal slope of 2.8 watts/cm? deg was ob- 
tained when a thin film of ice coated the copper surface. 

The “‘cooling’’ curve represents heat transfer through a 
gas film. The surface temperatures represented by points 
on this curve were not as stable as those for contact cooling, 
but always the ‘‘contact” type of transfer took over when 
the power was reduced below 26 watts/cm?. 

The word ‘‘contact”’ transfer suggests the unsatisfactory 
picture of the liquid wetting a surface whose temperature 
exceeds the boiling point of the liquid. It is here proposed 
that so-called ‘‘contact”’ transfer is really transfer through 
a film of gas whose thickness is less than the mean free 
path of molecules in the gas. Into this gas layer molecules 
evaporate from the liquid surface and strike the hot solid 
surface without intermediate collisions. There they take 
up energy and return to the liquid. 

The theory of this type of heat transfer? shows the heat 
transfer independent of the thickness of the gas film, 
roughly proportional to the pressure and the temperature 
difference between the bounding surfaces, and dependent 
upon ‘‘accommodation coefficients” for the gas in contact 
with each surface. Accommodation coefficients are one or 
less and represent the effectiveness of the surface in im- 
parting energy to the molecules. 

If the accommodation coefficient for the liquid nitrogen 
surface is assumed to be one (it may be less), then a not 
unreasonable accommodation coefficient for the copper 
surface of about a=0.35 will account for the upper curve 
and a=0.54 for the ice-coated surface of the lower curve. 
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These values assume that the whole surface is conducting 
and no attached bubbles reduce the effective area. 

The “noncontact” heat transfer is effected through a 
film of gas much thicker than the mean free path. If 
ordinary conduction theory is applied to the last stable 
point on the ‘‘cooling’’ curve and the observed data used 
to calculate the average film thickness, one finds that the 
thickness is about 10 mean free paths. No account of 
temperature jumps‘ at the surfaces entered this estimate. 
With jumps included the average thickness would be less. 

It is easily conceivable that the switch from gas film heat 
transfer to “contact,” or more appropriately, thin-film 
transfer, occurs when a reasonable fraction of the surface 
is covered by a film only a few mean free paths thick. 

J. Etmer RuHopEs, Jr. 


Georgia Institute of Technology, 
Allanta, Georgia 


1M. H. Edwards, Am. J. Phys. 20, 313 (1952). 


2M. Jakob, Heat Transfer (John Wiley and Sons, Inc., New York, 
(1949), Vol. I, Chapter 29. 


3 E. H. Kennard, Kinetic Theory of Gases (McGraw-Hill Book Com- 
pany, Inc., New York, 1938), p. 316. 
4 See reference 3, p. 312. 


Frequencies Resulting from Distortion 


N a recent Letter to the Editor' Howard N. Maxwell 
asks for a physical explanation of the additional fre- 
quencies found by a Fourier analysis of a distorted mono- 
chromatic wave train. It is evident that such frequencies 
are not due to resonances in the distorting device itself, 
whether that is an electrical rectifier or a nonlinear me- 
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chanical device. Where then do the additional frequencies 
arise? 

This question reminds me of the hot controversy con- 
cerning the real existence of sidebands in amplitude- 
modulated radio transmission, fought out in the popular 
radio journals about twenty years ago. If a carrier wave 
of frequency F is modulated at a much lower frequency f, 
a suitable wavemeter will show responses at the sideband 
frequencies Ff; yet there is no circuit in the transmitter 
tuned to either of these sideband frequencies. The con- 
clusion then and now accepted is that we have two mathe- 
matical ways of representing the phenomenon, both 
equally valid; namely, an equation showing it as a modu- 
lated wave, and another showing it as the carrier plus two 
sidebands. 

When Helmholtz analyzed the wave forms set up by 
musical instruments, he obtained harmonic responses in 
his cavity resonators; and he indicated the modes of 
vibration of strings and air-columns that would account 
for the presence of the harmonics in the wave form. 
Perhaps this is why we instinctively look for similar 
mechanisms when we observe a distorted wave of any 
kind. However, what is appropriate in acoustics is not 
necessarily so in other branches of physics. If Dr. Maxwell 
is asking whether we can identify an original source of each 
harmonic in a distorted wave, no matter what may be the 
cause of the distortion, I suggest that a wish to do so may 
be a survival from the too-mechanistic point of view of 


the nineteenth century. 


Cuas. WILLIAMSON 
Carnegie Institute of Technology, 
Pittsburgh 13, Pennsylvania 


1 Howard N. Maxwell, Am. J. Phys. 20, 310 (1952). 
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Book Reviews 


Industrial High Vacuum. J. R. Davy. Pp. 243+viii, 
Figs. 82, 14.522 cm. Sir Isaac Pitman & Sons, Ltd., 
London. Price 25s. 


If high vacuum is a technology which sometimes rises to 
the dignity of a science, it is also an art and a cult. Devotees 
chat and reminisce, and compare notes and recipes in the 
language of their cult, safe in the knowledge that they will 
be understood by their brethren. They even write books— 
indeed many of this kind have been written in the last few 
years. These volumes must be distinguished sharply from 
the masterworks of the giants—Dushman, Dunnoyer, 
Espé and Knoll, Strong—yet they have individual useful- 
ness as each devotee describes the rituals in which he has 
participated. 

The volume under consideration deals expertly, and in 
more than half its pages, with the production and resulting 
properties of thin films of metals, metalloids, and salts; 
the manufacture of mirrors, antireflection coatings, and 


selective transmission filters. Just how to set about these 
tasks is detailed from the ground up—the selection and 
fitting of the room, the components of the apparatus, 
choice of operating chamber with its demountable leads, 
control knobs, etc. The source-furnace and raw materials 
for evaporation are given in variety by one who is thor- 
oughly familiar with them. Particularly good are the sec- 
tioned sketches of gaskets, rotary seals, and other closure 
devices; and the directions for discovering and curing leaks. 

The basic or peripheral information, which would be 
required by those’ not yet of the cult, is scanty in the ex- 
treme. There is one picture of a condensation pump, of a 
type little known in the United States. The fractionating 
pump, the choice and properties of pump oils, and the art 
of making and maintaining a high vacuum are virtually 
absent. The chapter on measuring vacuum is equally 
sketchy—a cultist would, of course, know where to look 
for information! In this vein, the section on kinetic theory, 
which occupies the first four voluminous chapters of Dush- 
man, is relegated to the last eight pages of the present book. 
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A relatively comprehensive chapter on Molecular Dis- 
tillation is included, presumably because of the early im- 
portant contributions made in the field by British scientists 
and manufacturers. The story has been compiled from 
scientific articles and the patent literature but is, on the 
whole, well told. 

If, then, we rate this book on what it treats, it is a good 
book. If we rate it on the contents implied by its title, it 
is less admirable. But for anyone requiring a treatise on 
“The Art and Practice of Surface Film Deposition in 
Vacuo,”’ the volume will be indispensable. Also, it pro- 
vides the American operator with a first-hand description 
of British methods. 

KENNETH HICKMAN 
Rochester, New York 


The Astronomical Universe. WasLEy S. KROGDARL. 
Pp. 599+-xvii, Figs. 275, 16X24 cm. The Macmillan 
Company, New York, 1952. Price $6.25. 


Professor Krogdahl offers a distinctively new and differ- 
ent approach in the field of astronomy to the textbook 
problems created by the changed emphases in preparatory 
education and the newer evaluations of the objectives of 
college science courses for nonmajors. 

This is not a watered-down conventional text with 
frustrating inadequacy of explanations, nor a purely de- 
scriptive book, but a serious and skillful effort to meet both 
the proper objectives of a science course and the abilities 
of the typical liberal arts student as well as to provide a 
solid, stimulating account of astronomy. The preface indi- 
cates the author’s belief in its suitability for astronomy 
majors also. For the latter purpose it suffers serious inade- 
quacies. On its effectiveness for the primary purpose, there 
will probably be as many divergent opinions as there are 
teachers of astronomy; only experience with classes will 
give an answer. The subsequent remarks are merely one 
teacher’s preliminary reactions. Certainly the text merits 
serious examination by all teachers of astronomy for general 
students. 

The book’s distinctions are a new approach to the order 
of topics and to the method of discussion and explanation, 
and unusual skill in exposition and vitality of expression. 
Specific, articulate emphasis is placed on the scientific 
method, which is discussed in the introduction in a general 
way. Section I opens unconventionally with a complete 
account of Kepler’s and Newton’s laws and gravitation, in 
which the operation of the scientific method is brought out 
descriptively but with little analysis of the specific reason- 
ing. Other unusual orderings of topics appear throughout 
the book for the same reason. 

The historical order is ignored. Names and dates are 
almost completely omitted as superfluous burdens to the 
student. The order of discoveries is not brought out. 
While these omissions may make possible more internally 
consistent demonstrations of the scientific method, they 
involve a loss of perspective of the science and leave one 
with a curious sense of lack of moorings. A more historical 
order, developing topics from general, historically presented 
observations, would deepen the student’s appreciation of 
the nature of science and need not interfere with, but 


would often enhance, his comprehension of the scientific 
method. 

The subject is developed with rigor and thoroughness in 
spite of avoidance of algebraic formulations. This avoid- 
ance is skillfully and successfully accomplished without 
loss of quantitative reasoning in the treatment of the solar 
system; gravitational problems are presented clearly in 
successive steps without the use of complete formulas. In 
physical problems, there is less success in this; radiation 
laws are used in algebraic form where a similar approach 
would be more real to average students. 

The writing is exceptionally clear, efficient, lively, and 
characterized by an admirably light touch. There is a 
tendency to lose some clarity, occasionally, by use of what 
most students would call (deplorably) “big words.’’ In 
discussions of many topics in modern astrophysics ambition 
for completeness and thoroughness has outstripped the 
typical student’s digestive ability for brief, almost purely 
verbal explanation of complicated concepts. His ability and 
need to comprehend the topics in such full detail is highly 
questionable. 

The book is very up to date. Brief accounts, rather 
cursory, are given of such topics as radio astronomy, shell 
stars, and recent thinking in cosmogony. The speed of 
advance in the science has caused a few inconsistencies; 
the account of Jupiter’s satellites includes the twelfth, but 
“Jupiter’s outer three” are referred to in the discussion of 
Saturn’s satellites, for example. 

Illustrations are well chosen and, in the main, well re- 
produced. Line drawings are meticulous and original, 
though not always distinguished by ready clarity, and often 
lacking where needed. Extensive and comprehensive ques- 
tion and problem lists accompany the chapters, many of 
them quite original and thought provoking, and some of 
them real teasers. Star charts (transverse Mercator pro- 
jections), with white stars on blue background, are well 
drawn with minimum distortion, but suffer from small 
scale and inclusion of too many faint stars, looking too 
bright. 

Puitip S. RiGGs 
Drake University 


Ultrasonic Physics. E. G. RicHARDsON. Pp. 285+x, 
Figs. 100, 15.5X23.5 cm. Elsevier Publishing Com- 
pany, Houston, Texas, 1952. Price $5.00. 


The field of ultrasonics has expanded during the past 
quarter century to a point where it does not seem possible 
to cover the entire subject comprehensively in any single 
volume. E. G. Richardson has wisely chosen to restrict his 
new book to those phases concerning the propagation of 
ultrasonic waves in relation to molecular physics. In several 
instances, associated material involving ultrasonic waves 
has also been included. 

Prime emphasis is placed on the construction and use of 
acoustical interferometers. After introductory chapters on 
the production and detection of ultrasonic waves, the 
author describes experimental techniques for the deter- 
mination of velocity and absorption in gases, liquids, and 
solids. Only sufficient theory is given to show generally the 
relation of these parameters to molecular structure. One 
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chapter is devoted to summarizing very briefly the theories 
concerning anomalous velocity and absorption in various 
media. The final chapter describes the effects associated 
with the propagation of ultrasonic waves through dispersed 
systems. A bibliography follows each of the eight chapters. 

In this book the author is at his best in the discussion of 
research with the acoustical interferometer, particularly 
those aspects with which he has had personal experience. 
In several cases, however, he seems to have sacrificed 
clarity for brevity to an undesirable extent. To cite a few 
examples, the discussion of cavitation (page 56) not only 
is inadequate but is actually misleading in terms of gener- 
ally accepted explanations of cavitation. In the treatment 
of the effect of electrolytes on ultrasonic absorption, data 
are presented without any statement as to the frequency. 
The author uses in excess of three pages to discuss the 
Debye effect involving ionic vibration potentials and yet 
fails to indicate that the effect has been experimentally 
verified for simple ionic solutions (as of 1949). It is regret- 
table that there is not a discussion of the spherical resonator 
technique used by Tamm in Germany and Leonard and 
his co-workers in this country for the determination of 
ultrasonic absorption. This method seems to be the only 
one capable of yielding results for liquids of moderate ab- 
sorption at frequencies from 30 to 500 kc/sec at the present 
time. 

In general, the organization of this book is adequate. 
With the possible exception of the first chapter, the ma- 
terial is presented in such a fashion as to constitute rela- 
tively ‘‘easy reading.” 

This book should prove interesting to the person who 
has relatively little experience with ultrasonic techniques 
and would like to learn of their potentialities as a tool for 
the study of molecular structure. It should be within the 
grasp of any advanced undergraduate student in physics 
and physical chemistry. 


EARNEST YEAGER 
Western Reserve University 


Elements of Wave Mechanics. N. F. Morr. Pp. 156+ ix, 
Figs. 30, 14.522 cm. Cambridge University Press, 
New York. Price $3.75. 


To those acquainted with the author’s Outline of Wave 
Mechanics, published in 1930 and now out of print, Pro- 
fessor Mott’s Elements of Wave Mechanics will come as a 
welcome addition to the literature. The book, like its 
predecessor, is intended for students who would presumably 
correspond, in this country, to undergraduates in their 
senior year. Assuming that the training of experimental 
physicists at this level should differ from that of students 
of theoretical physics, the book is aimed at the former 
group. It is, of course, to be regarded as an introduction to 
more advanced treatises, as the author is careful to point 
out. With this in mind, it is understandable that the book 
is very brief and that many subjects must be treated in an 
almost cursory fashion. 

As in the Outline the approach is confined to the wave 
treatment, and no mention can be made of the more power- 


ful methods and more general point of view which become 
available in a more advanced treatment. Continuing the 
comparison from the point of view of subject matter, the 
present book not only brings the discussion up to date 
but it covers a somewhat wider range of topics than does 
its predecessor. Beginning with a mathematical introduc- 
tion, which includes a brief discussion of differential equa- 
tions and Fourier’s theorem, the author proceeds to the 
quantum-mechanical description via the Schroedinger 
equation. After a very much simplified discussion of the 
wave properties of matter and of the uncertainty principle, 
a more detailed account of stationary state and dynamical 
processes follows. As is to be expected in a book with the 
avowed purpose of serving as a brief introduction to the 
subject, the choice of topics must be a conventional one. 
For example, under the heading Stationary States, one 
will find discussions of the hydrogen atom, harmonic oscil- 
lator, perturbation and variational methods, effects of 
electric and magnetic fields, and electron spin. Of course, 
the treatment of these various problems in some cases may 
serve only to whet the appetite, and this may very well be 
all to the good. The many-body problems also receive 
attention, and here a simplified version of the problems 
of molecular structure, chemical binding, and the solid 
state is presented in some 37 pages. The chapter on quan- 
tum dynamics has been extended and amplified consid- 
erably as compared to the Outline. The book closes with 
brief, largely qualitative comments concerning relativistic 
wave mechanics and nuclear physics (alpha- and _ beta- 
decay, in particular). Throughout the text a number of 
appropriate problems have been interspersed. 

There are a number of pedagogical devices which appear 
in the book and which were particularly pleasing to this 
reviewer. For instance, the early introduction of the WKB 
solution of the wave equation for the purpose of providing 
a (qualitative) understanding of how a wave function 
behaves is to be applauded. On the debit side the excessive 
use of reasoning by analogy is perhaps subject to adverse 
criticism. Thus, the wave equation is introduced by the 
conjunction of the Davisson-Germer wavelength-mo- 
mentum relation and the scalar wave equation describing 
sound waves (or superficially, electromagnetic waves in 
vacuum). This, it seems to this reviewer, confuses the 
issue since DeBroglie waves are like neither of these 
classical waves physically. This is exemplified by the 
well-known fact that the procedure adopted does not lead 
to the correct time-dependence in a natural way. It seems 
that, without violating the requirements of brevity and 
simplicity, the less objectionable procedure of utilizing 
Fourier integrals, the DeBroglie-Bohr relationships, and 
the classical Hamiltonian could have been followed. All 
the machinery for this is, in fact, introduced before the 
Schroedinger wave equation (time-independent) is written 
down. This may be a matter of taste, but it would appear 
to place the postulational basis of wave mechanics in 
evidence in a clear way. 


M. E. ROosE 
Oak Ridge National Laboratory 
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Mechanics and Properties of Matter. REGINALD J. 
STEPHENSON. Pp. 371+ x, Figs. 175. 1523.5 cm. 
John Wiley and Sons, Inc., New York, Chapman and 
Hall, Ltd., London, 1952. Price $6.00. 


This book is a new textbook for the intermediate me- 
chanics course. The student is expected to have studied 
general physics and calculus for at least one year each; the 
treatment is designed so that the prior knowledge of 
physics which is assumed is quite modest, and students 
who have had a thorough two-year course in general 
physics will find a certain amount of familiar material. The 
subjects covered are the mechanics of the particle and the 
rigid body including oscillatory motion and gyroscopic 
motion and, in addition, there are short treatments of elas- 
ticity, hydrostatics and hydrodynamics, and wave motion. 

The author has tried to emphasize the physical content 
of the theory rather than the mathematical formalism. 
With that purpose the vector algebra has been introduced 
as needed through the book so that the scalar product 
appears first in the chapter on work and energy, and the 
vector product is not introduced until the author is well 
along in the chapter on rigid body motions. For some of the 
more complicated subjects the author has given a first 
approximate treatment by way of introduction before 
taking up a more thorough-going discussion. Mathematical 
points with which the student may be unfamiliar are 
usually explained carefully as in the discussion of conic 
sections in the chapter on gravitation. There are many 
illustrative problems whose solutions are given in detail, 
and each chapter concludes with a large collection of prob- 
lems of various kinds. 

There are several deficiencies in the book which seem 
important to the reviewer. After a careful discussion of 
the definition of mass, the author introduces the concept 
of force without definition and does not discuss the logical 
status of Newton’s second law in his development. The 
treatment of rotational motion is lengthy, including a 
discussion of the motion of the top, but the author does 
not always emphasize the basic physics of the problem. 
Thus, the vector nature of the angular velocity is not 
adequately discussed: it is not made clear that although 
angular velocities are vectors, angular displacements are 
not. The relationship between angular momentum and 
angular velocity in not stressed, so that the student may 
not grasp the basic point that these two vectors are not, 
in general, parallel to one another. The inertia tensor is 
mentioned only briefly and on page 216 tensor is mis- 
defined and wrongly exemplified. Finally, there is no dis- 
cussion of the important subject of coupled oscillations and 
normal coordinates. 

One last point on which the reviewer would like to 
comment is the author’s general approach to the teaching 
of intermediate mechanics. The author’s method is to 
continue the elementary approach of the general physics 
course in dealing with the subject matter of this book. 
The reviewer would prefer a treatment which highlights, 
to a greater extent, the general ideas and methods of the 
subject. The student of this book will learn a good deal of 


mechanics, but he may not acquire an adequate apprecia- 
tion of the structure, power, and beauty of this science. 
MartTIN J. KLEIN 
Case Institute of Technology 


The Nature of Some of Our Physical Concepts. P. W. 


BRIDGMAN. Pp. 64. Philosophical Library, New York, 
1952. Price $2.75. 


The reconstruction of physical concepts which has oc- 
curred through theories of relativity and quanta has 
furthered the establishment of philosophy of science as a 
special field of knowledge. In recognition thereof a Depart- 
ment of the History and Philosophy of Science was created 
at City College, University of London, under the leader- 
ship of Professor Herbert Dingle; recently there was 
founded the British Journal for the Philosophy of Science. 
In April, 1950, Professor Bridgman, who has contributed 
significantly to the analysis of concepts and methods of 
science, delivered three lectures in London at the invitation 
of Professor Dingle. These lectures, which were then 
published in the just cited British Journal, have now 
appeared as the book which is the subject of the present 
review. 

In the Logic of Modern Physics, 1927, Professor Bridg- 
man expounded a theory of concepts which under the name 
of Operationism has risen to the status of a philosophical 
doctrine. In the present book he endeavors to develop an 
awareness of the operational content of several of the 
important concepts of physics. In the first lecture he states 
that the fundamental idea back of an operational analysis 
is simple enough, namely, that we do not know the meaning 
of a concept unless we can specify the operations which 
were used by us or our neighbor in applying the concept 
in a concrete situation. The first lecture also discusses the 
concepts of field and empty space with application to the 
propagation of light. In the second lecture he examines the 
concepts of internal energy and entropy in thermody- 
namics, and especially he finds need for flow vectors and a 
creation function for entropy. In the third lecture the 
author examines concepts for the description of electrical 
and thermal phenomena in massive conductors. 

Professor Bridgman’s original statement of his opera- 
tional theory appeared inadequate to specialists in logical 
analysis. The theory was not formulated in terms of stand- 
ard distinctions for modes of meaning, such as intension 
and extension, sense and denotation: Further, inasmuch as 
rationalist philosophers had stressed the creative element 
of thought in the constitution of science, Professor Bridg- 
man’s apparent emphasis on physical operations seemed to 
minimize the rational constituent of science. This apparent 
emphasis on physical operations encouraged psychologists 
to interpret mental states in terms of modes of behavior. 
From another point of view, Professor Bridgman’s empha- 
tic recourse to immediate experience appeared to imply 
acceptance of phenomenalism in metaphysics. 

In response to discussion of his work, Professor Bridg- 
man subsequently declared that he recognized mental as 
well as physical operations to be constitutive of meaning. 
In the present work he stresses mental operations; among 
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such he singles out for special attention mathematical 
manipulations which he characterizes as ‘‘paper-and-pen- 
cil’ operations. He further cites verbal operations and 
emphasizes the function of verbal demands in the develop- 
ment of science. In the present book there is less emphasis 
on immediate experience, but the language is weighted in 
favor of empiricism as represented by nominalism in the 
theory of concepts and behaviorism in psychology. 

In his discussion of the concepts of field and empty 
space, Professor Bridgman argues that there is no instru- 
mental operation by which action at a distance can be 
distinguished from action through a field. In both cases the 
propagation of influence is to be detected by a test body. 
He recognizes, however, that propagation of light by action 
at a distance would require delay in time. In the reviewer’s 
opinion, such delayed action in empty space would not 
meet the requirements of the principle of sufficient reason, 
which principle is satisfied by the determination of a finite 
velocity by the properties of an intervening medium. 

Physical properties essentially are manifested in inter- 
actions as dispositional properties which stand in certain 
relations. Thus, a physical concept signifies a property that 
is constituted by relations which are expressed by laws. 
Operations may be required to set the stage for the dis- 
covery of a property, but operational procedures merely 
enable one to isolate a property so that its relation to some 
other property may be discerned. But this analysis, which 
was adequately expressed by Helmholtz in Zaehlen und 
Messen, erkenntnistheoretisch betrachtet, 1887, does not con- 
tradict operational analysis, which Professor Bridgman 
characterizes as a necessary but not sufficient condition of 
meaning. Indeed, he recognizes that the first and second 
laws of thermodynamics furnish definition of internal 
energy and entropy. He finds that by instrumental supple- 
mentation a definition develops into a natural law; the 
reviewer, however, would argue that initially the transition 
occurs in the reverse direction. Regardless of this question, 
acknowledgement that functional relations may serve as 
definitions of concepts implies that functional relation con- 
stitutes meaning. 

V. F. LENZEN 


University of California, Berkeley 


Symmetry. HERMANN WEYL. Pp. 167, Figs. 72, 14.5 22.5 
cm. Princeton University Press, Princeton, New Jer- 
sey, 1952. Price $3.75. 


Symmetry is not only basic in the forms which man 
creates but it is one of the ideas by which he tries to under- 
stand the order and relations of nature. This little book 
aims to show the wide application of this idea as it cuts 
across various fields and it further aims, after making the 
definition precise, to outline its mathematical essence. 

Any writing by so distinguished a mathematician as 
Professor Wey! must deserve attention, and here is a semi- 
popular summary of the work of a penetrating mathemati- 
cal mind in viewing the world through the one small but 
fascinating window of symmetry. The volume has been 
described by Weyl himself as a “‘swan song on the eve of 
my retirement from the Institute of Advanced Study.” 


Is there an intrinsic difference between right and left? Can 
a uniform hexagonal network be made to cover a sphere? 
Are all a priori statements in physics based on symmetry? 
Weyl answers these questions and many others. 

The book is richly illustrated with many striking ex- 
amples from the arts and from organic and inorganic 
nature. These range from the cells of a beehive to living 
cell division, from the geometric intricacies of Arabian 
ornament to Laue spots, from the mosaics of Ravenna 
and from Sumerian and Babylonian design to the triple 
pole symmetry of the iris and the hexagonal symmetry of 
the snow flake. 

The four chapters are a slight revision of the Louis Clark 
Vanuxem Lectures given at Princeton University in 1951. 
The first two lectures discuss the subjects of (1) Bilateral 
Symmetry, and (2) Translational, Rotational and Related 
Symmetries with illustrations chiefly from the arts and 
biology. The third lecture on Ornamental Symmetry begins 
a somewhat more systematic attempt to analyze the parti- 
cular kind of geometric symmetry which in two dimensions 
leads to the seventeen possible types of symmetry involved 
in ornament. In the last lecture this line of thought is 
extended to the three-dimensional case basic to crystallo- 
graphic analysis. The expressed aim in this lecture is to 
exhibit the principle of symmetry in both physics and 
mathematics and to arrive at a final general statement 
(mathematical) of that principle. 

The theories of relativity and quanta are referred to in 
the final lecture, the former particularly for the symmetry 
of the transformations of the Lorentz group, and the latter 
through brief reference to the kinds of symmetry involved 
in atomic structure and spectra. 

There is much for the nonmathematical reader who may 
also glimpse the workings of the mathematical mind. De- 
tailed steps are omitted, but the final mathematical con- 
clusions would be clearer to those with some knowledge of 
modern geometry. Many readers will benefit more from 
seeing how the principle applies in fields other than their 
own. There are two appendices containing mathematical 
proofs. 

Rocers D. Rusk 
Mount Holyoke College 


New Members of the Association 


The following persons have been made members or junior members 
(J) of the American Association of Physics Teachers since the publica- 
tion of the preceding list [Am. J. Phys. 20, 598 (1952)]. 


Alexander, William Everett (J), R.R. 6, Marion, Ind. 

Bade, William Lemoine, Brace Laboratory, University of 
Nebraska, Lincoln, Neb. 

Basye, Clifton A., 103 Bellevue Dr., Richmond, Ky. 

Bounds, Nell M., Midwestern University, Wichita Falls, 
Texas 

Bradley, George Edgar, Western Michigan College, Kala- 
mazoo, Mich. 

Dittrich, William J., 515 Lakeway Dr., Bellingham, Wash. 

Dunkerley, Herbert Spencer (J), 102 N. Coler St., Urbana, 
Til, 
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Hart, Edward Hutson, 1055 Pine St., Menlo Park, Calif. 

Hawley, Fred T., 187 Hamilton, Muskegon, Mich. 

Jewell, Paul Vernon, Maryland State College, Princess 
Anne, Md. 

Johnson, Carl Valdemar, 285 Main St., Portland, Conn. 

Johnson, Halbert Herring, 7346 Arnim, Houston 17, 
Texas =, 

Kriessman, Charles Joseph, 1411 Viers Mill Rd., Rock- 
ville, Md. 

Leonard, Robert S. (J), P.O. Box 9342, University Station, 
Reno, Nev. 

Marion, Jerry Bascervil (J), Physics Department, Rice 
Institute, Houston 1, Texas 

Morley, Arthur L., Rt. 5, Hendersonville, N.C. 

Pollack, Edward (J), 681 E. 181 St., New York 57, N. Y. 


AMERICAN JOURNAL OF PHYSICS 


VOLUME 21, 


Rapp, Thomas Louis (J), Apt. B-8-B College View, College 
Station, Texas 

Schrader, Carlton Davis, Physics Department, State 
University of Iowa, Iowa City, Iowa 

Schriever, William W., Jr., 910 Chautauqua, Norman, 
Okla. 

Stecher, Milton, 916 Union St., Brooklyn 15, N. Y. 

Stephenson, Harold P., Illinois Wesleyan University, 
Bloomington, III. 

Stubbs, George R., 3126 23rd Ave., S., Minneapolis 7, 
Minn. 

Vernon, Robert Carey, Department of Physics, Williams 
College, Williamstown, Mass. 

White, James Frederick, 1949 Woodard St., 
Texas. 


Abilene, 


NUMBER 1 JANUARY, 1953 


RECENT MEETINGS 


Western Pennsylvania Section 


The Western Pennsylvania Section of the American 
Association of Physics Teachers met at Carnegie Institute 
of Technology on December 1, 1951. Forty-three members 
and guests were present. PROFESSOR CHARLES WILLIAM- 
son, Carnegie Institute of Technology, presided. Dr. J. C. 
WARNER, president of Carnegie Tech, welcomed the mem- 
bers and presented an invited paper entitled The present 
status of the scientific manpower problem. The following 
contributed papers were presented: 


Discussion of the report of the Coulomb’s Law Com- 
mittee of the AAPT. W. H. MICHENER, Allegheny College. 
The history of the committee and the background of the 
report were given. Some of the important recommenda- 
tions of the committee were explained. The speaker be- 
lieves that, in the teaching of magnetism, the traditional 
model of magnet poles will gradually be supplanted by 
something like the Amperian current theory. A coming 
generation of physicists may be as unfamiliar with Cou- 
lomb’s law for magnets as our present younger generation 
is with stresses in the ether. 


A decade scaler using neon tubes. MERLE N. Hirsu, 
University of Pittsburgh. A decade scaler using neon tubes 
is now in development as an undergraduate thesis project. 
Based on a circuit by Manley and Buckley,! it offers 
economy of power consumption and panel space, dependa- 
bility at high counting speeds, and relative insensitivity 
to variation in supply voltage. To increase dependability, 
but not initial cost, we have replaced the germanium diodes 
of the original circuit with 6ALS double diodes, one tube 
to two counting stages. Eventual counting speeds of 1500 
cps should be possible. 


1J. Manley and E. Buckley, Electronics 23, 84 (1950). 


Teaching the connection between the physical world and 
abstract mathematical analysis. S. DE BENEDETTI, Car- 
negie Institute of Technology. The mind of the student 
should be kept focused on what happens and not on what 
is written in textbooks. It is important that, chiefly at the 
beginning (mechanics), problems should be of an experi- 
mental nature and involve some measurement. As a second 
step the student should be trained in the process of abstrac- 
tion which leads from the direct experience with a concrete 
physical situation to the formulation of a ‘‘model” capable 
of mathematical analysis. Only when the student is thor- 
oughly familiar with this process should one introduce 
“metaphysical” entities such as booms 5 ft long (exactly), 


weights of 100 lb (exactly), frictionless planes, and weight- 
less strings. 


Wanted: Generic term for a class of units. Ceci O. 
Riccs, Waynesburg College. There seems to be no name to 
distinguish the correct unit for any quantity from auxiliary 
units of the same quantity. For example, where an equa- 
tion requires or gives a result in erg rather than joule, 
ft-lb per sec rather than horsepower, etc., there is no 
general statement with which the student can be warned. 
The author of one textbook uses the term fundamental 
unit for the proper unit and derived unit for the auxiliary 
unit, but these terms have of course been pre-empted for 
use in another sense. Some writers seem to use the term, 
absolute unit, for the proper unit, but this is ambiguous 
as absolute unit is usually used for any unit of an absolute 


system. Common agreement upon a general term would 
be a convenience. 


Science is fun at the Buhl Planetarium. RicHarp C. 
Hitcucock, Buhl Planetarium. Our primary purpose is 
to show science principles to people of all ages. For young 
audiences, school work is supplemented rather than sup- 
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planted. For adults, the science of everyday living is 
stressed. It is essential that exhibits be safe. In addition, 
action exhibits must have sturdy mechanical construction 
to permit continued operation. This latter is most im- 
portant, as a museum may be regarded as a laboratory run 
by inexperienced students. All-day, every-day operation 
is a real problem for the exhibit designer and the main- 
tenance staff. 


The Muller microscope. BERNARD L. BRINKER, St. 
Vincent College, Latrobe, Pennsylvania. With blackboard 
sketches and verbal descriptions, the Muller microscope 
was reported to the meeting. This new instrument is capable 
of very high resolution (about 0.3 millimicrons) which 
allows the imaging of details in the crystal structure of 
metals down to about atomic diameters. The microscope 
was devised by Dr. Erwin Muller of the Kaiser Wilhelm 
Institute at Berlin, Germany. It was demonstrated by 
Dr. Muller at the Ninth Annual Pittsburgh Conference on 
X-ray and Electron Diffraction (November 30, 1951). 
With the instrument and a projection lens, the structural 
patterns of tungsten and barium crystals and the crystal- 
lizing process itself were shown at the conference. 


The “cycle” in a physical interpretation of Planck’s 
constant (h) for circular dynamics. BERNARD L. BRINKER, 
St. Vincent College, Latrobe, Pennsylvania. In an analysis 
of periodic movement, the cycle or its equivalent has 
physical significance. It cannot be omitted from an expres- 
sion without changing the physical concept explicitly 
stated. By an explicit transformation of the physical quan- 
tities involved in Planck’s constant (h=6.62 10-27 erg 
sec/cycle), an equivalent torque-impulse (#=1.05 1077 
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cm dyne sec) can readily be derived for circular (orbital) 
dynamics. Such a conceptual transformation had been 
implied by Dirac’s introduction of h=h/2zx for quantum 
mechanics generally. 

At the business meeting in the afternoon, the group 
voted to change the name of the organization to ‘“‘Western 
Pennsylvania Section of the American Association of 
Physics Teachers.’’ New officers are President, W. C. 
KELLY, University of Pittsburgh; Vice-president, R. L. 
Brown, Allegheny College; Secretary, R. C. HitcHcock, 
Buhl Planetarium. The secretary read the reply of the 
Pennsylvania State Council of Education to the recom- 
mendation of the Section that certification requirements 
for high school physics teachers include at least 6 semester 
hours of college physics. The Council feels that the present 
supply situation makes it unwise to take action that would 
result in more difficulty for school districts to obtain 
teachers and therefore withheld approval of the recom- 
mendation. The Section voted to continue an interest in 
this matter by including an item in the minutes of future 
meetings as a reminder for future action. 

FATHER BERNARD L. BRINKER, the national representa- 
tive of the Section, reported on proceedings of the business 
section of AAPT at its last New York meeting. He also 
called attention to the annual meetirgs of the Pennsylvania 
Physics Conference. The Section instructed the secretary 
to send a letter of sympathy to the family of Dr. W. N. 
St. Peter. A vote of thanks was given to the Carnegie 
Tech members for their hospitality. The spring meeting 
will be held at Allegheny College on May 10, 1952, and 
the fall meeting at Buhl Planetarium in December. 


W. C. KELLy, 
Secretary 


How to fecure Houfes, &c. from LIGHTNING 


It has pleafed God in his Goodnefs to Mankind, at length to difcover to them the Means of fecuring 
their Habitations and other Buildings from Mifchief by Thunder and Lightning. The Method is 
this: Provide a [mall Iron Rod (it may be made of the Rod-iron ufed by the Nailers) but of fuch a 
Length, that one End being three or four Feet in the moift Ground, the other may be fix or eight 
Feet above the higheft Part of the Building. To the upper End of the Rod faften about a Foot of 
Brafs Wire, the Size of a common Knitting-needle, fharpened to a fine Point; the Rod may be 
fecured to the Houfe by a few {mall Staples. If the Houfe or Barn be long, there may be a Rod and 
Point at each End, and a middling Wire along the Ridge from one to the other. A Houfe thus fur- 
nifhed will not be damaged by Lightning, it being attracted by the Points, and Paffing thro the 
Metal into the Ground without hurting any Thing. Veffels alfo, having a [harp pointed Rod fix'd 
on the Top of their Mafts, with a Wire from the Foot of the Rod reaching down, round one of the 
Shrouds, to the Water, will not be hurt by Lightning —Poor Richard’s Almanac, 1753. 
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